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NUMERICAL ANALYSIS OF THE ANDERSON LOCALIZATION 

P. Markos' 
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The aim of this paper is to demonstrate, by simple numerical simulations, the main trans- 
port properties of disordered electron systems. These systems undergo the metal insulator 
transition when either Fermi energy crosses the mobility edge or the strength of the disorder 
increases over critical value. We study how disorder affects the energy spectrum and spatial 
distribution of electronic eigenstates in the diffusive and insulating regime, as well as in the 
critical region of the metal-insulator transition. Then, we introduce the transfer matrix and 
conductance, and we discuss how the quantum character of the electron propagation influ- 
ences the transport properties of disordered samples. In the weakly disordered systems, the 
weak localization and anti-localization as well as the universal conductance fluctuation are 
numerically simulated and discussed. The localization in the one dimensional system is de- 
scribed and interpreted as a purely quantum effect. Statistical properties of the conductance 
in the critical and localized regimes are demonstrated. Special attention is given to the numer- 
ical study of the transport properties of the critical regime and to the numerical verification of 
the single parameter scaling theory of localization. Numerical data for the critical exponent 
in the orthogonal models in dimension 2 < d < 5 are compared with theoretical predictions. 
We argue that the discrepancy between the theory and numerical data is due to the absence of 
the self-averaging of transmission quantities. This complicates the analytical analysis of the 
disordered systems. Finally, theoretical methods of description of weakly disordered systems 
are explained and their possible generalization to the localized regime is discussed. Since 
we concentrate on the one-electron propagation at zero temperature, no effects of electron- 
electron interaction and incoherent scattering are discussed in the paper. 



1 Introduction 

Localization of electrons in disordered system has fascinated scientists for almost fifty years. In 
1958, Anderson [1] predicted that randomly distributed impurities in the crystal lattice can local- 
ize an electron in a certain spatial region. The localization is given by the quantum character of 
electron propagation: electron wave function is scattered on randomly distributed impurities, and 
mutual interference of scattered components cancel the wave function on large distances. Local- 
ization of the electron is responsible for a new kind of insulators - the Anderson insulator, which 
possesses zero electric conductivity, cr, in a part of the energy bands, where the density of states, 
p{E), is non-zero. Propagation of the electron in disordered systems is therefore qualitatively 
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different from that in periodic structures. Although the electron wave function is reflected also by 
a periodic potential, the interference of the reflected and the transmitted waves in a regular lattice 
gives rise to bands and gaps in the electron energy spectrum. In the bands, where the density of 
states is non-zero, the electron propagates freely throughout the structure. 

Of course, disorder is always present in the real world, and it influences the electric trans- 
port. Scattering on weak impurities causes diffusive propagation of electrons. The electrical 
conductivity can be expressed through the diffusive coefficient, D{E), and the density of states, 
PiEl 

a = e^D[E)p{E). (1) 

In the derivation of the expression Q, it is assumed that an electron on its travel through the 
sample is scattered on individual impurities. This assumption is correct only if the electron 
wavelength, Xp, which is determined by the wave vector, /cf on the Fermi energy, is smaller 
than the electron mean free path due to the coherent scattering on impurities, £, 

^ » 1. (2) 

[2]. The mean free path represents, in the first approximation, the mean distance between two 
impurities. Clearly, I is large in the limit of weak disorder, and decreases when disorder in- 
creases. Therefore, condition (|2j is violated for strong disorder and localization is expected 
when i ^ Xp- The last relation is known as the Lifshitz criterion for localization. 

Increase of disorder changes the transport regime considerably. In the limit of very strong dis- 
order, all electronic states become localized. The wave function of localized electron decreases 
exponentially as a function of distance from the localization center, rg, 

*(f)^exp-^^, (3) 

where A is the localization length. 

Anderson showed in his pioneering work [1] that all electronic states become insulating when 
disorder increases above the critical value. The transition from metal to insulator due to the 
increase in disorder is called the Anderson transition. Similar to the theory of phase transitions, 
it is believed that Anderson transition is universal and can be described by the one-parameter 
scaling theory [3]. The key parameter of the scaling theory is the conductance, g. Introduced by 
Landauer [4], the conductance measures the transmission properties of disordered systems both 
in the metallic and the localized regime. 

The scaling theory of localization analyzes the size dependence of the conductance in the 
limit of large system size. Only three transport regimes exist in this limit: the system is either 
in the metallic, localized or critical regime. In the metallic regime, the conductance increases 
with the system size, and the system possesses non-zero electric conductance. In the localized 
regime, all electronic states are localized, and g decreases exponentially due to further increase 
of the system size. The system is in the critical regime only at the critical value of the disorder. 

While the localized regime exists in any system, provided the disorder is sufficiently strong, 
the existence of the metallic regime is not guaranteed, especially for systems with lower spatial 
dimension. It is well-known that all electronic states are localized in any one dimensional (ID) 
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disordered system [5]. Spin-less electrons are always localized already in the two dimensional 
(2D) systems, with the only exception caused by an external magnetic field. 

The scaling theory of localization predicts that the metal-insulator transition is universal. 
The size and disorder dependence of the conductance in the neighborhood of the critical point 
is governed only by universal critical exponents. All the parameters, which define the micro- 
scopic structure of the model, become irrelevant when the size of the system is sufficiently large. 
Verification of the universality of the metal-insulator transition and the calculation of the critical 
parameters - the critical disorder and critical exponents - for systems of various dimension and 
physical symmetry are the main theoretical and numerical problems of the theory of localization. 

Contrary to the scaling theory of localization which discusses the transport regimes in the 
limit of L ^ oo, in everyday life we must deal with systems of finite size. Here, the transport 
regime depends on the relation of the system size, L, to the characteristic lengths. For instance, 
if 

^ < i < A, (4) 

then the system exhibits metallic behavior with the conductivity given by Eq. Q. This happens 
in 2D systems, where the localization length A is extremely huge for weak disorder [6]. Since 
the electron diffuses through the sample, we called the transport regime, defined by inequalities 
(|4j, the diffusive regime. Of course, an increase of the system size over the localization length 
causes localization of all electronic states and the conductivity decreases to zero. 

There are small quantum corrections to the conductivity Q when the conditions (|4} are 
fulfilled. For instance, if the size of the two dimensional (2D) sample is much larger than the 
mean free path, L ^ £, the mean value of the conductance decreases logarithmically when the 
size of the system increases. This effect - weak localization corrections to the conductance - 
is the first manifestation of the quantum character of the electron propagation. Similar weak 
localization corrections to the conductivity exist also in one and three dimensions and in quasi- 
one dimensional (quasi- Id) geometry. 

In the opposite limit of very small systems, L £ we observe the ballistic regime. In 
this regime, the electron, in its travel through the sample, is scattered only on a few impurities. 
Clearly, the sample is almost transparent and the conductance might be large akeady in ID 
systems. 

While the weakly disordered systems can be described analytically, for instance by the 
Dorokhov-Mello-Pereyra-Kumar (DMPK) equation [7], or perturbation Green's function meth- 
ods [8,9], the theoretical description of the critical regime is still not complete. The main problem 
is the absence of small parameter, since the critical disorder is of the order or even larger than the 
bandwidth in three dimensional systems. The critical disorder is small only when the dimension 
of the system, d, decreases to the lower critical dimension, dc = 2 [10]. Analytical theories in 
dimension d — 2 + e [10-12] calculate critical parameters of the model in powers of e. 

Since each sample contains randomly distributed scatterers, the measured quantities, like the 
conductance g, fluctuate from sample to sample and must be averaged over random disorder. 
One possibility is to average the conductance over the statistical ensemble of macroscopically 
equivalent samples, which differ only in the microscopic realization of the disorder The ergodic 
hypothesis [13] states that the same statistical ensemble can be constructed with the use of the 
single sample by varying the Fermi energy or the magnetic field. This enables us to compare 
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the theoretical or numerical results, which use the ensemble averaging, with experiments, where 
usually only a few samples are analyzed. 

The key parameter of the scaling theory, the conductance, g, is not a self-averaged quantity. 
Already in the metallic regime, the conductance becomes sample-dependent. Measured values 
of the conductance [14] fluctuate around the mean value. These fluctuations, of the order of 
e^/h, depend neither on the size of the system, nor on the mean value [9] and they lead to the 
absence of self-averaging of the conductance. In the localized regime, the fluctuations of the 
conductance [15] are so strong that the mean value, (g), is not a representative quantity [16]. 

The absence of self-averaging and huge fluctuations of the conductance which do not disap- 
pear in the limit of infinite size of the system must be included in the theoretical description of 
the localization. First, it is not clear how the averaging over the disorder should be performed. 
Moreover, not only the mean values, but also the higher cummulants of the conductance must 
be calculated. Contrary to the classical systems, the higher cummulants do not diminish in the 
limit of large system size. Also, the averaged quantity must be carefully chosen. In the localized 
regime, it is more suitable to average the logarithm of the conductance than the conductance 
itself. The average procedure is easier in the numerical simulations than in the analytical theory. 
This is the reason why numerical methods, based on the finite size scaling [17] provides us with 
the most reliable information about the Anderson transition in higher dimensions. 

In this Paper, we discuss the basic ideas of the localization theory. We concentrate on nu- 
merical methods of investigation of disordered electronic systems. Numerical simulations enable 
us to describe the metallic, the critical and the strongly localized regime. They verify theoret- 
ical predictions, and, last but not least, they provide us with data necessary for the analytical 
description of the localization. 

As already mentioned, both the weak localization and the localization are quantum effects, 
caused by the mutual interference of the scattered components of the electron wave function. The 
quantum coherence of the electron propagation is destroyed by the incoherent scattering. Since 
the incoherent mean free path, L^, increases when the temperature decreases [18], we expect 
that the best experimental conditions for the analysis of the localization effects are those in the 
limit of small temperature, when ^ L. In this paper, we restrict our discussion only to the 
limit of zero temperature, when ^ oo. For the case of non-zero temperature, when > L, 
we can obtain a qualitative estimation of the transport in the diffusive regime if the size of the 
system is replaced by L^. The transport in the strongly localized regime, > X requires more 
detailed analysis [19,20] which is above the scope of this paper. 

We will consider only the one-electron problem. Although the electron-electron interaction 
might play important role in the localization [21], numerical analysis of systems of disordered 
interacting electrons is too difficult to be able at present to answer the main questions of the 
scaling theory. 

Since the single electron localization is caused by interference of electron wave function, 
the theory can be easily applied also to propagation of classical waves in disordered media. Of 
special interest is propagation of electromagnetic waves through random dielectrics [22-24]. 

The Paper is organized as follows. 

In Sect. |2]we describe the Anderson localization and demonstrate the localization of electron 
in two dimensional (2D) strongly disordered system. In Sect. |3]we introduce several models, 
used in numerical simulations. For numerical convenience, all the models are based on the 



Numerical Analysis of the Anderson Localization 



5 



propagation of electrons on the d - dimensional lattice, with disorder represented by the random 
on site energies. The Anderson transition is explained briefly in Sect. |3 Spatial structure of the 
electron wave functions and spectra of eigenenergies is shown in Sect. [S] Transfer matrix and 
the conductance are introduced in Sect. |6l 

In the most simple case of the one dimensional (ID) disordered chain, it was proved [5] that 
all electronic states are localized even for weak disorder provided that the system is sufficiently 
long (longer than localization length. A). We discuss the transport in the ID systems in Sect. 
The one dimensional systems are not only easy to simulate numerically, but allow also to 
derive exact analytical results. The statistics of the conductance and the resistance is discussed 
in details, and the quantum origin of the localization is explained. 

Section |8] analyzes the electronic transport in the diffusive regime. The scaling theory of 
localization, formulated in terms of the conductance, is introduced in Sect. |9l It is argued that 
the metal - insulator transition is a universal phenomenon, and the conductance, g, is a function of 
only one parameter in the critical regime. This scaling hypothesis is verified numerically. First, 
in Sect. ^Jand^Jwe discuss statistical properties of the conductance in the critical and localized 
regime. The critical conductance distribution is presented and discussed in detail. Then, in Sect. 
[T2I we review the numerical scaling analysis of the Anderson transition in the Anderson model. 
Sect. Elpresents the critical exponent, obtained by numerical scaling analysis of systems with 
dimension 2 < c? < 5. Obtained numerical data are compared with theoretical predictions. 
The two dimensional critical regime is discussed in Sect. \^ Finally, Sect. discusses two 
possibilities of the analytical description of the localized regime. 

The Paper contains five Appendices, which present some technical details useful for the 
understanding of discussed numerical data. Properties of the transfer matrix are reviewed in 
Appendix IaI The next two Appendices introduce two successful theories of the transport in 
diffusive regime: the DMPK equation in App. Oand random matrix theory in App. |Cl Lyapunov 
exponents of the product of transfer matrices are introduced in AppendixIDl The last Appendix 
|E|discusses numerical algorithms for calculation of the conductance. 

Various aspects of the electronic transport in disordered systems were reviewed recently. 
The weak localization effects are discussed in [25]. Refs. [26-28] review the main ideas of 
the localization theory. Experimental and theoretical aspects of the locahzation are reviewed 
in [30]. The Quantum Hall effect is discussed in refs. [31, 32]. Supersymmetric field theory and 
its application to electron transport is explained in [33]. Statistical properties of spectra and wave 
functions are reviewed in [34]. Random matrix theory and its application to electronic transport 
is reviewed in [35]. The development of research is in conference proceedings [36-38]. 

Wave transport in disordered media is a subject of many textbooks and monographs. We 
want to mention the classical book of Economou [39] and Mott and Davis [40]. Transport of 
electrons in mesoscopic systems is discussed in [41^3]. The book of Mehta [44] presents the 
theory of random matrices. 

2 Localization 

Localization of an electron in disordered system was predicted by Anderson [1]. Anderson cal- 
culated the probability p that an electron, being at time t = at point fo with the wave function 
^'(r, i ~ 0) = (5(r — ro) returns to the same point in time t > 0. It is evident that p = Q for 
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Fig. 1. The time dependence of the wave packet in the two dimensional disordered system defined by 
Hamiltonian ^6) with Box disorder and = 6. This disorder corresponds to the localization length A ~ 50 
[45]. Different colors distinguish points where the wave function > 0.0001 (gray), [^^(r)! > 0.0010 

(dark gray), and > 0.0050 (black) at different times shown in legend. The time is measured in units 

of lOOh/V with hopping term V = 1. The size of the lattice is 512 x 512 (in units of lattice spacing a). The 
initial wave function, ^{f,t = 0), is the eigenfunction of the sub-lattice of the size A*';i x Nh {Nh ~ 24) 
of the system, located in the middle of the lattice. The eigenfunction belongs to the eigenenergy closest to 
the band center, E = 0. From the time development of the wave packet, one can see that the electron is 
indeed localized in the center of the lattice. 



regular lattice, since the electron propagates freely, provided that its energy lies in the allowed 
energy band The question is whether there are such disordered systems, where p ^ 1 in the limit 
t —^ oo. If yes, then the electron is localized in the space around tq. If not, then electron can 
propagate through the sample, in spite of the disorder. 

To demonstrate Anderson's idea, we simulated numerically the time evolution of a single 
electron wave function in two-dimensional disordered lattice. At time < 0, we add the electron 
to the center of the lattice and solve the time-dependent Schrodinger equation, 

^h^^^^^H^if,t) (5) 
ot 

where H is the Anderson Hamiltonian, 

n^Wj2^rclcr + Vj24cr'- (6) 
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Fig. 2. The wave function of the electron in time t — 2800 h/V. The size of the lattice is 512 x 512 (in 
units of lattice spacing, a). Note that the radius of occupied region is much larger than the localized length, 
which is in this case A « 50. The localization length determines only the exponential decrease of the wave 
function at large distances, not the size of the region in which the electron is localized. 



Following Eq. (|6jl, the electron moves on the d-dimensional lattice of size L"^. The first sum is a 
local term, where r counts sites on the lattice and the disorder is given by random on-site energies 
Er- The second term in Hamiltonian enables the hopping of electron between the nearest neighbor 
sites. Parameter V is given by the overlap of the wave functions localized on neighboring sites, 
and W defines the strength of the disorder. The distance between two neighboring sites, a, is 
used as a unit length, a = 1. 

Figure^shows the time development of the electron wave function in the strongly disordered 
two dimensional (2D) lattice. At the beginning, the electron wave function broadens with time. 
However, later this broadening becomes slower and finally stops. The electron is localized in the 
central part of the lattice and its wave function far away from the localization center is negligibly 
small. Figure|2]shows the detailed spatial distribution of the wave function at time t = 2800y/?i. 

The presented data show that indeed the electron might be localized for sufficiently strong 
disorder. Anderson calculated the critical strength of the disorder, Wc, such that the electron is 
delocaUzed (and the system is a metal) when W < Wc, and the electron is localized and the 
system becomes an insulator when W > Wc- The transition of the system from the metallic to 
insulating regime due to increase of the disorder is called metal-insulator transition. 

We will not reproduce Andersons analysis here, only present his main result: there is indeed 
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the critical disorder which separates locaHzed and delocahzed electron states. In the most simple 
approximation, critical disorder can be found as 

Wr 

-y Ki2eK\n{eK) (7) 

where if is a connectivity of the lattice. A detailed discussion of Anderson analysis is given in 
Refs. [46]. 

Note, relation Q does not contain the dimension d of the system. Later [10] it became clear 
that, similarly to critical phenomena in statistical physics, the dimension of the system is crucial 
for the existence of the Anderson transition. The dimension d = 2 is a lower critical dimension; 
there is no metallic state for d < 2. That means, all electron states are localized in space when 
d < 2. However, localization can be observed only when the system is large enough, 

L > A. (8) 

In 2D, the localization length is extremely large for weak disorder [6]. This is the reason why 
good metallic behavior is observed in numerical simulations performed on 2D weakly disordered 
samples of finite size, L. (Sect. |8}. 

Analytical estimation of the critical disorder is possible only when the dimension of the 
system is close to 2, d = 2 + e with e ^ 1. In realistic three dimensional systems, Wc can be 
obtained only numerically. Various methods of calculation of Wc and of other critical parameters 
will be explained in Sect. 



3 Models and symmetries 



The most suitable model for the numerical analysis of the localization is Andersons model, de- 
fined by the Hamiltonian (|6}. It represents an isotropic model, in which the hopping term, V, 
is the same in all directions. In the isotropic models, we consider V = 1. This also defines the 
scale of the energy. 

It is often suitable to consider anisotropic models, with different hopping terms in different 
directions. For instance, we will use the three dimensional (3D) model 

Ti. — W ^ ^ Erclcr + t ^ ^ C'xyz^x'yz + t ^ ^ '^'xyz^xy'z + V ^ ^ '^'xyz'^xyz' (9) 
r [rr'] [rr'] [rr'] 

where t < V. The hopping term, V, defines the scale of the energy. We often use V — 1. 
Random energies e„ are distributed either with the Box distribution, 

PB{e)^{2/W)e{W/2-\e\) (10) 

or with the Gaussian distribution, 

PG(e) = (27rT4^^)exp-(eV2W^^). (11) 

The parameter W measures the strength of the disorder. Other distributions of random energies 
(binary, Lorentzian) are often used in the literature too. 
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The Anderson model defined by Eq. (|6j will be used in the present paper to demonstrate 
numerically the basic ideas of the localization theory. For zero disorder, Sr = 0, we easily find 
all the eigenenergies of Hamiltonian (|6jl. For instance, for the 3D anisotropic model we have 

E = 2tcoska; + 2tcosky + 2Vcosk^. (12) 

The energy band spans between — V — 2t and V + 2t. We define the bandwidth, 

B = 2V + 4:t. (13) 

Andersons model, given by Hamiltonian (|6j, belongs to models with time reversal symmetry. 
It describes the propagation of a single spin-less particle in the disordered lattice. Such models 
are called orthogonal. Systems with orthogonal symmetry do not exhibit the metal-insulator 
transition in dimension dc = 2. This changes when the hopping between neighboring sites 
becomes dependent on the orientation of the spin of the electron [10]. Evangelou and Ziman [47] , 
and Ando [48] showed numerically that the two dimensional (2D) systems with spin dependent 
hopping (called symplectic models) exhibit the metal-insulator transition already at dc = 2. In 
this paper we will discuss the Ando model which is defined by the Hamiltonian 

r [rr'] [rr'] 

with spin dependent hopping terms given by Ando [48], 

y^-[ Z and ), (15) 

with + — 1. Note how hopping depends on the direction of propagation. Also, the wave 
function, c^, has two components, 

) • 

for different orientations of the spin of electron. Another 2D model with symplectic symmetry 
is the Evangelou-Ziman model. In this model, the hopping term has the form 

y^( '^ + ivtz -V{ty-it^)\ 

\ V{ty + it^) 1 - IVt^ ^ ' 

where t^, ty and are real random variables chosen from box distribution, and v measures 
the strength of (random) spin-orbit hopping. In Refs. [49], SU(2) model is analyzed, with the 
random hopping term 

V=(Cr^ -^r^iV (18) 
\ e ^ sm /3 e cos /3 y 

where random phases a and 7 are uniformly distributed in the interval (0, 27r) and the distribution 
of random phase f3 is p{f3)df3 = sin(2/3)d/3 for < /3 < n/2. 
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An external magnetic field breaks the time reversal inversion of the electron propagation. 
Such systems belong to the unitary universality class. Magnetic field can be added to the two 
dimensional Hamiltonian (|6} by Peierls factor, 

^ A^xy ^ '^x+a^y'-^xy ^ 2^^y c '^x-a,y'^xy 

where = Bieaf' jfi is a magnetic flux through the elementary plaquette of the size . 

Hamiltonian is used for the numerical analysis of the critical quantum Hall regime of 
2D disordered system in a strong magnetic field. Although this model does not exhibit metal- 
insulator transition, there is, in each Landau band, the critical energy, E^, at which the electron 
is delocalized [31]. The existence of the critical - delocalized - state inside the energy band is 
crucial for the transmission of the system between two neighboring quantum Hall plateaus [50]. 
Another model, often used in the numerical analysis of the critical quantum Hall regime is the 
model of Chalker and Coddington [32,51]. 



4 Anderson transition 



The Anderson transition is the transition from the metallic (extended) regime to the localized 
regime. Before describing the scenario of the Anderson transition, we must understand how the 
disorder influences the density of electronic states in the conducting band. Then, we introduce 
the mobility edge, describe qualitatively the transition from metal to insulator, and define critical 
exponents. 



4.1 The density of states. 

For zero disorder, W — ^, the energy of the electron in the Anderson model is given by the 
dispersion relation 

= 2^008 fc^ + 2T/ cos fcy + 21/ cos fc^, (20) 
where fc = (/c^;, ky,kz) is the wave vector of electron. The density of states, 

p^E) = —1^. (21) 
' 2TTdE ' 

can be calculated analytically [39], and is shown in Fig. |3] The energy band spans from —6V < 

E < +6V and is symmetric, p{E) = p{—E). This symmetry is typical for tight-binding 

Hamiltonians given by Eq. (|6j. 

Note that the density of states, given by Eq. \2\l does not include the degeneracy of the 
electron system due to the two possible orientations of the spin. When spin of the electron is 
included, p{E) must be multiplied by a factor of 2. 

Figure|3]demonstrates the effect of disorder on the energy spectrum of the three dimensional 
Anderson model for Box and Gaussian disorder. As expected, the band becomes broader when 
disorder increases. This is more pronounced for the Gaussian disorder which allows larger fluc- 
tuations of random energies. Open circles indicate the position of the mobiUty edge, discussed 
in the next Section. 

For completeness, we also show in Fig. |4]the density of states of the disordered two dimen- 
sional Anderson model. 
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Fig. 3. The density of states, p{E), of the three-dimensional Anderson model with box (left) and Gauss 
(right) distribution of random energies e„, defined by Eqs. <10t and <11> . respectively. As expected, the 
energy band becomes broader when disorder increases. In the tails of the band, the electronic states become 
localized (Fig. |5}. Open circles show the position of three critical points, one for i5 = and W = Wc, 
and two other for fixed disorder W < Wc and the mobility edge Ed = —Ec2- Note, for weak disorder 
(W = 10 in the left figure) the mobility edges lie outside the unperturbed energy band, \Ec\ > B. 



PiE) 




Fig. 4. The density of states for the disordered two dimensional Anderson model. Note the typical singular- 
ity of the density of states at the band center for the unperturbed tight binding Hamiltonian (W = 0) [39]. 



4.2 Mobility edge 

Intuitively, we expect that the localized states appear first in tails of the energy band. These states 
are namely created by large random energies, Sr- One expects that the states with the energy close 
to the band center are less affected by randomness. Another argument for the locahzation in tails 
was given by Lifshitz [2]. The classical scattering of electrons on impurities, which led to an 
expression for the conductance, given by Eq. Q, is only possible when the mean free path, £, 
is much larger than the wavelength of the electron, Xp, given as the inverse of the Fermi wave 
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Fig. 5. The mobility edge, Ec, separates the localized states in the band tail from the conducting states. By 
changing the Fermi energy, Ef, the system exhibits a transition from the metallic regime (Ef > Ec) into 
the localized regime (Ef < Ec). 

vector, kp. When both lengths become comparable, 

Xf ^ £, (22) 

then the spatial extent of the electron wave function is comparable with the typical distance be- 
tween two impurities, so that we cannot discuss the propagation of electron in terms of individual 
scattering procedures. The electron wavelength, Xp, is small, comparable with the lattice dis- 
tance, a, at the band center, but increases with the distance of the Fermi energy from the band 
center. Consequently, the Lifshitz criterion for the localization, given by Eq. (I22> is fulfilled first 
in the band tails. 

Separation of the energy spectra into localized and delocalized intervals is schematically 
shown in Fig. [S] For a given strength of the disorder, W, there is an energy, Ec, called the 
mobility edge, which separates the localized states from the delocalized states. When Ep < Ec, 
the system is an insulator, since all states at the Fermi energy are localized (we remind the reader 
that the temperature T — 0). We call this insulator the Anderson insulator to emphasize the 
fact, that the system possesses the zero electric conductance, although the density of states is 
non-zero. When the Fermi energy crosses the mobility edge, the system undergoes the transition 
from insulator to metal. For Ep > Ec, the system possesses a finite electric conductivity. 

Note that there are two mobility edges, Ed and Ec2 which separate localized states in the 
upper and lower band tails from the delocalized states in the central part of the band. For a partic- 
ular disorder, we show the mobility edges in Fig. |3] Of course, the position of the mobility edges 
depends on the strength of the disorder. With increasing disorder, both mobility edges move to 
the band center. There is a critical value of the disorder, Wc, for which Ed — i?c2 approach 
the band center. For W > Wc there are no propagating states in the spectrum. The function 
Ec = E{W) defines the critical line, which separates the metallic and the localized states in 
the phase diagram. An interesting feature of the Anderson transition is that the disorder creates 
propagating states outside the unperturbed energy band. We see that for non-zero disorder, the 
metallic regime exists also for energies \E\ > 6V. 

The phase diagram of the three dimensional Anderson model is shown schematically in Fig. 
|6l For W = 0, all the electron states inside the energy band, < 6V are delocalized. In- 
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Fig. 6. Schematic piiase diagram for the 3D Anderson model with the Gaussian disorder. Metallic (conduc- 
tive) states are separated from the localized states by the critical line, which is an envelope of the shadow 
region. Open circles show the position of three critical points, one for E — and W — Wc ~ 6.1, and two 
other for fixed disorder W = 2 < Wc and the mobility edge Ed = —Eci = 6.58. Note, these mobility 
edges lie outside the unperturbed energy band. When disorder reaches its critical value, W = Wc, two 
mobility edges reach the band center, Ed = Ec2- No metallic states exist when W > Wc- 

creasing disorder broadens the density of states, and the interval Ed , Ec2 of the metallic states 
becomes broader, too. Only for rather strong disorder, close to to the critical disorder, Wc, the 
mobility edge starts to converge towards the band center For the disorder W > Wc, all electron 
states are localized. The phase diagram for the 3D Anderson model was calculated numerically 
inRef. [52]. 

4.3 Critical exponents 

The metallic states are separated from the localized ones by the critical line. Any crossing of 
the critical line is accompanied by the metal insulator transition. The common belief, supported 
by the scaling theory of localization [3] is that this transition is universal. This means that the 
transition from metal to insulator should not depend on microscopic details of the model, and on 
the position of the critical point, lying at the critical line. However, it depends on the dimension 
and on the physical symmetry of the model. We can, similarly to the theory of phase transitions, 
formulate the theory of the metal insulator transition in terms of the order parameter and of the 
critical exponents [10]. The critical exponents are defined in terms of the energy or disorder 
dependence of the metallic conductivity and localization length. 

Metallic conductivity, cr, characterizes the transport properties of the metallic regime. We 
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Fig. 7. Definition of thie critical exponents. At the metallic side of the critical point, the conductance a 
decreases to zero as cr oc {Wc — Wy, whereas at the insulating side the localization length diverges as 

X QC {W -Wc)'". 



have non-zero conductivity cr > in metal, and a decreases when the system approaches the 
critical point. For many years it was not clear whether the conductivity at the critical point, <7c, is 
zero or not. Mott [40,53] suggested that the conductance possesses non-zero value at the critical 
point and drops discontinuously to zero on the insulating side of the transition. At present, the 
common belief is that the critical conductivity CTc = in 3D orthogonal systems. The decrease of 
the conductivity in the neighborhood of the critical point is determined by the critical exponent 
s. 



{Wc - wy 



(23) 



On the insulating side of the transition, we can characterize the electronic states by the lo- 
calization length, A. By definition, the localization length which characterizes the exponential 
decrease of the localized wave function far away from the localization center ro, is given by Eq. 



4* (r) cx exp - 



^ol 



A 



(24) 



It is intuitively clear that A should increase when system approaches critical point, and become 
infinite at the critical point. 



X^{W- Wc)- 



(25) 



As discussed in the Sect. 14.11 the Anderson transition is induced either by a change of the 
Fermi energy with fixed disorder strength, or by an increase in the disorder, W , for fixed energy 
(E — 0). Therefore, we expect also that 



a^{Ec~Ey X^iE-E^y 



(26) 



for the case when the Fermi energy crosses the mobility edge. 
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Scaling theory of localization [3] assumes that the metal-insulator transition is universal. 
That means that critical exponents, s and i^, are universal, depending only on the dimension of 
the system, and on the physical universality class. They do not depend on the microscopic details 
of the model. The scaling analysis gives [10] 

s^{d- 2)v. (27) 

Proof of the universality of critical exponents and their dependence on the dimension of the 
system represents the main problem of the theory of Anderson localization. 

5 Wave functions and energy spectrum of disordered systems 

For the case of a regular lattice {W — 0), all the eigenenergies and eigenfunctions of the Hamil- 
tonian can be found analytically. Of course, disorder influences the spectrum of eigenenergies, 
and also the form of the corresponding wave functions. In this Section, we describe briefly the 
main qualitative and quantitative properties of the wave functions and of the energy spectra both 
in the limit of weak and strong disorder, and discuss the non-homogeneous spatial distribution 
of electrons in the critical regime. 

5.1 The wave function 

Following Anderson [1], we expect that an electron is delocalized if disorder is small, and local- 
ized if disorder is strong. To verify this assumption, we calculate numerically all the eigenener- 
gies and eigenvectors of the 2D Ando model, defined by Eqs. J14ll5t . We will see later that the 
Ando model exhibits MIT for critical disorder Wc — 5.838. Thus, for W <C Wc the system is 
in the metallic regime and we expect that the electronic wave function is almost homogeneously 
distributed throughout the sample. On the other hand, the electron should be localized in some 
small region, if » Wc- 

These expectations are confirmed by numerical results presented in Fig. |8]where we compare 
the spatial distributions of two electron eigenstates. The "metallic" eigenstate (W = 2 <C Wc) is 
almost homogeneously distributed within the sample^. Also, the localized eigenstate, calculated 
numerically for the system with disorder = 8 3> Wc, is spatially localized in a certain region 
of the sample, and is orders of magnitude smaller elsewhere. 

The critical point deserves special attention. Figure |9l shows the electron wave function for 
two eigenstates close to the band center for the Ando model with W = Wc- Detailed analysis 
shows that the wave function at the critical point possesses the multi-fractal spatial structure. To 
define multifractality, we introduce the quantities [34] 

r 

In Eq. ( I28> . £"„ is the eigenenergy of the Hamiltonian, and \I>„ (f) is the corresponding eigenfunc- 
tion. It is evident that /i(i?„) = 1 for all eigenstates. The quantity I2 is the inverse participation 
ratio, often used for the characterization of the eigenvectors in disordered systems (Fig. I10> . 

^Note, similar result can be obtained also in the 2D orthogonal systems, if the localized length is larger than the size 
of the system, L <g A. 
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Fig. 8. The spatial distribution of the electron, given by the absolute value of the eigenvector, |'I'(r) with 
the eigenenergy closest to the band center E — for the two-dimensional Ando model. In the metallic 
regime (left, W = 2), the electron is homogeneously distributed in the system, while in the localized 
regime (right, W = 8), the electron is localized in a small part of the lattice and its wave function is very 
small elsewhere. The size of the system is 64 x 64 lattice sites. The critical disorder of the Ando model is 
Wc ~ 5.838. The data was obtained by numerical diagonalization of the random 2D Hamiltonian J6j using 
the LAPACK subroutines. 



We can easily estimate the size dependence of the parameters Iq. If the eigenstate n belongs 
to the metallic part of the energy spectra, then the wave function is homogeneously distributed 
throughout the sample. Assuming that |4'„(r)| ^ L^'^/'^ for all lattice sites r, we immediately 
see that in the metallic regime 



^ ^-(g-l)d_ (29) 



On the other hand, in the localized regime, we expect, in agreement with the right panel of 
Fig. 1^ that the wave function is localized in a certain small region and is almost zero in the rest 
of the system. Then, the wave function |5'n(r)l ~ 1 in the localization region and is almost zero 
otherwise. Inserting this Ansatz into the definition ( I28> we obtain that Iq ^ 1 for all q> \. 

At the critical point, the wave function is multifractal [55]. That means, by definition, that 

~ L-^^-i)-^' (30) 

[34]. In Eq. ( I30I I. < d are fractal dimensions. Note, dq depends on q. 

To understand the meaning of multifractality of the wave function, note that the wave function 
|4'(r^| < 1, for all r. Therefore, the value of |^(r)p^ decreases when q increases. Hence, the 
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Fig. 9. Two electron eigenstates, |\E'(r)p, in the critical regime of the two-dimensional Ando model. The 
size of the system is 64 x 64, and W = Wc ~ 5.838 for the electron energy E — [54]. For a 
given realization of disorder, we calculate two eigenvalues closest to the band center (E = —0.028 and 
E — —0.01476 for the left and right figure, respectively). The spatial distribution of the electron is more 
complicated than in the metallic or localized regime. Detailed analysis of the wave function (discussed in 
Sect. 1141 shows that the wave function ^{f) possesses the multi-fractal spatial structure. 



higher q projects out such sites on the lattice, where is large. For each value of q, these 

sites create a fractal structure. Since different q projects different fractal structures, the complete 
description of the spatial distribution of the wave function requires the knowledge of dq for all 
values of q. 

Figure^! shows values of I2 calculated for all eigenstates of the 3D Anderson model with 
Gaussian disorder W — 2. Data confirm our qualitative estimations, namely that I2 is small in 
the metallic part of the spectra, \E\ ^ \Ec\ ~ 6.58, but possesses values of order of 1 in the 
band tails, where localization is expected. We also see that I2 is a statistical variable, which 
might fluctuate from one eigenstate to another. It is therefore useful to calculate the mean value 
over all eigenvalues lying in a given narrow interval of energies, SE. Also, since I2 might 
fluctuate in many orders of magnitude, we analyze In l2- This is useful especially in the metallic 
regime. Figure^Jshows the probability distribution of In I2 for the three dimensional Anderson 
model. In the metallic regime, the maximum of the distribution p(ln/2) is around the values 
^ L~^, in agreement with Eq. ( I29l l. Also, the distribution decreases exponentially on both 
sides of the maxima, indicating, that the probability to find I2 ^ I rapidly decreases when L 
increases. This agrees with the commonly accepted paradigm, that there are no localized states 
inside the metallic phase. In the localized regime, I2 possess values of ~ 1, as expected. Again, 
the probability to find delocalized state in the tail is very small. 
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Fig. 10. The inverse participation ratio, I2, as a function of energy for tlie 3D Anderson model with 
Gaussian disorder W = 2. The size of the lattice is L = 16. The shaded area indicates propagating 
states, l^l < \Ec\ = 6.58 [56], where I2 is very small. In the tails of the energy band, the electronic states 
are localized and I2 is close to 1. Note, I2 wildly fluctuates and see Fig. llll for the probability distribution 
of /a. 
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Fig. 11. The probability distribution of the logarithm of the inverse participation ratio. In I2, defined by Eq. 
<28t for the 3D Anderson model with Gaussian disorder W — 2 and for system size given in the legend. 
Left: data for all eigenstates with energies \En\ < 0.2 (band center). From the phase diagram, shown in 
Fig.|6| we know that these eigenstates are metallic. We indeed see that the typical values of I2 decrease as 
L~^, in agreement with Eq. ^3^. The distribution decreases exponentially for large values of I2 3> 
indicating that there are no insulating states in the metallic phase. Right: the data for the eigenstates with 
7.5 < \En\ < 7.7. These states are localized (see Fig. |6|. The typical values of I2 are of the order of 1, as 
expected. 
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Fig. 12. Eigenenergies of disordered three dimensional Anderson model with Gaussian disorder. The size 
of the system L = 10. Left: metallic regime. The spectrum is almost equidistant, there are no degenerate 
energies. Right: localized regime. The degeneracy of the energy spectra is more probable than in the 
metallic regime. Middle: spectrum at the critical point. 



5.2 Level statistics 

Other important information about the transport regime can be obtained from the analysis of 
the spectrum of the eigenenergies of the disordered Hamiltonian. In Fig. [21we show a typical 
spectrum of the eigenenergies for the metallic, critical and insulating regime. One sees that the 
three spectra are qualitatively different. This difference can be quantitatively measured by the 
statistical distribution of differences 

Sn = En+1 — (31) 

of two neighboring eigenenergies of the random Hamiltonian (|6j. 

Because of the randomness of the Hamiltonian, s is a statistical variable. It turns out that 
the probability distribution p{s) converges in the limit L oo to three characteristic universal 
functions, depending on whether the system is in the metallic, critical or localized regime. 

If the system is in the metallic regime, then the hopping term between sites causes level 
repulsion, typical for random matrix theory. (Random matrix theory is discussed in Appendix 
O . Random matrix theory states that the distribution p{s) of normalized differences s, defined 
by Eq. i31\ is universal, and depends on the physical symmetry. For orthogonal systems, we 
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Fig. 13. The probability distribution p{s) of the normalized differences, s, between two neighboring 
eigenenergies in the 3D Anderson model. Only eigenstates close to the band center, \En\ < 0.5, were 
considered. Left panel shows the distribution for the model with Gaussian disorder W — 2, which is 
smaller than the critical disorder, Wc ~ 6.1. The system is in the metallic regime and the distribution p(s) 
is close to the Wigner distribution pi (s), given by Eq. <32t . (solid line). The size of the system L = 14 and 
A^stat = 450. Right panel shows data for the box disorder W = 32, L = 10. System is in the localized 
regime, and p{s) is close to the Poison distribution, p(s) = exp —s. 



have 

Piis) f « exp-^s^. (32) 

For unitary and symplectic systems, we have the Wigner distributions given by 
32 4 

P2{s) = ^ exp s^, (33) 

and 

2i8 64 

Pi{s) = -g-^ exp- — s^, (34) 
respectively. Note, all three distributions can be obtained analytically if we assume that 

Pl3{s) = Co s'^ exp-cis^, (35) 

where (3 = 1,2 and 4 for orthogonal, unitary and symplectic symmetry, respectively, and coeffi- 
cients cq and ci are given by the requirement of normahzation. 



dspf}{s) = 1, (36) 
and by the normalization, 

d s spi3{s) = 1. (37) 
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A characteristic property of the Wigner distribution is level repulsion. We indeed see in left panel 
of Fig. [21 that the spectrum in the metallic phase is almost equidistant, and the probability to 
find two degenerate energies is very small. The absence of the degeneracy can be understood 
already in the simple case of the 2x2 matrix, 

A V 



V A 



(38) 



which can represent the two site Hamiltonian. Without the overlap of the wave functions (v — 0), 
the spectrum is degenerate, Ei = E2 = A. However, any non-zero overlap removes the degen- 
eracy, and the eigenenergies become Ei,2 = A ± v. Similarly, in more complicated systems, 
we expect that the overlap of wave functions prevents the eigenenergies being degenerate. In 
random matrix theory, the absence of degeneracy follows from the probability distribution of the 
eigenvalues of random matrices, which is derived in Appendixicl 

The energy spectrum of the insulator is different. The simplest model of localized systems 
is the Hamiltonian given by Eq. (|6jl with V ^ 0. In this limit, the energy spectrum consists of 
random energies e. It is well known [44] that the differences of random uncorrelated numbers 
are distributed with the Poisson distribution. 

Plods) = e-^ (39) 

We therefore expect that in the locahzed regime, the distribution p{s) will be close to the Poisson 
distribution. This is confirmed in the right panel of Fig. \^ Note, however, that the overlap 
F 7^ in any disordered system. Therefore, the level repulsion is always present in the energy 
spectra, independent of the strength of the disorder, and p{s) when s also in the 
strongly localized models [57]. 

In numerical analysis of the level statistics, it is important to note that the mean difference, 
(s), depends also on the density of states. Indeed, (s) would posses smaller values in that part 
of the energy band where the density of states, p{E), is larger. Therefore, we restricted our 
analysis of the statistics p{s) to the energy interval E,E + SE, where the density of states is 
approximately constant. The analysis of the entire energy band requires rescaling the differences 
shyl/piE). 

When the system undergoes the metal-insulator transition, the distribution p{s), transforms 
from the Wigner distribution, given by Eq. i32\ . to the Poisson distribution, given by Eq. J39> . 
There is the third universal distribution, characteristic for the critical point [58]. We do not know 
the analytical form of the critical distribution, but we can calculate it numerically. In Fig. I14lwe 
show the critical distribution Pc{s) for the 3D Anderson model. Data confirms the theoretical 
expectation that Pc{s) must decrease to zero when s ^ 0. From the symmetry considerations, 
it follows that Pc{s) ^ for s ^ 1. There is no agreement about the form of exponential 
decrease for large s ^ 1. Exponential decrease, Pc{s) ^ exp(— s"), with a = 1 + l/{di^), was 
found in Ref. [59], while the semi-Poisson distribution, 

pds) = 4.se-^\ (40) 

was proposed for the orthogonal critical regime by other groups [60]. Numerical analysis [61] 
did not distinguish between these two distributions. 

Figure also confirms the universality of the critical distribution which does not depend 
on the size of the system and on the distribution function of random energies. However, it was 
found that ps depends on the choice of the boundary conditions [60] . 
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Fig. 14. The critical distribution Pc{s) of the normalized differences s. Both the Gaussian and box disorder 
were considered. The size of the system L = 10. The solid circles show data for the Gaussian disorder and 
the system size L = 14. The data confirm that the distribution p{s) does not depend on the microscopic 
details of the model. Also, p(s) is independent on the size of the system. Only the eigenenergies \E„ \ < 0.5 
were considered so that the density of states can be considered as constant. For comparison, we show also 
the Wigner distribution, given by Eq. <32t (solid line) and the Poisson distribution, pioc(s) = exp— s 
(dashed line). 



5.3 Boundary conditions 

The sensitivity of the electron eigenstates to the boundary conditions provides us with another 
criterion for localization [62,63]. Consider a disordered d dimensional system of size L. We can 
calculate all the eigenenergies of the system with periodic boundary conditions, 

+ = V(a;), (41) 
and then repeat the same calculation with the anti-periodic boundary conditions, 

+ = -V'la;). (42) 

The change of the boundary conditions influences the spectrum of eigenenergies, 

En^E„+SE„, (43) 

where SEn is a change of the eigenenergy due to the change of the boundary conditions. Clearly, 
SEn depends on the character of the electron eigenstate. 

Suppose first that the electron in the nth eigenstate is delocalized. Then its wave function 
is spread over the sample. We expect therefore that the change of the boundary conditions in- 
fluences strongly the position of the eigenenergy, £"„, so that SEn might be comparable or even 
larger than the typical spacing between two neighboring eigenenergies. Since the wave function 
is delocalized for any system size, SE will not decrease when L increases. On the other hand, 
If the electron is localized in a certain region inside the sample, then its wave function decreases 
exponentially as a function of the distance from the center. If the size of the system is larger 
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Fig. 15. Dependence of eigenvalues, En{0) on the boundary conditions, given by Eq. <46> . in the metallic 
(left), critical (middle) and the localized regime. Only eigenvalues close to the band center are shown. The 
size of the system is L = 6. 



than the localization length A, the localized eigenstate is almost insensitive to the boundaries and 
th localized electron does not know what happens on the boundary of the system, so that 5En 
decreases exponentially when the size of the system increases, 5E ~ exp — L/A. 

To measure the sensitivity to the change of the boundary conditions quantitatively, consider 
a small interval of energies, E ± 5, and calculate the parameter 



where Ai? is a typical difference between two neighboring eigenenergies. 



and {5E) is a typical change of eigenenergies lying in the interval E ± 5. Then, from the size 
dependence of g, we can distinguish between the metal and insulator. 

In Fig. [21we show how the eigenenergies £"„ of the disordered system depend on the bound- 
ary conditions. We consider the 3D Anderson Hamiltonian, given by Eq. Q with the random 
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on-site energies, £„, and the boundary conditions determined by the real parameter, 0, 

^ix,L,z) = e'''^{x,0,z) (46) 
^ix,y,L) = e'««'{x,z/,0). 

Clearly, 6 — corresponds to the periodic boundaries, Eq. J4U . and 9 — n gives us anti-periodic 
boundary conditions, J42> . 

It was shown in Ref. [62] that if the energy E belongs to the metallic part of the spectra, then 

gr = ctL'^-^ metal (47) 

where a is the conductivity of the sample. When the eigenenergies around the energy E are 
localized, gx decreases exponentially with the system size, 

gj, ^ e"^/^ insulator. (48) 

The quantity gx is called the Thouless conductance. It played an important role in the formula- 
tion of the scaling theory of localization [3]. 

The mean value, {SE), is called the Thouless energy, Et- By comparison of the expressions 
J44I47> with the formula for the metallic conductivity, a ~ e^Dp, we find that 

hD 

^T^jj- (49) 
The corresponding Thouless time, 

TT = h/Er, (50) 

represents the typical time, /D, which the electron needs to diffuse from one side of the 
sample to the opposite side. 

To the best of our knowledge, the relation ( I44t has so far not been used for the numerical 
analysis of the metal-insulator transition. The reason is that both SEn and AE fluctuate not only 
as a function of the microscopic realization of disorder in a given sample, but also as a function of 
the energy within one given sample. Also, it is not clear which averaging procedure - arithmetic 
or geometrical - is more appropriate for the calculation of mean value, (£'„). Nevertheless, an 
introduction of the variable gx provided the first step toward the scaling theory of localization. 
Since it is rather easy to estimate how both differences, 6E and AE, depend on the system size, 
L, we can estimate the size dependence of the Thouless conductance, gr- 



6 Transfer matrix and conductance 



As discussed in the introduction, the electron locahzation affects the transport properties of dis- 
ordered systems. Therefore, the theory of localization concentrates mostly on the electron propa- 
gation. A key role in this study is played by the conductance, a quantity defined by Landauer [4] 
and Economou and Soukoulis [64]. The conductance is expressed in terms of the transmission 
and reflection amplitudes of the electron. 

Consider scattering experiment shown in Fig. ^] The sample is connected on both sides to 
the semi-infinite ideal leads. We are interested in the probability that an electron, coming either 
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from the left or from the right side of the sample, can transfer to the opposite side. Thus, we want 
to characterize the sample in terms of four matrices: transmission of the wave from the left to the 
right, t~^, from the right to the left, t^, and by the reflection coefficient from the right to the right, 
r+, and from the left to the left, . In the most simple case of one-dimensional conductors, all 
the above parameters are complex numbers. 

Transmission and reflection amplitudes define the scattering matrix S: 

Figure[T6lshows that the scattering matrix, S, expresses the wave functions of outgoing waves B 
and C, in terms of the wave functions of the incoming waves, A and D, 

^^^Sfl) ,52, 



Linear relation i52\ can be re-written into the form 

T(t]^ (53, 



D J 

where T is the transfer matrix, which determines the fields on one side of the sample with the 
fields on the another side. An explicit form of the transfer matrix, derived in Appendi?JAlreads 



<+-r-(t-)-V+ r-(t-)-i 



(54) 



Some useful properties of the transfer matrix are given in Appendix IXI If the system possesses 
time reversal symmetry, then the transmission amplitudes t+ and are related by 

t+ = {t-f. (55) 

In the special case of ID system, we have that t+ = Therefore, we omit the superscripts in 
the transmission amphtudes in the following discussion. 




Fig. 16. Definition of amplitudes A - D. 
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Fig. 17. Experimental setup for the measurement of the conductance g. Two semi-infinite leads are attached 
to the sample. Electrons are emitted from the left reservoir, propagate through the left lead and scatter on the 
sample. The transmitted electrons are absorbed in the reservoir on the opposite side, the reflected electrons 
propagate back. There is no scattering in the leads. There are two possibilities of the measurement of the 
voltage. Voltmeter Vl measures the voltage on leads; resulting conductance is then ql = J/Vl, given 
by Eq. <57> . The voltmeter Ves measures the voltage on reservoirs. In this measurement, we obtain 
Economou-Soukoulis conductance Qes, given by Eq. <58> . We show also the voltage for the case when the 
sample is totally transparent to show the voltage drops due to the contact resistances between the reservoirs 
and leads. 



6.1 Conductance 

Consider the one dimensional experimental setup shown in Fig. ^] The sample is connected to 
two semi-infinite ideal leads, which transfer electrons from two reservoirs. The current through 
the sample is proportional to the voltage difference, AV, and to the conductance, g, 

J = g^v. (56) 

To find the conductance, we have to measure the voltage, AT^. Two possibilities were considered. 
Landauer [4] proposed to measure the voltage difference between the leads, and obtained the 
conductance 



_ T 
~ Xl-T' 



(57) 
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On the other hand, Economou and Soukoulis [64] considered the voltage difference on reservoirs 
(Fig.ll7>. and derived an alternative formula, 

gss = ^T. (58) 
h 

Both formulas are equivalent to each other only in the limit of small transmission, T ^ 0. 
However, they lead to different results when the sample becomes transparent (T 1). Landauer 
formula predicts that the conductance diverges (as it should be since the resistance is zero). On 
the other hand, expression ( I58> converges to e^/h. 

The origin of the difference between two formulas lies in the presence of a contact resistance 
between leads and reservoirs. We can write 

' ' l^^^ + l (59) 



h Qes T T h gL 

or, in terms of resistances, 

1 h 

PES = = PL + (60) 

9ES 

Thus, in the measurement of Economou and Soukoulis, the total resistance included in the mea- 
surement contains also contact resistance, pc = h/e^, measured on the contact of leads and 
reservoirs. 

In numerical simulations, we will calculate the conductance gss- When necessary, we can 
extract the effect of contact resistance with the help of Eq. ( I59> . 

Historically, both quantities, and gES^ are called "Landauer conductance" in the literature. 



6.2 Conductance of multi-channel system 

In the previous Section, we have introduced the conductance of the one dimensional systems. 
For a given energy, we have only one possible value of the wave vector. This is the reason 
why the maximum value of the conductance is 1 (in units of /h). In the real world, the leads 
might be quasi one dimensional, which means that the electron can propagate also in directions 
perpendicular to the propagation direction. Since the cross section of the leads is finite, the wave 
vector in the transversal direction, is quantized and possesses only discrete values. Each 
value of the transversal wave vector defines one channel. It might happen that transmission is not 
possible for some values of k^_. We call such channels evanescent, or closed. For simplicity, we 
do not consider evanescent channels in the present discussion. Evanescent waves are discussed 
in Appendix lA.51 The number of channels iVo, defines the size of the transmission matrices, t+ 
and t~ . The physical meaning of the matrix element is clear: t^j, is the transmission amplitude 
of the electron from the channel a to the channel b. Similarly, we define the matrix r+ which 
contains reflection amplitudes from channel a to channel h. Then, 

Tl,^K,? (61) 

is the probability that the electron, coming in channel a, is transmitted through the sample into 
channel 5, and 



^ab ^ \' ab\ 



(62) 
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Fig. 18. Transmission of electron through disordered sample. The sample is connected to two semi-infinite 
metallic leads. Channels are determined by transversal momentum, or, equivalently, by incident angle, 
a = tan~-^(fc„/fe2n)- 



is the probability that the electron is reflected back into channel b. 
It is useful to introduce the probability 

7^a+-E^a.(^:.r=E^.t (63) 
b b 

of the electron coming in the channel a to be transmitted through sample, and 

-T.^a ^T^^ab (64) 

b ah 

is the total probability that the electron transmits through the sample from left to right. Similarly, 

^+=El^aJ' (65) 
ab 

is the probability that the electron is reflected back. Since electrons can not be absorbed in the 
sample, we have 

r+ + i?+ = iVo, (66) 



as is proved in AppendixlAl Then, generalization of expression ( I58t to the multi-channel system 
is straightforward, 

g2 

9ES - -^Tr th, (67) 
[65-68]. 

Multi-channel conductance is fully determined by eigenvalues of the matrix tH. Using the 
parametrization of the transfer matrix, derived in appendix lA.3l we obtain 

g2 ^° 2 ^° 1 

a a 



Numerical Analysis of the Anderson Localization 



29 



where 1/(1 + Aq) is the ath eigenvalue of the matrix and parameters Xa are defined by the 
relation 

Aa = ^ (cosh^Q - 1) . (69) 

It was proved in Refs. [66,67] that in the metallic regime, the conductance g, given by Eq. (I67> 
is related to the conductivity a, given by Eq. ([0, by 

gES = ^L''-^- (70) 

6.3 Relation between the Thouless conductance and gss 

In Sect. I5.3l we introduced the Thouless conductance, g^, which measures the sensitivity of the 
energy spectra to the change of the boundary conditions. In the diffusive regime, gT is given by 
Eq. j47l l. Comparison with Eq. j7m . indicates that two quantities, gr and gss are closely related 
to each other. This equivalence was studied numerically in Ref. [69]. Instead of changing the 
boundary conditions from periodic to anti-periodic, the level curvature. 



c — 



(71) 



was studied numerically and compared with the conductance gss- The curvature, c, measures 
the sensitivity of the energy spectra to the change of the boundary conditions. In the metalhc 
regime, the curvature c is related to gEs by the relation [69] 

gEs = 7rp{E)L^{\c\). (72) 

Also, in the strongly localized regime it was confirmed in Ref. [69] that 

(ln.g£5) « (ln|c|) (73) 

Owing to Eqs. j72l l and j73> we identify the conductance gT with gEs and conclude that the 
conductance gEs not only measures the transmission properties of the sample, but also provides 
us with information about the sensitivity of the wave functions to the change of boundary condi- 
tions. In what follows we discuss only the conductance gss^ defined by Eqs. ( I58> . and use the 
notation 

g = gES- (74) 



7 One dimensional systems 

It is instructive to analyze first the most simple problem, namely the one dimensional disordered 
(ID) chain. The simplest system which exhibits localization is the ID Anderson model. If '^n is 
the wave function of the electron on site n, then the Schrodinger equation reads 



(75) 
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It is well known [5] that all electron states in the disordered ID system are localized. Local- 
ization is characterized by the localization length, which defines the exponential decrease of the 
wave function , 

= exp -L,/\ (76) 

were = an is the length of the system. 

We can calculate the conductance, given by Eq. 065t . It is more useful to start with the 
analysis of statistical properties of the variable x, related to the conductance by 

5 = T (77) 

cosh2(a;/2) 

(in units of e^/h). 

Since aU states are locaUzed in ID, we expect that conductance decreases exponentially when 
the length of the system increases, so that 

x{L,) = 2L,/X, ^ oo (78) 

as is shown in Fig. ^] This follows also from the Oseledec's theorem, discussed in AppendixIDl 
Conductance g is a statistical variable. To describe the transport properties of the systems 
with given disorder, we must study the statistical ensemble of samples, which differ in the mi- 
croscopic realization of random energies, and calculate the mean, (g), and higher cummulants, 
or the entire conductance distribution, p{g). We prefer to study the probability distribution of the 
parameter x, shown in Fig. |20|for a finite system length Lz- We see that the distribution p{x) is 
similar to Gaussian. However, p{x) differs from Gaussian when x is small. Indeed, p{x) 
when a; since no negative values of x are allowed. Detailed numerical analysis confirmed 
(inset ofFig.Hni [20]) that 

p{x) ^ X x < 1. (79) 



x(L) 20 - 




600 
L 



Fig. 19. The size dependence of parameter x defined by Eq. <77> . Energy E = 1 and disorder W = 1. By 
combining the linear fit, x{Lz) — 1.94 + 0.0282L2, shown by solid line, with Eq. <78> we estimate the 
localization length, A = 70.92. This agrees with the analytical estimation, A — [Re 7]^^ for the real part 
of the Lyapunov exponent. Re j{E) = 1/24(4 ~ E'^) = 1/72, given by Eq. <33tl . 
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X 



Fig. 20. The probability distribution, p{x), of the parameter x for the one-dimensional disordered chain of 
the length Lz = 300, disorder W — 1 and E = 1. For comparison, we show also the Gaussian distribution 
with the same mean value, (x) — 10.41 and variance, var x = 14.45 (dashed line). Inset shows how the 
two distributions differ for small x. Note, p{x) ~ x for small x. To calculate the distribution, the statistical 
ensemble of TVstat ~ 10^ samples was considered [20]. 



Note, no negative values of x are allowed. Thus, to first approximation, p{x) can be written in 
the form 



p{x) =(x X e-^^-'^)'/^^ 



where 



and 



a = (x) + 0(1) 



b — a{x). 



(80) 



(81) 



(82) 



Here, a is a constant of order 1. In the limit of — > oo, a ^ 2 but a < 2 for finite L^- For 
instance, a — 1.388 for system analyzed in Fig. |20l 

From known distribution of p{x), we can calculate the distribution of conductance. 



pig) 



d X p{x)S 



1 



9 



or the mean values 



(g") = / dgpig)g^^ 



cosh x/2. 



d x p{x) 



1 

cosh^"a:/2 



(83) 



(84) 
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We evaluate the integral in the r.h.s of Eq. ( I84> by the steepest descent method. The function 

.se 

d r (x-a)^ 



p{x) cosh ^" x/2 possess a sharp maximum around Xn, which solves the equation 



dx 2b 



+ Ina; — 2n lncosh(a;/2) 



= 0. (85) 



X — Xr, 



It is easy to find that a;„ ~ 0(1). Then, neglecting in Eq. ( I84> a;„ with respect to a which is 
oc ^2 3> 1 we find that 

with the constant c„ independent of the length Lj. 

Equation ( I86> agrees for n = 1 with the analytical result derived in Refs. [70,71]. For ti > 1, 
we recover the universality of the moments of the conductance, derived in Refs. [71,72], namely 
that the ratio 

^ = ^ (87) 

id) ci 

does not depend on the system length. 

Also, note that since (g^) ~ (g) 3> (5)^, we obtain that the ratio 

^ g+(.)/4a ^ ^ (88) 

{9) 

increases exponentially when the length of the system increases. Therefore, the mean value, (g), 
is not a good representative of the statistical ensemble. It follows from the derivation of Eq. (1861 
that all the moments of the conductance are determined only by a vanishingly small number of 
samples of the ensemble. These samples with small x are by no means representative. The same 
conclusion can be drawn out from numerical data in Fig. I^which shows the conductance calcu- 
lated for iVstat = 200 samples which differ from each other only by the microscopic realization 
of disorder We see that the conductance fluctuates from sample to sample in many orders of 
magnitude. The mean value, (g), is determined by a few samples with 5^1, while the most 
probable value of the conductance is in many order of magnitude smaller. 

Since x possesses a good probability distribution, and the transmission, T = cosh^^ x/2, it 
is evident that for long systems it is more convenient to analyze the distribution of the logarithm 
of the transmission In T. We introduce a typical conductance, 

gtyp = e<'"3> cxe-'^cxe~<"\ (89) 
Comparing the typical conductance with the mean conductance, 

ln(g) ~ {x)/{2a) (90) 
we obtain that 

Ingtyp - (ln5> = ^ln(g). (91) 



^^The factor '^^^ arises form the normalization constant of the distribution p{x). 
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Fig. 21. Conductance calculated for A'stat ~ 200 realization of disorder on the ID system of the length 
L = 200 and box distribution of the disorder with W = 1. The localization length A « 100. 



The typical conductance is orders of magnitude smaller than the mean conductance. This is typ- 
ical for the localized regime. The conductance, g, is therefore not a good variable for description 
of transport in the localized regime. 

The same holds for the resistance. Using p — R/T, we have 



p = sinh^(a;/2) 



(p") = / dxp{x)smh'^{x/2) « e°"+^"'/^ 



and 



In particular, for n = 1 we obtain, "* 



(PL 



1 



(92) 



(93) 



(94) 



Using the composition law or transfer matrices, derived in the next Section, the mean value of the 
Landauer resistance, (pl), increases exponentially with the length of the system, but the typical 
resistance. 



Ptyp 



increases much slower than the mean resistance: 



Ptyp 



(95) 



(96) 



Again, the reason for the difference between the mean and the typical resistance is that the value 
of (p") is determined by very specific samples which have extremely huge resistance. Detailed 
analysis of the statistics of the resistance can be found in Ref. [16, 73-77] 



7.1 Role of quantum coherence 

To understand the physical origin of the localization and huge fluctuations of the conductance, 
we calculate the transmission of the quantum particle through two barriers, shown in Fig. 1691 

"^We remind the reader that b = aa and assume a = 2 in the limit of — > oo. 
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In Appendix lA.il we derived the transmission amplitude for transmission through two barriers, 
given by Eq. (I228t . 

ti2^tl[^-rtrir't-. (97) 

From Eq. ( I97> we obtain the transmission probability, 

2 T\T^ 

1 + R1R2 - 2/R^cos( 



T12 = |tr2l' = ^-TT^ i /T^,., . (98) 



where we use Ti^2 = \t1.2?, = /RTe^'^i, = v^e'"^^ and ( 
Equation ( I98> can be written in the form 



lnTi2 = InTi + In^, - ln(l + R1R2 - 2^/RjRacos<j>). (99) 

We can introduce the disorder by considering the phase factor, (p, being random with uniform 
distribution in the interval (0, 27r) [16]. Then, averaging over all realization of phase gives the 
relation 

(lnTi2) =lnri+lnT2, (100) 

since the integral of the last term on the r.h.s. of Eq. ( I99l l is zero [16]. 

Relation (fTOOl i can be generalized for the N scattering centers. We obtain 



(InTw) = iVlnTi (101) 



where In Ti is the average of the transmission probability through one scatterer, 

lnTi = -^lnTa. (102) 
Equation dlOU gives immediately the exponential decrease of the typical transmission, 

where the localization length is given by 

X-^ =h[j\/2. (104) 

With the use of the expression pl + I = 1/T (Eq. |94} we recover the relation ( I96> . 

On the other hand. Ohm's law claims that the resistance of the ID system increases linearly 
with the length of the system, or, equivalently, with the number N of scatterers. To obtain this 
result, we have to assume that electron is a classical particle. Then, instead of combination of 
transmission amplitudes, we combine the transmission and reflection probabilities. We obtain 
that 

T T 

T12 = (105) 

I-R1R2 

With the use of relations T = 1 — R and p = R/T,we obtain from Eq. ( I105> that 

Pl2=Pl+P2, (106) 

which is nothing but Ohm's law as we wanted. 
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Fig. 22. Left: Distribution of the logarithm of the conductance, p(lnp) for the ID disordered Anderson 
model. The energy of the electron is i5 = 1, and the disorder W = 0.25, 0.5, 1 and 2. To keep the same 
value of the ratio L/X, we use the analytical result, A ~ W~^, for the localization length, discussed in 
Appendix IdI and choose the length of the system Lz = 4800, 1200, 300 and 75, respectively. Note the 
drop of the distribution at In (7 = 0. This is because the conductance never exceeds 1 in the ID systems. 
Right: p(lng) for disorder W — 2 and two lengths of the system, Lz — 200 and = 400. For a longer 
system, the probability to have g close to 1 is negligible and p{ln g) is Gaussian. 



7.2 Distribution of tlie conductance 

Figurel22lpresents numerical data for the distribution p(ln g) for a strongly disordered ID system. 
The distribution is similar to a Gaussian, in agreement with theoretical expectation discussed in 
Appendix IdI and with our estimation of the distribution p{x). The main difference from the 
Gaussian distribution is that the distribution drops to zero at In g = since the conductance 
never exceeds 1. 

Figure|23]shows the probability distribution, p{g) for ID chains shorter than the localization 
length. For very short systems, <C A, the electron "almost always" propagates through the 
sample. The chance to be scattered is very small. Interestingly, the distribution is almost flat 
when Lz/\ ~ 0.5. 

Numerical data shown in Figs. |S] and confirm that the distribution depends only on the 
ratio Lz/\- Although the system is defined by many parameters - energy, disorder, length of the 
system - the only parameter which really determines the transport properties is the ratio Lz/\. 
This observation simplifies considerably the theoretical investigations of the localization, and it 
is the key assumption in the scaling theory of localization. 

7.3 Ergodic hypotliesis 

As discussed above, the conductance of a ID system is a statistical quantity which wildly fluctu- 
ates as a function of distribution of random energies, Sr- This is demonstrated in Fig. |^ which 
shows the conductance for iVgtat — 200 realization of disorder. These data should be compared 
with those shown in Fig. |23which plots the conductance of a given sample as a function of en- 
ergy E of the electron. We observe that the conductance fluctuates similarly as in Fig.|^ More 
detailed analysis proved that the values g{E) create the same statistical ensembles as the values 
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Fig. 23. Conductance distribution for the ID Anderson model with short system length L^. Figures contain 
data for W = 0.5 (the localization length A = 360) and L = 25, 50, 100, 160, 200 and 300. These data 
are compared with data for W = 0.25 (A = 1440) and a four times longer system, to show that systems 
with the same ratio L/X have the same distribution p{g). The energy of the electron is i5 = 0.5. 




Fig. 24. The conductance of the single disordered one dimensional chain as a function of the energy E of 
the electron. The length of the sample is L = 200, disorder W = 1. Note, the conductance fluctuates 
within the same orders of magnitude as in Fig. 1211 The energy dependence of the conductance is studied in 
details in Ref. [78] 



of the conductance calculated for different samples with fixed energy. This result is known as the 
ergodic hypothesis [13]. 

The ergodic hypothesis plays a crucial role in the experimental studies of the transport in 
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disordered systems. Contrary to numerical simulations, it is impossible to study a huge number 
of microscopically different samples in experiments. Fortunately, due to the ergodic hypothesis, 
it is sufficient to use only one sample and vary the Fermi energy. Similar equivalence holds 
also in the metallic regime. Here, not only Fermi energy, but also the magnetic field can be 
varied [79]. The ergodic hypothesis in higher dimensional systems was numerically tested in 
Ref. [80]. 



8 Diffusive regime 

By definition, the diffusive regime is characterized by the diffusive propagation of the electron 
in the sample. This happens when the mean free path is sufficiently large so that electron is 
scattered separately on two successive impurities. Also, the size of the system must be large 
enough, 

L > ^, (107) 

so that the electron scatters many impurities before it leaves the sample. 

The diffusive regime can be obtained in the 3D systems with disorder W < Wc- Although 
no metallic phase exists in 2D systems, we can observe the diffusive regime if the size of the 
system is smaller than the localization length, 

L < A. (108) 

Similarly, in weakly disordered quasi one dimensional systems, described by the DMPK equation 
(Appendix 151. the diffusive regime exists when the length of the system does not exceeds 
localization length. Note, there is no diffusive regime in one dimensional samples since ^ = A 
[84]. 

The diffusive regime is characterized by the conductivity, a, given by Eq. Q, which is related 
to the conductance, g, by Eq. ( I70> , which we now write in more accurate form, 

(g) = (109) 

In Eq. ( I109> , (. . .) means ensemble averaging. The averaging is crucial in Eq. ( I109> . Contrary 
to the conductivity, a, the conductance g, is a sample dependent statistical variable, which is not 
self-averaged. In Sect. I8.3l we show that fluctuations of the conductance in the diffusive regime 
are universal, independent of the size of the system and the mean value, {g) . 

Owing to Eq. (I107> . the electrons are scattered many times inside the sample, and their 
wave function becomes totally randomized. Thanks to this randomization, electron transport is 
totally universal in diffusive regime. We discuss two main transport properties of the diffusive 
regime: the weak localization correction to the mean conductivity, and the universal conductance 
fluctuations. 



8.1 Weak localization 

In the derivation of the conductivity, a, given by Eq. Q, no quantum interference processes were 
considered. Experimental and numerical results, however, show that these processes play an im- 
portant role already in the diffusive regime, and lead to the so called weak localization corrections 
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Fig. 25. The conductance (in units of /h) in tlie 2D ortliogonal system. Weak localization corrections 
are given by a logarithmic decrease of (g) when the size of the system increases. Symbols show numerical 
data for the conductance, while the lines are logarithmic fits (note logarithmic scale on the horizontal axis). 
Corresponding slopes are given in the legend. Theory predicts that the slope is universal, equal to n^^ — 
-0.318 (Eg.fTTsl. 



to the conductance. We demonstrate the weak locahzation effects in Fig. |53]which shows that 
the mean conductance, (g), of the 2D disordered Anderson model depends logarithmically on 
the system size. This is a consequence of coherent backscattering of quantum electrons. 

To understand the origin of this correction, we calculate the probability that the electron, 
propagating through the sample, is scattered back. Figure|26lshows the typical closed trajectory 
of the electron. The probability amplitude associated with the nth closed trajectory is A„. Since 
the electron is a quantum particle, the probability to return back is given by (in units of e'^/h) 

2 



^9q 



(110) 



where n counts all possible closed paths. The corresponding contribution for the classical parti- 
cle is 



(111) 



Note, 6gc is included already in the conductance, cr, given by Eq. Q. Therefore, the logarithmic 
corrections to the conductance, shown in Fig.|25]originate form the difference 



^.9 = ^9q - 59c- 
The r.h.s. of Eq. ( II lOt can be written in the form 



El 



E 



(112) 



(113) 
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In random systems, the sum of the off-diagonal terms is zero because of random phases of com- 
plex amplitudes A. However, there is an important exception, when A„ and An' correspond to 
the same closed path, but opposite in direction to the electron propagation. Consider one loop, 
as shown in Fig. |26l The electron can travel this path in both directions; both possibilities give 
the same contribution, — = l^-P- Moreover, if the system possesses time reversal 

symmetry (which is the case of 2D Anderson model), then also Aj^ = A^, and the off-diagonal 
terms, 

A+A*_ + A*^A^ = 2\A\'^ (114) 

are real, independent of phase. Then, in the quantum case, a given loop contributes to backscat- 
tering by 

dgg^\A+\^ + \A^\^ + A+A*_+AXA^=A\Af. (115) 

In the classical case, this contribution is only 

6gc^2\A\^, (116) 

The difference between classical and quantum backscattering, given by Eq. dl 12> . gives rise to 
the negative weak localization corrections to the conductance. 
We do not calculate dg here, only give a final result [30], 

Sg - Sg, Sg, ^ -A^'l'^-^-L- f ^. (117) 



h {27r)d J 

Here, q is the wave vector of the electron propagating on the loop. To avoid the divergences 
in integral dl 17l i. we introduce the lower, qmin = ^/L, and the upper, gmax — ^/i, integration 
limits [81]. Then, integration gives for dimension d = 2 the expression 

Sg = --\nL/£, d = 2, (118) 

TT 

and for d = 1 and d = 3 the formulas 

( -7r-^{L/£-l) d = 3 
Sg^l (119) 

[ e/L-1 d = i. 

Since the backscattering is larger for the quantum particle, the correction to the conductance 
must be negative. Note, however, that the above results hold only in systems with time reversal 
symmetry. When a strong magnetic field is applied to the system, the time reversal symmetry is 
broken and Sg = 0. Indeed, the electron acquires on the closed loop the phase 

^n=± j A{r)dr, (120) 

where ^ is a vector potential and the sign ± depends on the propagation direction. Then, the 
product AnA*-^, is not real, but possesses the phase factor e^*"^". Since the phase 0„ depends on 
the length of the loop, in the summation over all loops, all complex contributions cancel. 

In symplectic models, when the hopping between two neighboring sites becomes spin de- 
pendent, the quantum contribution to the backscattering changes the sign. The reason is that the 
wave function of the particle with spin 1/2 transforms into itself after rotation in Att, contrary 
to spin-less particle [25]. The weak anti-localization corrections to the conductance for the 2D 
Ando model is shown in Fig. |^ 
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Fig. 27. The weak anti-localization corrections to the conductance for the 2D Ando model. The conductance 
is given in units /h. The legend presents the calculated slopes. Theory predicts the slope +7r~^ = 0.318. 
More accurate data for the weak anti-localization correction were obtained within the SU(2) model in 
Ref. [82]. 



8.2 Weak localization in quasi-ld systems 



Of particular interest is transport in weakly disordered quasi one dimensional systems. This 
problem is exactly solvable. In Appendix 151 we introduced the DMPK equation for the joint 
probability distribution of parameters A of the transfer matrix. From the DMPK equation, the 
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Fig. 28. The mean conductance, (g) , as a function of N j Lz for the 2D systems (left) and the 3D anisotropic 
model ^9) with t — 0.4 (right). The slope of the linear fit determines, due to Eq. <121> . the mean free path, 
I. Left: The 2D system, iV = L = 50. We found I = 17.0, 4.9, 2.30 and 1.25 for Vf/ = 1, 2, 3 and 4, 
respectively. For W = 1, we add also data for L — 100. Right: The 3D systems of the size L x L x Lz- 
For W — 4: (circles) we found £ = 1.92, while £ ~ 9 for = 2. Note, the weak localization corrections 
in both systems are close to —1/3, as predicted by DMPK equation (Appendix 151. The only exception is 
the 3D system with W = 2, since £ ~ Lin this case. 



following expression for the mean conductance can be derived [83] (in units of /h), 

, . Nl 1 

(.> = ^-3, (121) 

(Eq. I282> . The weak localization correction is constant, independent of the system length, 
provided that ^ <C L < Nl. 

Note that expression (I121> was derived in the limit of infinite number of transmission chan- 
nels, N ^ oo. Therefore, it cannot be applied to the ID models, where = 1. 

Figure|28lverifies both the dependence of the mean conductance and the weak localization 
correction, 5g = l/i for two different quasi-ld systems. Note, relation (I121> can be also used 
for the numerical calculation of the mean free path, I. Numerical data for the conductance in 
quasi-ld systems are presented in Fig. |28] The obtained values of the mean free path are rather 
small, of the order of the lattice period. Only when disorder is very weak {W = 1 in 2D), then 
t K, n. Note that in the Born approximation [84], the mean free path decreases with the disorder 
as 

£(xW-'^. (122) 



The numerical data for the mean free path of 2D systems, obtained in Fig. |28lare in very good 
agreement with results of Ref. [84], shown in Fig. |29] 
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Fig. 29. The numerical data for the mean free path, £, from Fig. l28l (diamonds\ compared with results of 
Ref. [84] (circles). 



8.3 Conductance distribution. Universal conductance fluctuations 

Because of the randomness, the conductance is a statistical variable, and can vary from sample to 
sample as a function of the distribution of random energies e„. In the diffusive transport regime, 
the conductance distribution is Gaussian, as shown in Fig. 1301 

Since the disorder W is weak, the statistical properties of the conductance can be derived 
theoretically, either by the perturbation Green's function method [9] or within the framework of 
the DMPK equation, introduced in AppendixlHI The distribution of the conductance is Gaussian 
with universal width. The variance. 



var.g ^ (g^) - {gf 



(123) 



is a universal number, independent on the disorder strength (provided that inequalities jl07lll08l 
hold), var g depends only on the dimension of the system and on boundary conditions [9]. This 
phenomenon is known as universal conductance fluctuation. 

Also, var g depends on the physical symmetry of the model. It was proved in Ref. [9] that 



1 

var q (X —. 

13 



(124) 



We remind the reader that (3 — 1, 2 and 4 for the orthogonal, unitary and symplectic systems, 
respectively. For the quasi-ld systems, the same universality relation was derived in Ref. [83,85]. 

In systems with higher dimension, var g depends also on the boundary conditions in the 
directions perpendicular to the propagation [9]. This dependence was theoretically investigated 
in Ref. [86] and numerically verified in [87]. 

The deviations of p{g) from Gaussian form can be measured by the third cummulant. 



{g')c^{g')-Hg'){g) + 2{g)' 



(125) 
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Fig. 30. The probability distribution p{g) of tlie conductance for weakly disordered systems. Left: The 2D 
samples of the size 300 x 300. Disorder W — 3 (circles) and W — 2 (squares). Solid lines are Gaussian 
fits with variance 0.168 and 0.179 for Vl^ = 3 and W — 2, respectively. This is close to the theoretical 
prediction, 0.185, obtained for hard wall boundary conditions in Ref. [87]. Right: The 3D system with 
box distribution of random energies, W — 6 and L = 10 (circles) and L — 14 (squares). Solid lines 
are Gaussian fits with mean value (g) = 7 (L — 10) and = 10.03 (L = 4). These values confirm the 
relation (g) — aL with the conductance a — 0.7. The width of distributions is var g = 0.38 and 0.366 for 
L = 10 and L = 14, respectively, are close to the theoretical prediction 0.314 [87,88]. Data were obtained 
numerically with statistical ensembles of A'atat = 10 samples. 

It was shown analytically [89] that {g^)c — for the quasi-one dimensional systems and is 
~ 10~^{g)^^ for the 3D system. Such small values are not observable with today's computer 
facilities. 

The conductance fluctuations for the 2D weakly disordered systems are plotted in Fig. |2] 
Note how the variance of the conductance is universal only when the size of the sample is much 
larger than the mean free path, L > 10 x £. However, L must be much smaller than the local- 
ization length, L/X ^ 1. Also, for small system size, when condition ilOli is not satisfied, the 
system is in the ballistic regime, in which the scattering of electrons is not sufficiently strong 
to randomize the wave function sufficiently. In this regime, the conductance is large, g ^ No, 
and the fluctuations of the conductance increase. Of course, the variance decreases to zero in the 
limit of a regular lattice (W = 0). 
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HW 


periodic 


1 


2/15 


2/15 


2 


0.185613 


0.154078 


3 


0.314054 


0.2194393 



Tab. 1. Universal values of the conductance fluctuations, var g = (g^) — (g)'^ for the d-dimensional 
disordered systems with the hard wall and periodic boundary conditions [88]. For other symmetry classes, 
var g must be divided by a factor of /3, (Ea. ll24t . var g ^ in the dimension d — 4 ~ e. 
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Fig. 31. Conductance fluctuations for the 2D Anderson model. Different symbols correspond to the disorder 
W = 1, 2, 3, 4 and 5. The size of the system increases from L = 10 to L = 1000. The horizontal dashed 
hne is the theoretical prediction var g ~ 0.1856 of the universal conductance fluctuations for 2D orthogonal 
system with the hard wall boundary conditions [87]. For W = 1, the mean free path ^ ~ 17 is already 
comparable with L for some smaller samples. The system is in the ballistic regime and var g exceeds 
universal value. A more detailed analysis of the universal conductance fluctuations has been done in [87]. 



8.4 Other universal relations 

Universal conductance fluctuations are the consequence of the universality of the diffusive trans- 
port. In Appendices IbI and Icl we discuss the statistical properties of the eigenvalues Aa of the 
matrix [t^t\ . It turns out, that in the diffusive regime the joint probability distribution of the 
eigenvalues A is universal and depends only on the length of the system and the number of the 
transmission channels. No- The disorder, W, influences the transmission parameters only in 
terms of the ratio L^/A of the system length and the mean free path. 

From the universality of the probability distribution p{X) we can obtain other universal rela- 
tions for the transmission parameters. In Appendix|Clwe show that the spectrum of parameters 
Xa, defined by Eq. ( I68> is linear. 



Xa oc a, 



(126) 



as shown in Fig. ^61 Also, the differences, Xa+i — Xa are distributed with the Wigner probability 
distribution (Fig.l74>. 

Here we want to show that not only the conductance, g, but also the transmission probabili- 
ties. Tab and Ta, introduced in Sect. l6.2l as 



Tab — \tab 



(127) 



and 



Ta — ^ \tabf — Tab 



(128) 
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Fig. 32. Probability distribution p(sa), given by Eqs. jl31ll32> for various values of the mean conductance, 
(9). 



exhibit universal statistical properties. We remind the reader that Tab is the probability that the 
electron, coming in channel a is transmitted through the sample and leaves the sample in channel 
b and Ta is the probability that the electron, coming in channel a, transfers through the sample. 
It is clear that 



9 



(129) 



Note, Ta is not the eigenvalue of the matrix tH. 

The universal probability distribution for the normalized transmission probability, 

Ta 



(Ta) 



(130) 



was derived in Refs. [94,95]. The distribution p(sa) is determined only by the mean conductance, 
and is given by the following analytical formula. 



P{Sa) 



dx 



-e 



where 



^x) = {g) [in (v/l + x/(<?) + 



(131) 



(132) 



The probability distribution p{sa) is shown in Fig. |32]for a few values of the mean conductance. 
Distribution p(sa) possesses universal variance. 



var Sa = (si) - (s„,)2 = 



3(5)' 



(133) 
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Fig. 33. Upper panels show the probabihty distribution of the parameters Sa, defined by Eq. jl30t for the 
2D Anderson model, with the random binary potential given by Eq. <134> . The energy of the electron, 
E = 0.14 (left upper panel) and E = 0.31 (right upper panel) is measured in this particular case from 
the bottom of the conductance band. Different symbols represent the data for channels a — 1, 2, 3 and 
10. The solid line is the theoretical prediction, given by Eqs. <131ll32t . Two lower panels show the mean 
transmission, (Ta) and the variance, var Sa, given by Eq. <133t . Note that all channels are almost equivalent 
and they give the same contribution to the conductance. The number of open channels is A'^o = 23 and 
A''o = 34 for £ = 0.14 and E = 0.31, respectively. The size of the system is 192 x 192 [93]. 



Universality of the distribution il3H was confirmed experimentally in experiment with mi- 
crowave electromagnetic waves [24,95]. Numerically, it was studied in Ref. [93] for the 2D 
system with correlated binary disorder 

p{e) = {l~x)S{e)+x5{e + Vb). (134) 

Spatial correlations were introduced so that random energies create randomly distributed poten- 
tial wells of the size of 3 x 3 lattice sites. 

Figure|33]shows the numerically obtained data for the system of size 192 x 192 and for two 
energies of the electron. Data confirm that indeed the channels are equivalent to each other, and 
they give the same value of the transmission, (Tq). This is a typical property of the diffusive 
regime: transmission in a given channel does not depend on the incident angle, since electron, 
being many times scattered inside the sample, forgets the initial direction of propagation. 

8.5 Diffusion 

Conductivity, a, can be expressed through the diffusion coefficient, D, a = e^Dp, given by 
Eq. n The density of states, p{E) was calculated in section ITTI The diffusion coefficient can 
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Fig. 34. The time dependence of tlie wave packet in the 2D disordered system, defined by Hamiltonian 
Shown are points where |*(r)| > 0.0005 (gray), [^(r)! > 0.0010 (dark gray), and |*(r)| > 0.0050 
(black) at time t — 100, 500 and 900 (from top to bottom). The time is measured in units h/V with V — 1 
and the size of the lattice is 512 x 512. For weak disorder, the wave function of the electron diffuses and 
)r'^{~ 2Dt where D is the diffusion coefficient. Clearly, D decreases when disorder increases. For W — 6 
(right column), diffusion already stops and the electron is localized (see also Fig.|2j. 



be calculated numerically solving the time dependent Schrodinger equation, Eq. (|5|l. Figure l34l 
shows the time evolution of the single electron wave function in disordered 2D samples with 
different strength of the disorder. We use the same data to calculate the width of the wave packet, 

{r^{t))= J dr-^*{f,t)r'^^{f,t) (135) 

in time t. When disorder is small, the electron diffuses from the center, and we obtain that 

(r^) = 2Dt. (136) 

In the case of strong disorder, (r^) should converge to the time-independent quantity, when the 
electron is localized. This was shown already in Fig.^for the case of strong disorder 

Figure |35l shows the time dependence of (7-^) for various strengths of the disorder. The 
diffusion constant, D, is obtained as the slope of the linear dependence (r^) vs t. If the density 
of states, p{E) is known, then we can obtain the conductance, 

5 = a = "^Dp. (137) 
h 
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Fig. 35. Left: (r^) vs time t for the 2D Anderson model witli box disorder The slope of the linear 
dependence determines the diffusive coefficient D. The data for disorder W — 2 and W = 4 corresponds 
to the wave function shown in Fig. |34| For weak disorder, W = 1, we see already the saturation of (r^) 
due to the finite size of the system and reflection of the wave packet from boundaries. Solid lines are linear 
fits, (r^) — 2Dt. Time is measured in units h/V. The diffusion coefficient, D, measured in units a^V/Ti, 
is given in legend. 

Using numerical data for the disorder W = A: p{E) « 0.15/a^y (Fig.|4|l, and D 1.69 a?V/h 
(Fig. I35> . we obtain that g[W = 4) = 1.59(e^/ h). This is consistent with numerical data shown 
inFig.Ea 

8.6 Beyond the diffusive regime 

As will be discussed in next Section, there is no true metallic regime in 2D orthogonal systems. 
What we have observed so far is the diffusive propagation of the electron on a finite square 
lattice of size L x L. When L increases over the localization length, A = X{W), the electron 
becomes localized and the L dependence of the conductance changes from the logarithmic to the 
exponential. 

In Fig. |36lwe demonstrate how the conductance distribution, p{g), of the 2D disordered 
systems changes when the size of the system, L, increases. For disorder W = 3, p{g) is Gaussian 
even for L = 1000. However, for the disorder W = 4, the localization length A w 50 and the 
probability distribution changes its form considerably. 

In the opposite limit of small L, the system is in the ballistic regime. In the limit of i < ^, 
there is almost no scattering inside the sample, so the conductance (g) is given by the number of 
open channels. In 2D, this means that (g) cx L. Increase of the conductance with the system size 
is shown in in Fig. |^ Note, the ballistic regime can be observed also in the ID systems. Fig. 
|23]shows that for sufficiently short systems, the conductance of ID disordered system is close to 
1 when L ^ ^. 

Since the electron is scattered only on a few impurities during its travel through a ballistic 
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Fig. 36. The system size dependence of the conductance distribution, p{g) for the 2D systems. Left: 
W — 3. The conductance is almost independent to the system size, since system is in the diffusive regime 
(£ ~ 2.3, and A — 5046). Right: W = 4. The distribution p{g) is Gaussian only for very small systems 
(L — 10), and it changes considerably when L increases. The mean free path, ^ ~ 1.3 and localization 
length, A = 481. The estimation of the localization length is taken from Ref. [6]. 
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Fig. 37. The length dependence of the mean conductance, (g), for the 2D disordered systems. When 
disorder is weak, (g) increases for small L since the system is in the ballistic regime. Naively, this might be 
interpreted as a metal-insulator transition: since the mean conductance increases for W = 1 and decreases 
for W — 5, one might conclude that the 2D system possesses the critical point, Wc < 1, i. e. that disorder 
W = 1 corresponds to the metallic regime. 
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sample, the values of the conductance for a given sample depends on the actual distribution of 
impurities. Consequently, fluctuations of the conductance increase in the ballistic regime [96]. 
This is confirmed in Fig. |^ 



9 Scaling theory of localization 

As discussed in the Introduction, the metal-insulator transition resembles the phase transition in 
statistical physics. We would like to develop the theory, similar to the renormalization group 
theory of second order phase transitions. 

The first step in this direction is to determine the relevant order parameter It turns out that 
the most suitable candidate for such parameter is the Thouless conductance, gr- First, it can be 
defined both in the metallic and localized regime. Second, the analysis of the sensitivity to the 
boundary conditions tells us how the conductance develops when the size L of the system 
increases. Since gx is equivalent to qes [69] both in the metallic and in the insulating regime, 
we will consider the conductance gEs ™d use the simpler notation, g. 

The main assumption we have to accept is that for the system size sufficiently large, the 
length dependence of the conductance is given only by the conductance itself: 

The "sufficiently large" system size L means that L exceeds all "natural" lengths which might 
determine the transport in the system. Some examples of such scales are the mean free path, 
£, which determines the coherent scattering on impurities, the coherence length of the random 
potential, or magnetic length, which is important if the magnetic field is present. If L is much 
larger than all these scales, we expect that microscopic details of the model become irrelevant. 

Now, we want to derive some consequences from the Eq. ( I138t . First of all, we must admit 
that the form of the function f3{g) is unknown. However, we can derive its limits for g — > oo and 
g ^ 0. Consider for simplicity only the orthogonal symmetry, so that the symmetry parameter 
(3 — 1. In the limit of large conductance, we can use the L-dependence of the conductance, 
derived in Sect. [S] The leading term of the conductance behaves as g = gU^^^. Inserting into 
Eq. (I138> . we obtain that 

lim I3{g) = d-2. (139) 

g-*oo 



Since the first correction to the conductance, derived in Sect. 18. H is negative, we immediately 
see that /3(ln g) reaches the limit ( I139> from below. In particular, for d = 2, we obtain from Eq. 
(ins that 

9 In Q 1 do 1 

/3 5 -TT^ = -^= ■ (140) 

a In L g omL ng 

In the opposite limit, g <C 1, we have exponential localization, 

g ^ e-2VA^ (141) 



so that 

2L 

f3{g) = -—=lng. (142) 
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Now, we can interpolate between the limits ( I139t and (I142t . This is straightforward if we 
assume that the function P{g) is always continuous and monotonous. Although we do not have a 
rigorous proof that P{g) really fulfills these conditions, there is no physical reason not to accept 
them. Then, connecting both limits, we obtain that the function i3{g) behaves as shown in Fig. 



The form of the function f3{g) has some important consequences. First, we see that for d < 2, 
the function /3(g) is always negative. So starting with a sample of finite size, Lq, and conduc- 
tance, go, the conductance, given by Eq. ( I138> decreases and develops to a smaller conductance 
when the system size, L, increases. Consequently, an infinite system with dimension d < 2 does 
not exhibit the Anderson transition since all electronic states are localized, independent of the 
strength of the disorder 

For d > 2 {d ~ 3, for instance), there is a critical point gc such that 



When starting exactly with g = gc, the conductance remains constant, independent of the system 
size even in the limit of L oo. Thus, a disordered system in dimension d > 2 possesses 
a critical point. This critical point is unstable: when starting with g = gc + Sg, we end up, 
in the limit of L ^ oo, in the metallic regime with finite conductivity a. Similarly, starting 
with g = gc — Sg, the system will develop into the insulating regime with exponentially small 
conductance. 

What remains is the calculation of the function f3{g). Analytically, this is possible only for 
systems close to the critical dimension: = 2 + e, (e <C 1) [97,98]. We will discuss these results 
in SectionHjl 

The scaling theory of localization, formulated first in Ref. [3] represents the main milestone 
in our understanding of the Anderson transition. First, it estimates the lower critical dimension. 



f3{gc) = 0. 



(143) 



= 2. 



(144) 



/3(5) 
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\ng 



Fig. 38. The function f3{g) for the disordered orthogonal systems with dimension d = 1, 2 and 3. Notice 
the critical point, gc for d = 3. 
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There is no metallic regime in the disordered system with dimension d < 2. However, this 
statement is true only for orthogonal systems, i.e. for systems with time reversal symmetry. We 
know akeady that in symplectic systems the first correction to the conductance in the diffusive 
regime is positive (Fig. I27> so that (3{\ng) 0+ when g ^ 1 [10]. Consequently, the 2D 
systems with spin dependent hopping exhibit the metal insulator transition. This was for the first 
time predicted in Refs. [47, 48] and confirmed numerically by many authors [49, 54, 99, 100]. 

Soon after formulation of the scaling theory [3] it became clear that 5 is a statistical quantity, 
which is not self-averaged in the limit of the infinite system size. In the metallic regime, the 
probability distribution of the conductance, p{g), possesses an universal width, given by the 
universal conductance fluctuations. In the insulating regime, g does not represent the statistical 
ensemble, as discussed in Section^ and we have to use its logarithm, In g. Even then, ^(In^) is 
Gaussian with the mean value (In g) = —2L/X and the variance, (In^ g) — (In 5)^, being of the 
same order as — (In g). We do not know the analytical form of the critical distribution, Pc{g), but 
both theoretical [75, 101] and numerical [102] results confirm that Pc{g) is independent on the 
size of the system. This is consistent with the size independence of the critical conductance, given 
by Eq. ( I143t . However, since the width of the critical distribution is non-zero, the conductance 
is not the self-averaged quantity at the critical point. The statistical character of the conductance 
opens new problems in the scaling theory. First, we have to prove that both mean values, (g), and 
(In g) obey scaling relations ( I138I I. Then, it would be useful to prove the same for all cummulants 
of the conductance. This is, of course, an unsolvable task. 

The first step in the verification of the single parameter scaUng theory is to understand the 
statistical properties of the conductance in the critical and localized regime. 



The statistical properties of the conductance at the critical point were discussed by Shapiro, 
[75, 101]. With the use of the Migdal-Kadanoff renormalization, he studied the size dependence 
of the conductance distribution and he proved that the critical conductance distribution, Pc{g), 
is universal, independent of the system length. The later works [103] however showed that the 
Migdal-Kadanoff renormalization overestimates the conductance fluctuations, and is therefore 
not suitable for a quantitative description of the critical conductance distribution. 

In dimension d = 2 + e, the non-universality of higher order conductance cummulants has 
been found analytically [12] to be 



Expression J145t states that the higher order cummulants, n > no = e^^, depend on the size 
of the system. This seems to be in contradiction with universality of the critical distribution. 
This discrepancy was explained in Ref. [104]. Starting from the known cummulants, given by 
Eq. ( I145> . the critical conductance distribution was derived. It was shown that Pc{g) is indeed 
universal in the limit of the infinite system size. The non-universality of higher cummulants 
follows from the form of the critical distribution. For small e, the bulk of Pc{g) is approximately 
Gaussian near the mean value (g) . The parameters of the Gaussian peak. 



10 Statistical properties of the conductance in the critical regime 




n < no = e 
n > e^^. 



1 



(145) 



(5) - 



(146) 
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Fig. 39. The critical conductance distribution for tiie 3D Anderson model. The data for the samples L'^ 
with 10 < L < 22 are shown. Statistical ensembles of A^atat = lO'' (20000) sample for L = 10 
(L = 22), respectively were used to create the distribution. Left: the distribution of the conductance, right: 
the distribution of the logarithm of the conductance. The legends present the mean values, which should be 
independent of the system size at the critical point. Note the non-analytical behavior of the distribution at 
5=1- 



and 

vai-5-e°, (147) 

but Pc{g) possesses long tails for large values of g, 

Pc{g) ^ (148) 

The power-law behavior of Pc{g) explains the non-universality of higher cummulants, (^g"), 
which are not defined for n > 2/e. 

These analytical results were derived only for dimensions close to 2, for e ^ 1. Any attempt 
to apply them to the 3D system fails. For instance, Eq. ( I148> predicts that pdg)^'''^^^ ~ g 
for large g, which is clearly not realistic. The most relevant information about the form of the 
critical conductance distribution was therefore obtained from numerical simulations [102, 105- 
108] based on the formula 

g^Ti-th^y^ , (149) 

^cosh2(xj2)' 

already used in the diffusive regime. 

In Fig. |39]we plot the critical conductance distribution for the 3D Anderson model. We see 
that Pc{g) does not depend on the system size. This is what we expect if the scaling theory really 
works. 

The form of pc{g) is rather unusual. Numerical data enables us to understand its main proper- 
ties. It turns out that it is useful to study also the distribution of the logarithm of the conductance, 
p(ln g), shown in the right panel of Fig. |39] 



54 



P. Markos 




Fig. 40. Details of tiie critical probability distribution for the 3D Anderson model. The left figure proves 
that Pc{g) — > when g 0. The right panel shows both the distribution, Pc{g) and its first derivative, 
dpc{g)/dg in the neighborhood of g = 1. Numerical data indicate the non-analytical behavior of the 
distribution at g = 1. 




a 



Fig. 41. The index dependence of the parameters Xa at the critical point in the 3D Anderson model. The 
data confirm that (xa)'^ oc a. Left inset shows that the variance, var Xa decreases as a function of index a. 
The solid line is a fit ao + aix°' with a = 0.51. This gives an estimate of the variance, var Xa ~ a~^^^. 
Right inset shows the probability distributions p{xi), p{x2), p{x3) and p{xe)- Note, the distribution p{xi) 
is similar to the Wigner distribution, pi(x), shown by the shaded area. 



First, we need to know the form of Pc{g) for small g. The behavior of Pc{g) for small g can 
be easier obtained from the distribution of the In g. In the right panel of Fig. |39]we see that 

lnp,(ln5)~-(lng)2 g « 1. (150) 
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Since pc{g)dg = Pc{ln g)d{\n g), we immediately obtain that [105] 

oc exp [-(ln.g)2 - Ing] 5 < 1, 



(151) 



so that Pc{g) when 5 — > 0. The left panel of Fig. I^shows details of the critical distribution 
for g < 0.1. The data, collected from the ensemble of the A'gtat — 10^ samples of the size 
i = 10 confirm that Pc{g) indeed decreases to zero when g 0. However, the probabiUty to 
obtain small values of g is very small. 

Other properties of the critical distribution, namely the form of the tail for large values g 
and the non-analytical behavior in the neighborhood of g — I, are more convenient to analyze 
in terms of parameters Xa [102, 105], introduced in Eq. ( I68t . We remind that Xa determine 
the eigenvalues of the matrix tH and that the conductance is expressed in terms of Xa by the 
following formula 



Numerical data [91, 102, 128] showed that the spectrum of Xa consists of two qualitatively dif- 
ferent parts. In the lower part of the spectra, for a < L, (xa) are independent on the size of the 
system 



This is what we expect, since the mean conductance, (g), is size-independent at the critical point. 
The size independence of Xa was confirmed numerically in Ref. [128] and will be discussed in 
Sect. 112.31 Contrary to the metallic regime, where (xa) oc a (Eq. I320t . at the critical point we 
obtain (Fig. ITlTl that 




(152) 



{xa) = const 



a < L. 
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Fig. 42. Comparison of the critical conductance distribution with distribution p{gi) of the contribution of 
the first channel. 
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It is important to note that Eq. il54\ is valid only for a < L. The upper part of the spectra, with 
a > L, is L dependent with (xa) oc a. Since this upper part of the spectra does not contribute to 
the conductance, we will concentrate only on the lower part, given by Eq. ( fT54t . 

Insets of Fig.lTTIshow that the distribution of xi is similar to the Wigner distributional, and 
all higher parameters Xa, with a > 2 are distributed according to the Gaussian distribution with 
decreasing variance, 

var Xa ^ a^", (155) 

where the exponent a is close to 1/2 [91]. 

We use the statistical properties of Xa to estimate the behavior of the critical distribution for 
large g. From the expression for the conductance, given by Eq. ( I152l i we see that we can obtain 
a large value of the conductance, 

gr^N:^l, (156) 

only when all a;Q, a = 1, 2, . . . , iV are small. However, since all parameters Xa are distributed 
with the Gaussian distribution, we obtain that the probability to have x^ <C 1 is 

exp-i^^~exp-gi+". (157) 
var xn 

In Eq. ( I157> we used that {x^y (x N — g and we estimated the variance var x^ N~°' with 
the use of Eq. ( I155t . Eq. ( I157> confirms that Pc{g) decreases faster than exponentially for large 
g. This is confirmed by our numerical data. The probability to have g > 1 is very small (less 
than 3%). The higher values of g appear with marginal probability. For instance, in an ensemble 
of iVstat = 10^ samples, we found only 470 samples with conductance g > 2 and only one with 
g>3. 

The most surprising property of the critical conductance distribution is that Pc{g) seems 
to be non-analytical at g = 1. The right Fig. 1401 shows that the first derivative, dpc{g)/dg, is 
discontinuous atg— 1. The origin of the non-analyticity can be again explained from the spectra 
of the parameters Xa- The numerical data, for mean values, (xi) — 3.42, (X2) = 5.52, and 
{xs) = 7.07, indicate that the conductance is given mostly by the first term, gi ~ cosh~^(xi/2). 
Indeed, the distribution p{gi) is a very good approximation to Pc{g) for g < 1, as shown in 
Fig. El Since gi < 1, the function p{gi) has a cutoff when g ^ 1^. It is difficult to estimate 
the contribution of higher channels but it is reasonable to expect that the non-analyticity of the 
distribution survives also when the contribution of higher channels is included. This is supported 
by analytical calculations of Muttalib et al [ 1 09, 1 1 0] , who reported the non-analytical form of the 
conductance distribution in the weakly disordered quasi- Id systems with the mean conductance 

(g) - 1. 

The above properties of the critical conductance distribution can be found also in other mod- 
els, where the critical conductance is close to or less than to 1 . For instance. Fig. l43lshows Pc{g) 
for the 2D Ando model. The non-analyticity at g = 1 is clearly visible, too. 



10.1 Properties of the critical conductance distribution 

The properties of the critical conductance distribution are important for the formulation of the 
scaling theory. Besides the system size independence of Pc{g), we need to understand its univer- 
sality, i.e. how the form of Pc{g) depends on microscopic details of the model, on the physical 
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Fig. 43. The critical conductance distribution for thie 2D Ando model. The different symbols correspond to 
the size of the system L = 82 and L = 200. The data prove that Pc{g) is system size invariant. Note the 
logarithmic scale on the vertical axis. Inset describes the non-analyticity of the distribution at p = 1. The 
mean value of the conductance is (g)c = O.Yle^/h. [54] 



symmetry and on the dimension of the system. Especially the dimension dependence of Pc{g) is 
important for the comparison of the theoretical predictions with numerical data. 
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Fig. 44. The critical conductance distribution for the 2D Ando model with the hard wall and periodic 
boundary conditions in the direction perpendicular to the boundaries. The critical conductance is (g) = 
0.655 (0.71) e^/h, and the variance vea g = 0.43 (0.36) (e^/h)^ for the hard wall (periodic) boundary 
conditions, respectively. 
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Fig. 45. Definition of three fractals used in the numerical analysis of the conductance distribution of fractals. 
Shown are the 2nd and the 3rd generations. All three fractals have the same fractal dimension, df — 
In 3/ In 2 ~ 1.58. Fractals A and B have spectral dimension d'g — 1.365, and fractal C has d'g — 1.226. 
Note the different number of nearest-neighbor lattice sites. The bifractal lattice is created by combination 
of the fractal with the linear chain in the perpendicular direction. The spectral dimension of the bifractal is 
da = d's + 1. The number of lattice sites on fractal is 3" in the nth generation, and the length along the 
linear chains is 2" lattice sites. 



As expected, the critical conductance distribution is independent of the size of the system. 
This is confirmed by numerical data shown in Figs. |39]and|^ For a given universality class, 
it does not depend on microscopic details of the model. This was confirmed in Refs. [111,112] 
where the critical distribution, Pc{g) was calculated for the 2D and 3D models with various 
distributions of random disorder. However, the form of Pc{g) does depend on the symmetry 
class [107]. Also, the form of Pc{g) depends on the topology of the lattice. For instance, Pc{g) 
for the 3D Anderson model with triangular, hexagonal and rectangular lattice in the transversal 
direction possess different critical distributions Pc{g) [113]. 

Surprising was the observation that pdg) depends also on the boundary conditions in the 
directions perpendicular to the propagation [86, 114, 115]. This result is, however, consistent 
with our understanding of the conductance as a measure of the sensitivity of eigenenergies to the 
boundary conditions [60]. We remind the reader that also the values of the universal conductance 
fluctuations depend on the boundary conditions. Nevertheless, the difference in the form of Pc{g) 
is surprisingly large for the 3D Anderson model [115], Also, the mean values (g)c, are consid- 
erably different, being 0.445 (0.280) (e^/h) for the periodic and hard wall boundary conditions, 
respectively [115]. In Fig. |^we demonstrate the sensitivity to the boundary conditions for the 
2D Ando model. 



Numerical Analysis of the Anderson Localization 



59 




Fig. 46. Left: The critical conductance distribution for three bifractals, shown in Fig. 1451 The legends 
give the spectral dimension, d^, of the lattice. Critical distribution is close to Gaussian distribution when 
d ^ 2. The power-law tail, predicted by the theory [104], is not observable in numerical simulations. Note, 
although two bifractals, A and B, have the same spectral dimension, the corresponding critical distributions 
differ from each other. This is because Pc(ff) depends on the topology of the lattice. Note also, the non- 
analyticity of the conductance distribution around g — 1 disappears, since more than one propagating 
channel contribute to the conductance in these models. Right: Distribution of pc(ln g) for the 4D Anderson 
model. 



10.2 Dimension dependence of the critical conductance distribution 

As discussed above, the analytical theories provide us only with the information about the critical 
conductance in systems with dimension close to the lower critical dimension, d = 2 + e. Since 
we do not expect that the analytical results could be applied also for e = 1, the numerical 
data for dimension d = 3 are not suitable for verification of the theory. Therefore, it might be 
interesting to calculate the critical conductance distribution on lattices with fractal dimension 
close to dc = 2. This was done in Ref. [113]. 

By definition, the bifractal lattice [116] is linear in the direction of propagation and possesses 
the fractal transversal structure. Three fractals, discussed in Ref. [113] are shown in Fig. |^ For 
the analysis of critical phenomena on fractals, it is important to note that the "dimension" which 
is important for description of the critical phenomena is the spectral dimension of the lattice, 
not the fractal dimension, df. We remind the reader that spectral dimension, ds, determines the 
low-frequency behavior of the phonon density, 

Pphonon('^) ~ UJ''" (158) 

[116, 117], while the fractal dimension, df, determines an increase b'^f of the "mass' (number 
of lattice points) when the scale of the system changes by a factor of b. Note that the spectral 
dimension of a regular system equals to its dimension, ds — d. 

All three fractals, shown in Fig. |^have the same fractal dimension, df = In 3/ In 2 = 1.58. 
Fractals A and B have the same spectral dimension, d'^ = 1.365, and fractal C has spectral 
dimension d'^ = 1.226, so that the spectral dimension is dg — d'^ — 2.365 (2.226) for bifractals 
A and B (C), respectively. In Fig. |60|we will show that the critical exponent, ly is indeed a 
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function of the spectral dimension only. 

Critical conductance distributions are shown in Fig. |46l We see that the mean conductance 
increases when dg 2+ and the distribution pdg) is similar to Gaussian, in agreement with 
theoretical prediction [104]. However, the width of the distribution, measured by the variance, 
var g, increases when e ^ 0. This seems to be in contradiction with the analytical formula (I147> . 
The numerical data cannot confirm the power law decrease of the distribution, given by Eq. ( I148> 
for large values of g. 

For completeness, we show in the right panel of Fig. |46lthe critical distribution pc{\ng) for 
the 4D Anderson model. The distribution is similar to that for the 3D system, including the non- 
linearity at g = 1. Similar to 3D models, the main contribution to the conductance is given by 
the first channel. 



11 Localized regime 

It is commonly beUeved that the distribution of the logarithm of the conductance is Gaussian, 
with the mean value 

(ln.g) = -2L/X. (159) 
The variance, 

var In 5 cx — (Ing) (160) 

should be proportional to — (In g) , with the coefficient of proportionality close to 2. This behavior 
is deduced from the analysis of the DMPK equation, discussed in Appendix 151 It is argued 
that since all parameters Xa oc L^, the difference between X2 and xi becomes large enough 
so that the contributions g2, gs ... of higher channels become negligible in comparison with 
gi = cosh~^(a;i/2) = exp— xi. Since xi possesses the Gaussian distribution in the limit of 
large L^, it is expected that the logarithm of the conductance, Ing w Ingi should posses a 
Gaussian distribution, too. 

This argument is supported by the numerical data for the normalized difference, S21 — 
X2 — xi, shown in Fig. ^\ As expected, p{62i) differs considerably from the Wigner distri- 
bution, and is more similar (although not identical) to the Poisson distribution [20]. Although the 
deviations from the Poisson distribution indicate that X2 and xi are not statistically independent, 
we expect that their correlation is small. Then, the contributions gi and 92 can be assumed to be 
almost statistically independent and we can split the spectrum of the transfer matrix to the first 
eigenvalue, xi, which determines the magnitude of the conductance, and the rest of the spectra, 
which have almost negligible contribution to g. 

However, the above arguments are valid only in the quasi- Id weakly disordered systems (see 
Appendix lB.3t . As shown in Fig. |48l the probability distribution p{\ng) for the strongly disor- 
dered samples is not Gaussian. The origin of the difference between the strongly disordered 3D 
systems and weakly disordered quasi- Id systems lies in the spectra of parameters Xa- Contrary 
to the weakly disordered systems, parameters (xa) do not increase linearly with a, but fulfill the 
relation 



{Xa) = {Xi) + Aal- 



(161) 
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As is shown in Fig. |^ Aai depends neither on the system size nor on the disorder. 

Apart from the deviation from the Gaussian form, P(ln g) possesses also the non-analyticity 
at In g = 0, shown in Fig. |50l The origin of this non-analyticity is the same as in the case of 
the critical distribution, Pc{g)- The contribution to the conductance comes mostly from the first 
channel and the chance to have g > 1 is negligibly small. 

More important is the question whether the variance is an unambiguous function of the mean 
value. This must be so if the one parameter scahng holds. The verification of the ambiguity 
is difficult in 2D and 3D. In ID, analytical calculations [119] showed that there is no universal 
relation between the mean and the variance of x. This was confirmed by numerical simulations 
[20, 120] of strongly disordered ID systems, indicating the existence of the second relevant length 
scale in the strongly localized regime. 

The second relevant length scale was indeed found in ID systems [170] and estimated ana- 
lytically by the formula 



is = —^-7TK, (162) 

sm Trp[rj) 

where p{E) is the density of states. Single parameter scaling is expected to be valid only when 
X > £s- This is not fulfilled in the band tail, where the density of states is small and, consequently, 
£s is large. The second relevant length scale was reported also in 2D systems in Ref. [121]. 

Figurel^shows that the variance, var In g is not a linear function of the mean value, (In g). 
This was observed already in [120] for the ID system. Very recently [121, 122], the non-linear 




Fig. 47. The probability distribution of the normalized difference 521 = X2 — xi for the 3D Anderson 
model (L — 10) compared with the Wigner distribution pi (dashed line) and Poisson distribution (solid 
line). Inset shows the same data in the linear scale [20]. 
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Fig. 48. Comparison of the conductance distribution for tfie quasi- Id weakly disordered system in the 
locaUzed regime and strongly disordered 3D system. While the distribution of p(ln g) is Gaussian in the 
quasi-ld system, for the 3D system, the distribution is asymmetric and narrower. The size of the quasi-ld 
system is 8 x 8 x and disorder W — 4. The size of the 3D system is 8^ and the strength of the disorder, 
W, is chosen to have the same value of the mean logarithm of the conductance [118]. 
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Fig. 49. Left: The differences Aai = {xa — xi) for the 3D Anderson model [118]. Right: The probability 
distribution p( a; a) for the strongly disordered cube (W = 20, L — 18). These data should be compared 
with those in Fig. 1731 which shows the distribution p{xa) for long weakly disordered quasi-ld systems. 
Note, the mean values, (xa) do not increase linearly. Also, the width of the distribution p{xa) decreases 
when a increases, similar to the distributions of parameters x at the critical point, shown in Fie. 1411 



relation between the mean value and the variance of In g, 

\ {-{\nq)f'^ d^2 ,,,,, 
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Fig. 50. The statistical distribution of the logarithm of the conductance in the strongly disordered squares 
(left) and cubes (right). The solid lines are Gaussian fits with the same mean and variance. Left: The 
2D Anderson model, W — 6, L = 200, (Ing) = —9.51, varln;? = 11.78. Note the sharp decrease 
of the distribution at Ing = 0. Right: The 3D Anderson model, W = 32, i = 18. (Ingr) = -18.88, 
varln g = 14.57 and = 14 and L = 10 ( (In g) = -6.42, varln g = 6.63). 
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Fig. 51. The variance, var In g as a function of the mean value, {Ing) for the 2D Anderson model. The 
expected linear relation, <160t . is valid only for rather small values of (In g) [118]. 



was numerically observed. 

Another test of the single parameter scaling was reported in Refs. [123, 124], where the 
electron wave function of large 2D disordered system was calculated numerically. The statistical 
distribution of the logarithm of the wave function at the distance r from the center of the lattice 
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was discussed. The numerical data was fitted to the distribution function 

, ^ 1 (Inl^'l +r/A)2 

g(-hi|vl/(rO|,r) = -==exp- ^ ' ' ' ' . (164) 
\/2t:(7 2.a 



The Gaussian form of the distribution H seems to be natural, since we expect that the wave 
function decreases exponentially at large distance. This decrease is controlled by the localization 
length, A. Surprisingly, fitting the numerical data to the distribution led to the r-dependent 
localization length, A(r). This result was interpreted as a failure of the scaling theory. 

However, an assumption that the distribution function H, given by Eq. M6A\ is exactly 
Gaussian, is not correct. Indeed, the Gaussian distribution possesses a long tail for both negative 
and positive values of its argument, but the value of — In |^'(r)| can not be negative if the wave 
function, is normalized. Therefore, the distribution H{x, r) must converge to zero when 
a; = — In l'I'(r) I 0''", in the same way as the distribution p{x), shown in Fig. |20| We believe 
that the fit of numerical data to the correct distribution function will recover the single parameter 
scaling. 



11.1 3D versus quasi-ld systems 

As discussed above, the statistical properties of the conductance in the 3D systems require nu- 
merical simulations. On the other hand, similar transition from the metallic to localized regime 
can be studied in weakly disordered quasi-ld systems. Indeed, such a system exhibits metallic 
behavior if the length of the system, Lz, fulfills the relations 

£<^Lz^Ni (165) 

where £ is the mean free path and N is the number of channels. The conductance of such a 
system is Nl/L^ ^ 1. By increasing L^, the conductance decreases. There is an interval of Lz, 
where g ^ 1. Further increase of the system length draws the system into the localized regime, 
where g ^ exp — 2L/A, where A = is the localization length. 

This scenario is very similar to the metal-insulator transition. Of course, we do not have a 
critical behavior in quasi-ld system (there is no divergence of the correlation length), but it is 
reasonable to expect that the conductance distributions, obtained in all three regimes, g ^ 1, 
5 ~ 1 and g <C 1, might mimic the main properties of the conductance distributions in the 
metallic, critical and localized regime. The advantage of quasi-ld systems is that they might be, 
at least within some approximations, solved analytically. 

The above idea was developed by Muttalib et al. [109, 110, 125]. By solving the DMPK 
equation in the regime of g ~ 1, they indeed found that the "critical" conductance distribution 
possesses non-analyticity close to g = 1 [110]. Also, in the localized regime, they observed that 
the distribution ^(In^) drops down at Ing = 1. Similar results were observed numerically in 
Ref. [126] where the transport in the quasi-ld systems with a corrugated surface was investigated 
and the conductance distribution in the "critical" regime was studied. 

However, these results provide us only with a qualitative description of the conductance 
distributions in the true critical regime and in the insulator. In the previous Section, we showed 
that the localized regime in the 3D system differs qualitatively from the localization in the weakly 
disordered quasi-ld systems. The origin of this difference lies in the different form of the spectra 
of parameters Xa- 
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Fig. 52. Comparison of the conductance distribution for tiie quasi- Id systems (squares) and cubic systems 
(circles) [118]. Left: The metallic regime: p{g) is Gaussian in both systems, but the widths of the distri- 
butions differ. This is consistent with our data for the universal conductance fluctuations, given in Table 
\\\ since var g of quasi-ld systems differs from that of 3D. Right: The critical conductance distribution for 
the 3D Anderson model compared with the distribution of the conductance for quasi-ld weakly disordered 
(W = 4) system of the size 8 x 8 x 210. Although the mean conductance is the same, the distributions 
differ in the region g > 1. The reason for this difference is that the difference, {x2 — xi) = (a;i) in the 
quasi-ld system, while it is much smaller in the 3D system, as shown in Fig. 1491 The non-analyticity of the 
distribution for the quasi-ld system was explained and qualitatively described in Ref. [110]. 



To demonstrate the difference between the two insulating regimes, we showed in Fig. I48lthe 
distribution p{\ng) for the 3D and the weakly disordered quasi-ld systems. Here, we compare 
in Fig. |52]the conductance distributions for both systems in the metallic and critical regime (Fig. 
I52> . For a quantitative description of the critical and localized regimes, we need to study the 
systems with strong disorder. 



12 Numerical scaling analysis 

Since the analytical calculations are possible only in the limit of the dimension d = 2 + e (e ^ 1), 
numerical simulations provides us at present with the most relevant information about the critical 
regime of the Anderson transition. 

The main problem of the numerical scaling analysis is the calculation of critical exponent, v 
in various physical models. Universality of the critical exponent for a given symmetry class and 
dimension is interpreted as evidence of the validity of the single parameter scaling. 

Historically, the first numerical simulations were performed for quasi-ld system [17, 127]. 
The scaUng of the smallest Lyapunov exponent was proved and the critical exponent, i' was 
found. Later, scaling of other variables was analyzed, mostly with the motivation to treat the 
true 3D systems and 2D systems with unitary and symplectic symmetry. We will discuss in this 
Section the scahng analysis of the conductance, conductance distribution and level statistics. 
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12.1 Scaling of the smallest Lyapunov exponent 

Pichard and Sarma [127] were the first who proposed to use the finite size scaling analysis for 
the calculation of the critical parameters of the Anderson model. They considered the quasi-ld 
system of the size L"^^^ x (d — 2 and 3) and calculate numerically the smallest Lyapunov 
exponent, 71. As discussed in Appendix lD.2l the Schrodinger equation for the Anderson model 
in the quasi-ld geometry can be written in the form 

^ M„ r !" V (166) 



where Af„ is the transfer matrix for the nth shce of the system, 

M„=(f"^" "J). (167) 

By multiplication of the transfer matrices, we obtain that the exponential decrease of the wave 
function along is given by eigenvalues A^ of the matrix 

n=l 

Oseledec's theorem states that in the limit of 00, all Lyapunov exponents, 7^ of the matrix 

Ml^ posses the Gaussian distribution with the variance proportional to the mean value. The 
smallest (in absolute value) Lyapunov exponent, 71, determines the localization length in the z 
direction. 

Since 71 oc L^, it is more convenient to use the parameter 

A = i^. (169) 

L71 

In this paper, we discuss the scaling behavior of parameter zi, 

2 2L71 

= A = -l7' ^'''^ 

which is closely connected to parameter xi, used in the previous Section for parametrization of 
the contribution to the conductance. From Oseledec's theorem it follows that zi converges in the 
limit of 00 to the mean value, and var zi oc L/L^- Consequently, the difference between 

— 1/2 

the numerically calculated value of zi and the mean value is cx ' so that zi is free from any 
statistical problems provided that the system is sufficiently long. Of course, zi is a function of 
the disorder, W, energy, E and the system width, L. 

For a sufficiently large width L of the system, we can estimate the disorder and L dependence 
of zi from basic physical considerations. It is reasonable to expect that when the system is in 
the localized regime, W ^ Wc, then zi converges to the ratio 2L/X. On the other hand, in the 
metallic regime, W <C Wc, we have that the spectrum of Za is linear, as given by Eq. j320l l. 
of Appendix lC.il Therefore, zi = where N is the number of channels. Since oc in 
the case of the 3D Anderson model we obtain zi ~ in the metalhc regime. Finally, at the 
critical point, we expect that zi = const. 
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The above consideration can be summarized as follows: 

zi = 27i- <^ const W = W^ (171) 

[ L/X W > W^. 

Relation MIW enables us to identify the critical point from the numerical data. The left panel of 
Fig. |53]presents the numerical data for the quasi-ld Anderson model x L^. We indeed see 
that for disorder W > Wc, z\ increases when L increases, while for W < Wc, z\ decreases with 
L. From the L dependence of z\, we estimate approximately the value of the critical disorder, 
Wc ~ 16.5 and the critical value, zic = zi(W ^ Wc) ~ 3.44. 

The next step in the scaling analysis was done by MacKinnon and Kramer in Ref. [17]. They 
assumed that zi is a function of only one parameter, 

ziiL,W)^FiL/aW)) (172) 

where £,{W) is the correlation length. Relation illlt follows from the assumption of the validity 
of single parameter scaling. The right panel of Fig. l53lshows that indeed all the data zi{L, W) 
lie on the universal curve. 
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Fig. 53. Left: The L-dependence of the variable zi — 2Lji/Lz for the 3D Anderson model. The data was 
calculated for the quasi-ld systems, x Lz and for disorder increasing from W — 16 (bottom) to = 17 
(top) with step AW — 0.1. The length Lz is sufficiently long to guarantee the relative accuracy of 0.2% 
for small L, and of 1% for L = 22. The dashed line is the fit of the critical value zic ~ 3.44 [128]. Right: 
The same data plotted as a function of L/^{W), Eq. <172t . Note the logarithmic scale on the horizontal 
axis. The function F has two branches, the lower branch for the metallic regime, and upper one for the 
insulating regime. The data with L = 4 were excluded from the data sets, because they do not lie on the 
universal curve due to the finite size effects. Inset shows divergence of the correlation length, ^{W), at the 
critical point. Note, ^ is calculated up to a multiplicative factor. 
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Since zic = zi {W = Wc) does not depend on the system size, L, the correlation length must 
diverge at the critical point, 

^{W) (X \W - Wcl^" . (173) 

Figure |53 confirms that zi (x {W — Wc) in the critical region. Therefore, if we expand the 
function F{x) in the Taylor series, and keep only the first two terms, F{x) = F{0) + Ax", 
we have that exponent a = l/iy. Consequently, zi{L, W) is given in the critical region by the 
simple scaling equation, 

zi{L,W) = zic + AiW-Wc)L^/r (174) 

The fit of numerical data to Eq. ( I174> enables us to calculate both the critical exponent, v, and 
critical disorder, Wc- 

The easiest scaling analysis can be performed in the following two steps [129]. First, we can 
calculate the linear fit 

ziiL,W) = z^°\L) + Wz[^\L) (175) 
Comparing the rh.s. of Eq. (I175> with Eq. J174> . we have 

z[°^ = zie - AWcL^/" (176) 

and 

z[^\l) = AL^^" . (177) 

Now, we can use Eq. (^^} to calculate the critical exponent. Next, with the use of Eq. illll we 
can write Eq. ( I176> in the form 

zf^ =zic~Wcz['\ (178) 

so that the slope of the linear fit z[^'' vs. z[^^ determines critical disorder, Wc- 

The physical meaning of the correlation length, £,{W) can be estimated by comparing Eq. 
(I172t with Eq. ( I17U . Clearly, ^(W^) = A in the insulating side of the transition. Also, we find 
that F{x) a; in the insulator and F{x) ^ x^^ in the metal. Then, from the expression of 
zi — 2L/{N£) — 2{g)^^ = l/{2La) (a is the conductance) we find that in the metallic regime 
the correlation length ^(VF) cx [17]. 

The correlation length (,{W) was first calculated numerically in Ref. [6]. The critical disor- 
der, Wc ~ 16.5, and the critical exponent, v « 1.50, were calculated. These calculations were 
repeated by many other authors [130-132] with the use of various scaling analysis, and for in- 
creasing system size. Surprisingly, these new analysis did not bring any considerable corrections 
to the critical parameters, obtained in the pioneering work, [17]. At present, the most accurate 
estimation of the critical exponent is 

1/ = 1.57 ±0.02, (179) 

obtained by very detailed scaling analysis of numerical data for zi in Ref. [133]. In Ref. [52], the 
phase diagram in the energy-disorder phase space (shown schematically in Fig.|6} was calculated 
for the first time for various distribution of random energies. 
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Fig. 54. Left: The same data as in left panel of Fig. I53l but plotted as a function of the disorder W. Solid 
lines are linear fits, given by Eq. J 1741 for 4 ^ L <C 22. In the ideal case, all lines have a common 
crossing point. This never happens for numerical data, as shown in the detailed plot in the right panel. The 
deviations from the universal scaling relation, given byEq. (T74j are either due to the finite size effects, 
discussed in Sect. 112.21 or due to the insufficient accuracy of numerical data. 

The scaling analysis [17] solved also the problem of the existence of the Anderson transition 
in the 2D orthogonal systems. It was shown, that there is no metallic phase, but the localization 
length, given by ^{W), is extremely large for small disorder (of orders of 10® lattice sites for 
W — 1. This result also explains why various previous works identify the Anderson transition 
also in 2D systems: these works analyzed only small systems [134]. The numerical data for 
the correlation length, ^{W) of the 2D Anderson model are given in Refs. [6,45] and for 3D 
Anderson model in Ref . [6] . 



where 77, are further L-dependent parameters. They determine how zi depends on various mi- 
croscopic parameters of the model, for instance distribution of the disorder, correlation length of 
the disorder, magnetic field. The one-parameter scaling is valid only if in the limit of L ^ cx3 all 
T]i 0. Only under this assumption, the equation ( I173> can be recovered, with (i = L/£^{W). 
Thus, the assumption of the single parameter scaling requires that 



12.2 Finite-size corrections 



A more general formulation of Eq. ( I173> is 



Zl = F{Cl,TJl,7j2, ■ ■ .), 



(180) 



(181) 



where jji are irrelevant scaling exponents. 

Although these parameters play no role for sufficiently large systems, they might influence 
the scaling analysis for smaller L. For instance, the one parameter scaling, given by Eq. ( I174> 
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Fig. 55. The critical disorder, Wc, and the critical exponent, u (inset) calculated from the finite size scaling 
analysis of the first 9 Lyapunov exponents [128]. The numerical data for 4 < L < 22 were used. The 
critical disorder, Wc, converges to the "correct" value, Wc ~ 16.5 when numerical data for smaller system 
size, L < Lmin are not considered and Lmin increases. The finite size effects are larger for the higher 
Lyapunov exponents, and may lead to wrong estimation of the critical parameters, when too small systems 
are considered. On the other hand, increase of Lmin influences the accuracy of the numerical fits, especially 
in the case of the power fit <177> for the critical exponent. Note, the finite size effects are very small in the 
scaling analysis of the smallest Lyapunov exponent. 



requires that all linear lines in Fig. |33]cross in one common point for W ~ Wc and zi = zic. 
This is evidently not true for the smallest L = 4. 

From Eq. ( I180> it follows that when L is not sufficiently large, one can generalize the scaUng 
analysis by the inclusion of an additional term on the r.h.s. of Eq. J173> . In the most simple case, 
when only one irrelevant parameter is considered, we obtain 



W) = zic + A{W - Wc)L^''' + BL- 



(182) 



[135]. We can estimate the critical parameters by fitting the numerical data to the function 
(I182> . A more detailed scaling analysis might consider also higher order terms in powers of 
Ci = ~ Wc)L^^'^ ■ The most detailed scaling analysis was performed in Ref. [133], where 
the numerical data were fitted to a function of 11 parameters. 

Another possibility to eliminate the finite size scaling effects is to perform the scaling analysis 
with reduced data sets by omitting data for the smallest system size [54, 128]. Using only the 
numerical data for L > Lmin, one can study the Lmin dependence of the critical parameters. We 
demonstrate this method in Fig. |55l 
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Fig. 56. Estimation of the critical exponent, u, from the power fit, given by Eq. <177L Data for the smallest 
Lyapunov exponent, zi, as well as for Za, a — 2 ~ 5 and 9 are shown. We remind the reader that zi^' is the 
slope of the W dependence of Za, defined by Eq. <176> . 



As discussed in Sect. ^| the critical conductance distribution is not a self averaged quantity. 
Therefore, the scaling theory might, at least in principle, contain an infinite number of inde- 
pendent relevant parameters. To verify this possibility, it is worth to verify whether the higher 
Lyapunov exponents, Za, scale in the same way as the zi. This was done in Ref. [128] and 
reproduced in [112]. 

As shown in Fig. |56l aheady the most trivial scaling procedure, given by Eqs. (I176I178> . 
provides us with reliable evidence that the nine smallest Lyapunov exponents indeed scale with 
the same critical exponent, v. This proves that all Lyapunov exponents can be expressed as a 
functions of only one variable, L/£^{W), with the same correlation length, ^{W). 

Scaling of the higher Lyapunov exponent serves also as a nice example how the finite size 
effects influence the scaling analysis. As was shown in Fig. |55l the finite size corrections are 
almost negligible in the case of zi, but they increase when a increases. The reason lies in the 
properties of the spectra of parameters Za, or, equivalently, of Xa- In Section^] we found that 
for a given size of the system, L, only parameters Xa with a < L are size independent. The rest 
of the spectra depends on L. The same must hold for Za- Therefore, only the data with L > a 
are relevant for the scaling analysis of the higher Lyapunov exponent. In the case of zg, only the 
data for i < 12 are relevant for the scaling analysis. 



Numerical analysis of the quasi- Id systems provides us with the most accurate estimation of the 
critical parameters. Nevertheless, the use of quasi-ld geometry is rather artificial. It would be 
more suitable to prove the scaling of the conductance of the true d dimensional systems. This 



12.3 Scaling of higher Lyapunov exponents 
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problem is rather complicated since, as we have seen in Sect. QO] the conductance is not self- 
averaged quantity in the critical region. We need therefore to calculate the mean value, (g) from 
the statistical ensemble of A^stat different cubes. 



\9) 



— H 5i 



(183) 



where gi is the conductance calculated for the ilh sample, and to estimate the accuracy of the 
mean value by using the relation 



acc g = 



var g 



(184) 



Since both (g) and var g are of order of unity, we need A'stat ~ 10^ to reach the relative accuracy 
-0.1%. 

The scaling formulas, presented in previous Sections, are valid for the scaling behavior of the 
mean conductance as well. Thus, we expect that the disorder and system size dependence of the 
mean conductance in the critical region is given by 



{g) = {g)c + A{W-Wc)L 



(185) 



where (g)c is the critical conductance. A similar equation can be constructed for (In g). 

For the 3D Anderson model, the scaling of the mean conductance, (5), and of the typical con- 
ductance, exp(ln(7) was numerically confirmed in Ref. [136]. The calculated critical exponent, 
V — 1.57, agrees with the result of the scaling analysis of the smallest Lyapunov exponent [133]. 
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Fig. 57. Left: The L dependence of the mean conductance of the 2D Ando model. For disorder W < Wc, 
(g) increases with L, and for W > Wc, (g) decreases with L, in agreement with scaling hypothesis. Right 
panel shows the same data plotted as a function of only one parameter, L/(,{W) with ^ oc \W ~ Wc\^'^ . 
Appropriate choice of the correlation length, £,{W), for each value of the disorder, W enables us to scale 
all numerical data to one universal curve. Note the logarithmic scale on the horizontal axis. 
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Fig. 58. Definition of percentiles. Note the logarithmic scale on the vertical axis. 



As an example of the scaling analysis of the conductance, we present the most recent numer- 
ical data for the mean conductance of the 2D Ando model [54]. The left panel of Fig. l57lshows 
the L-dependence of the mean conductance for fixed disorder. The metallic, localized and criti- 
cal regime can be estimated in the same way as for zi, with the only difference, that increasing 
(g) indicates the metallic phase in this case. The right panel of Fig. 1571 shows the same data 
but as a function of L/£^. The data confirms that indeed (g) ~ g{L/£_) is a function of only one 
variable. Two branches of the function g correspond to two different transport regimes, metallic 
and localized. 

The numerical proof of the scaling of the mean value is still not sufficient for the verification 
of whether or not one parameter scaling really works. Namely, we cannot exclude that the higher 
cummulants of the conductance do not scale. Of course, it is impossible to verify the scaling of 
all cummulants. Instead, we discuss in the next Section the scaling of conductance distribution. 



In the previous Section we verified the scaling of the mean conductance, (g) . However, we still 
cannot exclude the possibility that the length and disorder dependence probability distribution, 
p{g), is determined by an infinite number of parameters, for instance all higher cummulants of 
the conductance. To prove that the entire conductance distribution scales as a function of only 
one parameter, the scaling of the percentiles, was analyzed numerically in Ref. [137]. 
The percentile, gq, is defined by the relation 



By definition (I186> . the probability that g < gq, equals to q (Fig. I58> . 

In Ref. [137], the one parameter scaling of percentiles, gq, was proved for four values of q, 
q — 0.025, 0.17, 0.50 (median) and 0.83. All four variables obey the one parameter scaUng with 
the critical disorder Wc close to 16.5 and with the critical exponent 1.56 < v < 1.60. 

Suppose now that two percentiles, ga and gp, {a < f3) obey the single parameter scaling. 
Then, the percentile g^ [a < < f3) must scale, too. Also, if ga and gf^ scale, then the difference 



pig)- 



12.5 Scaling of the conductance distribution 
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9(3 — da scales. Therefore, for the proof of the one parameter scahng, it is sufficient to prove the 
scahng of only a few percentiles, which was done in Ref. [137]. We conclude that the numerical 
verification of the single parameter scaling of a few percentiles provides us with the evidence that 
the entire probability distribution, p{g), obeys the single parameter scaling in the critical region. 

12.6 Scaling of the level statistics 

As discussed in Sect. 15.21 the distribution p{s) of the differences between the neighboring 
eigenenergies depends on whether the system is in the metallic, localized or critical regime. 
In Refs. [58, 138, 139], the scaling analysis of the level statistics was proposed and studied. 

The critical parameters were calculated for the 2D Ando model [100], the 3D Anderson 
model [61, 140], and for the problem of quantum percolation in 3D system [141]. Recently, the 
scaling analysis of the level statistics was applied to the symplectic model on the fractal lattice 
with the aim to prove that the lower critical dimension for the symplectic systems is less than 
2 [142]. 

12.7 Scaling of the inverse participation ratios. 

The scaling of inverse participation ratios, Iq, defined in Sect. |5] was performed recently in 
Ref. [56]. Contrary to the conductance, which is defined only for energies inside the unperturbed 
energy band, \E\ < 6V, Iq can be calculated for the entire energy spectrum of the Hamiltonian, 
and can be used for the verification of the universality of the metal-insulator transition along the 
critical line. However, one has to keep in mind that Iq is not size-invariant at the critical point 
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Fig. 59. The fractal dimension, dq, obtained from the scaling analysis of the inverse participation ratios, Iq, 
for the 3D Anderson model with the Gaussian and box disorder. The dashed line is {q — l)d for d = 3. 
Inset shows the estimated critical exponent for three critical points. At the band center with the Gaussian 
(circles) and box (boxes) distribution, and for the Gaussian disorder W — 2 with critical energy Ec ~ 6.58 
(triangles). The electron eigenenergies and eigenfunctions were calculated numerically for cubes of the size 
8 < L < 54 [56]. 
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but behaves as 

as discussed in Sect. |5] 

In the scaUng analysis, the logarithm of Iq was used, since the values of Iq{En) might fluc- 
tuate in orders of magnitude within a small energy interval, SE (Figs. IIOIII It . The quantity of 
interest is then 

iq{L,W)^{{\nIq{E,-,))sE), (188) 

where the averaging is performed within the energy interval, 6E, and over statistical ensemble 
of microscopically different samples. Iq{L, W) is then fit to the scaling equation, 

iq{L,W) = A - {q - l)dqliiL + B{W -Wc)L^/'' (189) 

and critical parameters, Wc, dq and are calculated for the 3D Anderson model with the Gaus- 
sian and box distribution of random energies. 

Figure|59]shows that not only the critical exponent, v, but also the fractal dimensions, dq, are 
universal, independent of the microscopic details of the model and on the position of the critical 
point on the critical line. The fractal structure of the critical wave function was studied also in 
Refs. [143-146]. 



13 Scaling in the d-dimensional systems 

The numerical scaling analysis provides us with rather accurate estimation of the critical ex- 
ponent for the 3D Anderson model. However, the obtained results are in disagreement with 
expectations of the theory, which reports 1^=1 for d = 3 [147]. It is therefore important, for 
the detailed comparison of the theory and numerical data, to calculate the dimension dependence 
of the critical exponent, and check, whether or not agreement with theory is better for, d — ^ 2+ 
or for d > 3. We summarize here the very recent data for the critical exponent, calculated in 
Ref. [113] for 2 < c? < 4, and we present also new data for d — 5. 



13.1 Dimension d = 2 + e 

Figure|60|shows the disorder dependence (d = 2 + e) of the critical exponent, ly, calculated from 
the finite size scaling of the smallest Lyapunov exponent of quasi- 1 d bifractal lattices with fractal 
structure of the cross-section, shown in Fig. |^ These data are compared with the theoretical 
prediction, based on the analytical calculation of the function f3{g). 

In the limit of e <C 1, the critical disorder Wc ^ e ^ 1 and the critical conductance gc ^ 
^ 1. The function f3{g) can be expanded in power series of g^^. It is more convenient to 
use, instead of the conductance g, the parameter 




The size dependence of the parameter t is given by the equation 
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Fig. 60. The critical exponent, v, as a function of (e = d — 2). Shown are the data for three bifractals, 
discussed in Sect. ll0.2l (Fig. l45l . and for a cubic 3D system. The dashed line is the linear fit, 1.24+0.365/e. 
The dotted line is the analytical e - expansion of the critical exponent, given by Eq. <194> . Note, there is 
no agreement between the theoretical prediction (dashed line) and the numerical data. Note also that two 
systems with the same spectral dimension have the same critical exponent, as expected. 



and the critical exponent, ly, is then given by the relation 

-1 dm 



(192) 



dt 

For the orthogonal systems, the t-expansion of the function /3(i) reads [11] 

27 

Po{t) =et^ 2e - 12C(3)i5 + yC(4)i' + • ■ • ■ (193) 

In Eq. ([T93j, C(3) = 1.202 and C(4) = 7rV90. 

Since the expansion M931 is known only up to the 6th power in t, it is difficult to estimate 
the accuracy of the obtained results for the critical conductance and critical exponent, specially 
when e is not small. For instance, for the 3D system, (e — 1) one finds, by solving the equation 
Poitc) = 0, that tc — 0.395. This agrees qualitatively with the estimation of critical conduc- 
tance, gc — l/(27rtc) ~ 0.40 (we remind the reader that the numerically observed values of the 
critical conductance are 0.445 and 0.280 for periodic and hard wall boundary conditions, respec- 
tively). However, from Eq. ( I192I I we obtain i' = 0.67, which is far from the numerical result, 
v w 1.57. The agreement with the numerical data is not better for small e, as is shown in Fig. |60| 
which compares the critical exponent, calculated from the e expansion, 

^C(3)e' (194) 
e 4 

with our numerical data. Clearly, there is no agreement between the theory and results of numer- 
ical simulations. 

For completeness, we add the e expansion of the /3 function for symplectic systems. It can be 
obtained from expression ( fT93l with the use of the symmetry relation [97] 

/3s = -2/3o(-i/2), (195) 
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which gives 

4 d4 

Note, /3s (i) is positive for e — > which confirms the existence of the critical point in the 2D 
symplectic systems. 

13.2 Dimension d> 3. 

The numerical data for the critical exponent in higher dimensions does not agree with the self- 
consistent theory [147] which predicts that 

For instance, for the 3D systems, Eq. ( I197> predicts v = 1, which clearly disagrees with the 
numerical result, v — 1.57. 

In order to get insight into the dimension dependence of the critical exponents, the finite size 
scaling analysis for d = 4 was performed in Refs. [113, 116, 148], and for c? = 5 in the present 
paper. Figure|^presents numerical data for zi obtained for the quasi-ld systems x Lz and 
i** X Lz- Although only the data for small L can be calculated, obtained results confirm the 




Fig. 61. The first Lyapunov exponent, zi, as a function of disorder, W, for various L for the 4D (left) and 
5D (right) Anderson model. The critical parameters are Wc = 34.3 ± 0.2, u = 1.12 ± 0.05 for 4D and 
Wc = 57.26 ± 0.2, V = 0.93 ± 0.05 for 5D. Inset shows the L-dependence of zj^', defined by Eq. (TTsl . 
Power fit z^^^ vs. L, given by Eq. <177> determines the critical exponent, u. 
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existence of the critical points in both systems. The simple scaling analysis, given by Eqs. illSl 
I178t was used to estimate the position of the critical points and of the critical exponents. For the 
4D Anderson model, we find 

i^4D = 1.12 ±0.05. (198) 

This result agrees with previously obtained data [1 13, 1 16]. To the best of our knowledge, there 
have been no published data for the critical exponent of the 5D Anderson model yet. Our esti- 
mation of the critical exponent is 

z^5D = 0.94 ± 0.05. (199) 

Figure|62]summarizes all obtained numerical data for the critical exponent of the orthogonal 
Anderson model in d dimensions and compares them with predictions of analytical calculations. 
We conclude that there is no agreement of numerical data with the theory, neither for small 
dimension, nor for d > 3. 



13.3 Theory vs (numerical) experiment 

Disagreement between the theoretical predictions and numerical data might lead to the conclu- 
sion that the numerical scaling analysis is insufficient or even wrong [149-151]. It is not our 
aim to discuss these objections here (see, for instance, the comment to Ref. [151], published in 
Ref. [152]). We only concentrate on the discussion of advantages and disadvantages of numerical 
analysis. 




Fig. 62. Dimension dependence of the critical exponent. Numerical data differ considerably from the 
theoretical predictions, both for d — 2 <c; 1 and for d ^ 2. The dashed line is the e-expansion, given by 
Eq. <194t . Triangle is The estimation of Hikami, obtained by the Fade approximation of the e expansion of 
the /3 function, /^Hikami ~ 0.73. The solid line is the analytical prediction <197t 
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The main objection against the numerical methods is that they are always restricted to the 
systems of finite size. This limitation can be partially avoided by the scaling analysis, Eq. ( I173> . 
and by including the irrelevant scaling variables, Eq. (I182t . Also, the accuracy of obtained 
critical parameters can be estimated by elimination of numerical data for small system size, as 
discussed in Sect. 112.21 Fortunately, the finite size effects seem to play a negligible role in higher 
dimension. The first estimation of the critical exponent for the 3D Anderson model, ly w 1.5, 
obtained in pioneering work of MacKinnon and Kramer, [17], differs only in a few per cent from 
the best today's estimation, ly = 1.57 [136]. Also, for the 4D Anderson model, the estimation 
ly = 1.1 ± 0.1 was obtained already in Ref. [116], where really small lattices, typically of the 
size 3'^ X Lz were analyzed. Nevertheless, the analysis of larger samples, up to 7^ x L^, [113] 
brought no corrections to the critical exponent. As is shown in Fig. |^ the additional data for 
even large systems have no influence to the critical exponent, originally estimated for smaller 
systems. We conclude therefore that the numerical simulation provides us with reliable data in 
higher dimension. 

Our belief that the numerical data for the critical exponent are correct, is supported also by 
results of the scaling analysis of various groups. In the 3D systems, various models were stud- 
ied, isotropic and anisotropic [153], with diagonal or off-diagonal disorder [131, 132]. Scaling of 
various parameters was analyzed, inclusive the smallest Lyapunov exponent, higher Lyapunov 
exponents, conductance [136, 137], conductivity [154], level statistics [61, 141] and inverse par- 
ticipation ratio [56]. All these works report the critical exponent close to 1.5, with the accuracy 
which definitely excludes the possibility that 1^=1. 

As discussed in Sect. O finite size effects become stronger in lower dimension. The reason 
is that the mean free path, £ is larger (Fig. I28> . A typical problem caused by the finite size of the 
system is shown in Fig. I^which presents the latest numerical data for the mean conductance, 
(g) of the 2D Anderson model. Although we accept that there is no Anderson transition in d = 2, 
the numerical data seem to mimic the metallic behavior. For W — 1, the mean conductance (g) 
increases when L increases. One might argue that there is a critical disorder, Wc ~ 2. However, 
as discussed already in Sect. 18.61 the above conclusion is not correct. To determine the character 
of the transport regime, the size dependence of the conductance must be carefully analyzed. The 
increase of the mean conductance for = 1 is the manifestation of the ballistic regime. Indeed, 
the mean free path £ Ri 17 for W = 1 is comparable with the size of the system when L < 100. 
Also, the variance, var g, is much larger than expected universal value, typical for the metallic 
regime (Fig.l?ll. 

Also, the decrease of the mean conductance with the system size for disorder W = 3 and 
W = A cannot be associated with localization. Indeed, the same data, plotted in Fig. |25]in the 
logarithmic scale, show that the decrease of the conductance is not exponential, but logarithmic, 
and is caused by the weak localization. To see the exponential localization, one needs much 
larger system size. 

The big advantage of the numerical simulation is that it can relatively easy analyze statistical 
properties of any quantity of interest. No averaging is necessary in the course of calculations. 
All mean values can be calculated "from first principles". This cannot be done analytically. 
The analytical theory must solve the problem how to perform the average over the disorder 
Wrong averaging might lead immediately to wrong results [151]. In our opinion, the discrepancy 
between the numerical data and results of analytical theories is due to the inability of analytical 
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theories to analyze completely the statistical fluctuations in the critical regime. 



14 Two dimensional critical regimes 

The critical regime in the 2D models deserves a special attention. As discussed above, only 
systems with symplectic symmetry exhibits the metal insulator transition in dimension d — 2. 
The 2D systems with unitary symmetry posses, in the presence of strong magnetic field, the 
critical energies, Ec where the localization length diverges. 

On first sight it seems that the 2D systems can be easier simulated numerically since the 
lower dimension of the system allows one to calculate the conductance for much larger samples. 
However, this advantage is "compensated" by much stronger finite size effects. 

Besides the calculation of the critical exponents and critical conductance distribution, the 2D 
critical regime is suitable for the verification of the general relation between the conductance and 
conductivity, 

(g) = <J, (200) 

given by Eq. ( I109t . Contrary to the orthogonal 2D systems, where Eq. ( I200t holds only in the 
diffusive regime, i.e. when the size of the system is smaller than the localization length, L ^ X, 
Eq. ( I200t holds also in the limit of infinite system size at the critical point of the unitary and 
symplectic models. 

Note, equation ( I200> compares two different quantities. The conductance, g, is given, by 
definition, by the transmission properties of the disordered sample at zero temperature. The 
value of the conductance depends on the actual distribution of disorder inside the sample. Owing 
to the quantum character of electron propagation, the conductance is not a self-averaged quantity. 

On the other hand, the conductivity, a, is a material parameter. It characterizes the transport 
properties of an infinite system. When calculated for the system of finite size, L, ^ a fluc- 
tuates around the mean value, but the fluctuations decrease when L increases. Contrary to the 
conductance, the conductivity, a is a self-averaged quantity. 

For the critical 2D regimes, Falko and Efetov [158] derived the following relation between 
the fractal dimensions of the critical wave function, dq, and the critical conductivity, a, derived 
the relation 



_q_ 

f3Tra{h/e^)' 



dq = 2- ,^ (201) 



Here, (3—1,2 and 4 determines the symmetry of the system^. Eq. ( 120 11 1 holds for small q, when 
terms proportional to higher powers of q can be neglected. Since both the wave functions and 
the conductance can be calculated numerically in the critical point, we can verify the relation 
(120 H by direct numerical simulations. It will be shown in next two Sections that indeed both 
Eqs. ( I200> and ( 120 H are satisfied within the accuracy of numerical data. 



'Numerical algorithm for the calculation of the conductance is described in Ref. [154, 155] and in Ref. [157] for the 
case of critical quantum Hall regime. 

*Note, that /3 = 1/2, 1 and 2 is used in original papers, which adds an additional factor of 2 in the denominator on 
the rh.s. of Eq. 12011 . Also, note that factor of 2 for two orientation of electron spin is not included in the definition of 
the conductivity, a. 
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14.1 Symplectic models 

The difficulty of the analysis of the 2D symplectic models is manifested by a wide variety of 
values of the critical exponent, 2 < ly < 2.88, reported in the literature within the last 15 
years [47,49,99, 100]. These discrepancies are due to the strong finite size effects. Recently, 
the finite size scaling of the smallest Lyapunov exponent on the SU(2) model, [49] provided the 
following estimate of the critical exponent. 



V = 2.75 ±0.01. 



(202) 



This value can be considered as the most accurate estimation of the critical exponent. The analy- 
sis of the scaling of the mean conductance for the 2D Ando model, [54] led to the similar value. 



2.80 ±0.04. 



(203) 



The scaling behavior of the mean conductance, ( g), is shown in Fig. |^ The critical conductance 
distribution for the 2D Ando model is shown in Fig. |^ Since the mean conductance. 



(3)c « 0.71, 



(204) 



is close to 1, the distribution possesses all characteristic properties of the critical distribution for 
the 3D Anderson model (Fig. I43> . Comparison of the critical distribution, calculated for the 
periodic and hard wall boundary conditions is shown in Fig. |^ 

The numerical data for the conductance, together with Eq. ( I200> enables us to verify the 
relation (120 1> between the conductance and the fractal dimensions, dg. To do so, the wave 
function of the 2D Ando model was calculated numerically for the disordered sample of size 
260 X 260. Then, the sample was divided into the small squares fta of size Lq x Lq, and the 




Fig. 63. The fractal dimension, dq, as a function of q, calculated for the 2D Ando model. For q < 1, the 
fractal dimension, dq, is a linear function of g, in agreement with Eq. <20H . The solid line is the linear fit 
with the slope 0.11202 which determines the critical conductance, CTc = 0.71, and which agrees with the 
estimation of the critical conductivity, {g}c = 0.71 [159]. 
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quantity 

P,iLo)^y^< V |*„(^ir y (205) 



was calculated for the iVstat = 10 different realizations of the disorder. The wave function '^n 
is the eigenfunction of the disordered Ando Hamiltonian, corresponding to the eigenvalue En, 
closest to the band center, E — 0. Since the wave function is expected to posses the multifractal 
spatial structure, the -dependence of Pq{Lo) is determined by the fractal dimensions, dq, 

Pq{Lo) ^ L-^'''^^\ (206) 

Numerically calculated fractal dimensions, dq are plotted in Fig. |63las a function of q. The results 
confirm that dq decreases linearly with q, for q < 1. The slope, given by Eq. ( 120 1> . determines 
the critical conductivity, Uc — 0-71e^/h, which is exactly equal to the critical conductance, {g)c- 

14.2 Critical quantum Hall regime 

The 2D disordered system in a strong magnetic field possesses, inside each Landau band, the 
critical energy, Ec shown schematically in the left panel of Fig. |63 When the Fermi energy 
crosses the energy Ec, the transmission from one Hall plateau to another one appears [31, 160]. 
The existence of the critical energy, Ec, was numerically proved in Refs. [50]. The correlation 
length, ^, diverges on both sides of the critical energy as 

aE)^\E-Ec\-'', (207) 

with the critical index, v w 2.33 [32,50]. 

The critical quantum Hall regime can be studied numerically within the Hamiltonian, ( I19> 
with the Peierls phase, (p = B{eaf' /Tl. Since we are interested only in the critical parameters, the 
periodic boundary conditions are used in the direction perpendicular to the propagation. Then, 
the Peierls phase, (p, must fulfill the relation (pL = 27r7i, where L is the size of the lattice in the 
transversal direction, and n is an integer. 

The numerical analysis of the critical regime is difficult, since the disorder W must be small 
in order to keep the Landau levels separated from each other Then, however, the mean free path, 
£ is large. This problem can be avoided with he use of the spatially correlated disorder [29]. 

In Ref. [159], the sample averaged conductance, (g), and the conductivity, a, were calculated 
for the first and the second Landau band. To eliminate the finite size effects, the random disorder 
was spatially correlated. Numerical data was fitted to the formulas 



y 

L 



{g{L)) =gc-g^[^\ (208) 



and 



L 

where y y' are irrelevant scaling exponents 



a{L) =ac-ao[-^] . (209) 
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Fig. 64. Left: The density of states of the weakly disordered 2D system in a strong magnetic field. The 
energy spectrum consists of the Landau bands. Inside each Landau band there is a critical energy, Ec, 
where the localization length diverges. Note, the position of the critical energy depends on the strength of 
the disorder. Right: The critical conductance distribution for the first Landau band for (fi = 2n/8a. The 
disorder is = 1.4 and the critical energy, Ec — —3.29. 



Numerical analysis, performed for various strength of the disorder proved that both the mean 
conductivity and conductance converge to the same values listed in Table|2] This data proves the 
relation ( I200t for the critical quantum Hall regime and they also shows that the critical conduc- 
tance (or conductivity) is universal, independent of the microscopic model and on the Landau 
level. 

Critical values listed in Table |2] are significantly larger than the commonly accepted value, 
0.5e^//i, found in Ref. [161], and confirmed in previous numerical simulations [162]. However, 
they are in agreement with numerical data obtained on the Chalker Coddington model [163, 
164]. Critical conductance, (Tc ~ Q.61e^/h, is also consistent with the recently calculated fractal 
dimension, di = 1.739, and with the relation ( 120 It . 

1 e2 

(210) 



27r(2 - di) h 



The difference between the theoretical and numerical data is probably due to the spatial inhomo- 
geneity of the electron distribution, which was neglected in the analytical calculations. 



critical conductivity 


CTc 


0.58 ±0.02 


critical conductance 






1st. Landau band 


9cl 


0.60 ±0.02 


2nd Landau band 


9c2 


0.61 ±0.03 



Tab. 2. The critical conductivity, CTc, and the critical conductance, {g), (in units of /h) for the two lowest 
Landau bands [159]. The data was obtained by fit of the numerical data to Eqs. J208> and <209t with 

y^y' ~ -0.4 
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15 Possible theoretical description of the localized regime 

In this Section, we present two possible descriptions of the locahzed regime. Both of them are 
based on the generahzation of the theoretical methods developed for the analysis of the transport 
in diffusive regime. In Sect. 115.11 we discuss the possibility to generalize the DMPK equation, 
[165] and in Sect. 115.21 a simple generalization of random matrix theory is proposed [91]. 



15.1 Generalized DMPK equation 



In AppendixIBIwe present the DMPK equation, which describes successfully the transport prop- 
erties of weakly disordered quasi-ld systems. The theory contains only one parameter - the mean 
free path £, which measures the strength of the disorder. 

The DMPK equation was derived under the assumption of the homogeneity of the matrices 
u and V which parametrize the transfer matrix (see Appendix l A. 3 I for details). Physically, homo- 
geneity of matrices u and v means homogeneous distribution of the electron on the opposite side 
of the sample. This is possible only if the electron has many paths to travel from one side of the 
sample to another side. Mathematically, this requirement is expressed by Eq. ( 127 6> . which can 
be written in the form 




\Ubc 



l + Sab 

N + 1 



(211) 



Here, N is the number of channels. Clearly, this assumption is fulfilled only in the limit of 
weak disorder, when the electron can choose, on its travel through the sample, many equivalent 
paths. This is not true when the disorder is strong and the electron hardly finds a single trajectory 
propagating from one side of the sample to opposite one (Fig. l7H . 

Recently, Muttalib and Klauder [165] proposed the generalization of the DMPK equation. 
Without any restriction to the value of the matrix elements Kab, they generalized the DMPK 
equation into the form 



d{L,/l) J ^ dXa 



1 X - O 

7 2^ TYT 



oXa 



with the Jacobian 



N 



J=Y[\Xa-Xb\''-\ 



"fab 



a<b 



2Ka 



(212) 



(213) 



In Eq. ( I212> . the symmetry parameter (3 — 1. 

For the weak disorder, one can substitute for Kab from Eq. ( 121 1> and obtain that Eq. M\2\ 
reduces to the "classical" DMPK equation. For strong disorder, the parameters Kaa and ^ab 
represent the additional free parameters which must be estimated from numerical experiment. 

To estimate values of ^ab in the locahzed regime (regime with strong disorder), we remind 
the reader that the probability distribution of the difference, 5^ — Xa+i ~ Xa is similar to the 
Poisson distribution in the localized regime (Fig. I47> . Therefore, it is natural to assume that the 
repulsion between the two levels in the Jacobian ( I213> is weaker than in the diffusive regime. 
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Fig. 65. Left: The disorder dependence of LKn for the 3D anisotropic Anderson model, given by Eq. <215> 
with t — 0.4. The symbols correspond to the different disorder. The data confirm that LKn decreases 
(increases) with L for < Wc {W > Wc), respectively. Here, the critical disorder Wc ~ 9.5. Note, 
LK\\ does not depend on the system size when W = Wc- The right panel shows that 712 1 when 
W < Wc and L — > 00, but 712 decreases with L when disorder W > Wc (insulating regime). At the 
critical point, 712 ~ 0.25 does not depend on the system size [167]. 



Consequently, we assume that ^ab ^ in the case of strong disorder. A similar conclusion was 
derived also in Ref. [91]. 

The second assumption made in Ref. [165] was that Kaa ^ Kab for a 7^ 6. In particular, in 
the localized regime, one can assume that 

Kaa (X L°. (214) 

To understand the physical meaning of J214t . note that 




(215) 



is nothing but the inverse participation ratio for the transversal wave function on the opposite 
side of the sample. In analogy with the properties of Iq, discussed in Sect. 15.11 we conclude 
that the condition ( I214t reflects the non-homogeneity of the distribution of the electron on the 





L X Kn 


712 


W ^Wc 


const 


1 

const 



Tab. 3. The size dependence of parameters LKu and 712 in the metallic, localized and critical regime. 
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Fig. 66. Probability distribution of a strongly disordered three dimensional system. Symbols show data 
form numerical simulations, solid line is the solution of the generalized DMPK equation which gives the 
same value of (In g) and with 712 = T/2. Dotted line is the distribution obtained for the same value of F 
but with 712 — 1. 



opposite side of the sample, which follows directly from the existence of only the single possible 
path through the sample (Fig. I71l i. 

Conjecture (I214> was numerically tested in Refs. [166, 167] by the direct numerical calcula- 
tion of the matrix Kai- The numerical method of calculation of the conductance, described in 
Appendix lE.2l can be easily generalized for calculation of the eigenvectors u. To avoid the prob- 
lem with evanescent modes, the calculations were performed for the anisotropic 3D Anderson 
model, defined by Eq. (|9j with t = 0.4. This model exhibits the metal-insulator transition for 
critical disorder Wc ~ 9.5. 

The numerical analysis confirmed that the parameters K indeed depend on the disorder. In 
Fig. |65]we see that the size dependence of parameters LKu and 712 differs depending on 
whether the system is in the metallic, localized or critical regime. Note, since Ai < A2 < . . ., 
and the matrix tH has eigenvalues (1 + Xa)^^, the parameter Ku contains information about 
the spatial structure of the first eigenvector, ui, which corresponds to the largest eigenvalue, and 
712 is a measure of the overlap of eigenvectors related to the two largest eigenvalues of tH. 

In the metallic regime, LKu oc L~^, and 712 converges to 1, as was assumed in the deriva- 
tion of the "classical" DMPK equation. However, in the localized regime, LKu increases with 
the system size, which indicates that indeed Ku is constant, in agreement with Eq. ( I214> . Also, 
7i2 decreases when W > Wc being cx L^^ in the limit of large system size. 

At the critical point, both parameters, LKu and 712, converge to the size independent con- 
stants, indicating that they could be used as the order parameters for the calculation of the critical 
parameters of the model by scaling theory. The size dependence of LKu and 712 is summarized 
in TablelH 

With known size dependence of Ku and 712, we need to estimate also other parameters, 
Kab- Detailed numerical analysis, performed in Ref. [167] confirmed that these quantities can 
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be expressed as a simple functions of indices a and b. It is also reasonable to assume that 
Kaa — Kii and jab ~ 7i2- This reduces the number of free parameters to two. The first 
one, 

measures the strength of the disorder, and the second one, 712, influences the mutual correlation 
of channels in the generahzed DMPK equation. In the limit of 

712 ~ ^ « 1' (217) 

the probability distribution p{g) can be obtained analytically by solving the generalized DMPK 
equation, given by Eq. ( I212> . by methods developed in Refs. [109, 110] The obtained conduc- 
tance distribution is shown in Fig. |66land compared with the numerically calculated distribution 
p{lng). Parameter F was chosen such that (In g), equals to the numerically obtained value. We 
see that the analytical model reproduces qualitatively correctly the numerical data. The quantita- 
tive difference is probably due to oversimplification of the model, which neglects the differences 
between Kab for higher channels. 

The generalized DMPK equation represents the most promising step toward the analytical 
description of the localized regime. However, the main assumption which allows the analytical 
solution is that 712 is very small. This is not true at the critical point, as can be seen in Fig. l65l 
Therefore, it is not known at present, how accurately the generalized DMPK equation describes 
the critical regime. 



15.2 Random matrix model of the Anderson transition 

We have shown in Appendix [C] that in the diffusive regime, the probability distribution for the 
eigenvalues A of the transmission matrix can be relatively easily obtained from the assumption 
that the matrix (tH)^^ belongs to the orthogonal class of random matrices. The form of the 
distribution was determined by the additional constrain that the density <t{x) of parameters x, 

a{x) = {Y,S{x-Xa)), (218) 

a 

is constant for x < L^/i- This constraint follows directly from the linear dependence, (xi) cx a, 
shown in Fig. EH 

Since the spectrum of x exhibits the universal properties also at the critical point and in the 
localized regime, it is tempting to try to generalize the random matrix analysis also for transport 
beyond the diffusive regime. This expectations are inspired by the numerical data for (t{x), 
shown in Fig. ^\ 

In Sects. [TOlwe found that contrary to the diffusive regime, the density <t{x) is not constant 
in the critical regime. Since (a^a)^ cx a at the critical point (Fig. I41> . the density must be linear, 

crcrit(a;) cx X, (219) 

at least in the lower part of the spectra. Note, only this part of the spectra is relevant for the 
description of the transport. The linearity of the density is confirmed numerically in the left 
panel of Fig. |^ 
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Fig. 67. Left: The density (7{x) calculated numerically for the 3D Anderson model with disorder W = 
6, 10, W ^ Wc = 16.5, and W = 25 and 32 (from left to the right). Note, at the critical point, the 
density is linear for x < 10. Right: (j{x) for the strongly localized regime with disorder W — 25, 32, 45 
and 55. To show the universality of the density, the data are shifted on the horizontal axis by the difference 

(a;i(W))-(a;i(W = 55))[168]. 



The right panel of Fig. I^shows that in the localized regime, the density cr(a;) is a function 
of the difference, x — 2L/X. This follows directly from Eq. ( I161> . 
It is therefore natural to postulate that the density (t{x) behaves as 

a{x) = c{W) x[x + a{W)] , (220) 




Fig. 68. The density, a{x), calculated from the random matrix model, given by Eq. j220> . Various values 
of parameters a are given in the legend [168]. 
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where the function a{W, L) changes the sign at the critical point [91]. 

The disorder and size dependence of functions a and c are not known. The analysis of the 
size dependence of of parameters x in the critical region [169] confirmed that 

r +2L/^{W) W-^Wc 
a{W,L) - <^ W = Wc (221) 

[ -2L/^{W) W>Wc 

and that c{W,L) is a size independent constant for W = Wc. In Eq. (122 U . £,{W) is the 
correlation length introduced in Sect. 112.11 In strongly localized regime, £,{W) equals to the 
localization length A. 

We can verify that the density (t{x), given by Eqs. ( I220> and ( I22U reproduces numerical data 
for parameters x. Indeed, inserting in Eq. ( I220> we find that in the metallic regime, a{x) ^ const 
since x ^ 2L/£^{W) (we remind the reader that (xa) — L/ {N£)a in the metallic regime). Since 
a — cx at the critical point, we immediately recover that (xq)^ cx a. In the localized regime, 
a = c{x — 2L/^), so that the entire spectrum of parameters x shifts by 2L/S^{W) for higher 
values of the disorder W. 

Using the method explained in Appendix lC.il it was found in Ref. [91] that the distribution 
of parameters x at the critical point is given by the universal distribution Eq. J310l l and J312> . but 
with a cubic one particle potential, 

V„it = cx\ (222) 

Then, applying the methods of orthogonal polynomials, [44], the density <j{x) can be calculated 
within the random matrix model. Results, shown in Fig. |68l agree qualitatively with results of 
the numerical simulations, shown in the left panel of Fig. |^ 

An interesting consequence of the expression \22Q\ is the following scaling relation for higher 
Lyapunov exponents z^: 

h - a{W, L)f - [zj - a{W, L)f = \{^-]), hJ<L (223) 

which enables us to find the critical parameters of the 3D Anderson model from numerical data 
for only one system size [169]. It also provides us with numerical evidence that the model J22U . 
based on the change of the density at the critical point, is correct. 

The simple form of the potential J222t indicates that the statistical description of the critical 
regime might be as simple as the analysis of the diffusive regime. However, we are not aware of 
any microscopic model which confirms the above phenomenological model. 

16 Conclusion 

We have discussed in this paper the main transport properties of the single electron disordered 
electronic systems. The quantum character of the electron propagation is responsible for the 
wide variety of interesting transport phenomena, which we demonstrated numerically. 

In the limit of weak disorder, the weak localization and anti-localization corrections to the 
conductance and universal conductance fluctuations were observed and the conductance was 
found to exhibit the universal conductance fluctuations. In the localized regime, the exponential 
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decrease of the wave functions and the conductance with the system length has been demon- 
strated and discussed. 

One of the most important phenomena in the locaUzation theory is the absence of the self 
averaging of the conductance. In the metallic regime, the electron wave function is spread over 
the entire sample, and is very sensitive to any change of the disorder configuration, as well as 
to the change of the boundary conditions. This remarkable property survives also in the limit 
of an infinite system size. The sensitivity of the energy spectra to the change of the boundary 
conditions can be used as a measure of the electron locaUzation, and defines the main parameter 
of the localization theory, the conductance g. 

The statistical properties of the conductance in the strongly localized regime are even more 
complicated. The transport is possible only by tunneling between localized centers, and the 
conductance decreases exponentially with the size of the sample. Still, a few samples exist in the 
statistical ensemble which possess rather high conductance. These samples determine the mean 
conductance, which might be in orders of magnitude larger than the typical (the most probable) 
conductance of the ensemble. For a given sample, the conductance is also very sensitive to the 
distribution of random impurities. Equivalently, a small change of the Fermi energy might cause 
a huge change in the conductance. 

The critical regime of the metal-insulator transition has been discussed in detail. First, we 
have analyzed the statistical properties of the conductance, and we have presented numerical data 
for the critical conductance distribution. Then, the scaling theory of localization is introduced and 
verified by numerical scaling analysis of Lyapunov exponents and conductance distribution. In 
order to compare the numerical data for the critical parameters with predictions of the analytical 
theories, the critical electronic transport on lattices with dimension close to the lower critical 
dimension, dc — 2 as well as with dimension d = 4 and d = 5 was simulated. 

We have tried to convince the reader that the numerical data, collected within the last 25 years 
supports the validity of the one parameter scaling theory of locaUzation. We beUeve that the 
numerical data provides us with the most reliable estimation of the critical exponents in various 
disordered models. We hope that the disagreement between the results of numerical simulations 
and theory wiU motivate further theoretical research [155]. 

Finally, we want to mention some other aspects of the single electron localization. We have 
not discussed the transport in systems with correlated disorder. The mobility edge in the ID 
system with correlated disorder has been found in Ref. [171] (see Ref. [172] for review). Special 
attention is deserved also for disordered systems with chiral symmetries [173]. Also, application 
of the concepts of the electron localization to the propagation of the electromagnetic [22, 24, 
93, 174], and acoustic [175] waves in disordered media opens a new field for studies of the 
locaUzation and promise further development of the localization theory. 
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Fig. 69. The composition law for the transfer matrix. 

A Properties of the transfer matrix 

We summarize here the most important properties of the transfer matrix. For details, we refer the 
reader to Refs. [77,177]. 

First, we derive the expression for the transfer matrix, given by Eq. ( I54> . From Eq. ( I5U we 
obtain 

C = t+A + r-D 

B = r+A + t-D ^^^^^ 

If we have N incoming channels, then the amplitudes A, B, C and D are the (complex) vectors 
of length , and the transmission and reflection amplitudes are the matrices of size N x N. 
We express D from the second equation (I224t as 

D = {t-)-^B - {t-y^r+A (225) 

and insert it in the first equation (|A|i. We obtain that the transfer matrix is given by Eq. ( I54l i. 

+ -r-(t-)-V+ r-(t-)-i 



(226) 



Note, and are the transmission and reflection amplitudes of the wave coming from the 
right side of the sample, and <+ and r+ are the transmission and reflection amplitudes of the 
wave coming from the left side of the sample. 

A.l The composition law 



The transfer matrix given by Eq. (I226> fulfills the composition law. If the sample consists of two 
subsystems, then the transfer matrix T12 of the whole sample can be calculated from the transfer 
matrices of the two subsystems as 



L12 = 



T2T1. (227) 



The resulting transfer matrix T12 has again the form ( 154k This composition law enables us to 
calculate the transmission of a complicated structure from the transfer matrices of its parts. In 
particular, the transmission through two scattering centers is given by the relation 

ti2=tl [i-rtr^r't^- (228) 
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Similarly for the reflection we obtain 

r+ =r++t+[l- r+r^] r+q . (229) 

A.2 Symmetries 

The elements of the transfer matrix are not independent. Various physical symmetries provide 
important relations between the transmission and reflection amplitudes. Here, we investigate the 
consequences of the flux conservation and of the symmetry with respect to the inversion of time. 

The flux conservation requires that the flux on the right-hand side of the sample equals to the 
flux on the left-hand side. This gives the following constrain to the transfer matrix [77], 



Inserting 

T=(l'' I'' ] (231) 



T21 T22 
we obtain from Eq. (I230> that 

TLT22 - r/2Ti2 = 1 (232) 

and 

TliTii - rJiTsi = 1. (233) 
We also have 

TIiTi2 - T^iT22 = 

Tl2Tn-T^2T2i = 0. 



With the use of the explicit form of the transfer matrix, given by Eq. (I226> . we obtain from Eq. 
(I232> the relation 



{t-)h- + (r-)V- = 1. (235) 
The second equation ( I234> gives 

r+ = -{t-)-^^r-)h+, (236) 
which can be written in the form 



r 



-{t-)-' = -{t+)-'Hr+y. (237) 



Inserting ( I236> into Eq. ( I233> gives, after some algebra, the relation 

it+)h+ + (r+)tr+ = 1. (238) 



The equations ( I235> and ( I238t are a direct consequence of the flux conservation. They tell us 
that the electron either transmits through the sample or is reflected back. 
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The relation ( I237> enables us also to simplify the expression for Tu, 

Til = t+-r-(t-)-V+=t+ + (t+)-it(r+)tr+ 

= (t+)-it [{t+)H+ + (r+)tr+] (239) 

The time reversal symmetry implies that the transfer matrix fulfills the relation 

^ T(*? i^=T*. (240) 



1 / V 1 



From Eq. ( I240t we obtain that for systems with time reversal symmetry the transfer matrix 
can be written in the form 



2^11 Ti2 
T*2 T*i 



(241) 



where Tu and T12 are the NxN complex matrices. Also, from Eq. ( I230> we see that |det T| — 1 
and from Eq. ( I240t we get (det T)* = det T. Consequently, det T = 1 when the time reversal 
symmetry is preserved. 

From Eq. ( 124 1> we also find that in the case of time reversal symmetry the matrices t+ and 
t~ are related by the relation 

t+ = {rf. (242) 

A.3 Parametrization of the transfer matrix 

Consider now the system with time reversal symmetry. Then, the transfer matrix, T, is deter- 
mined by two complex matrices, Tu and T12, which are completely determined by 4A^^ real 
parameters. However, these parameters are not independent of each other, since the matrices Tu 
and T12 must fulfill the relations of flux conservation. 

We insert the expression JMTl into Eq. ( l230l and from Eq. ( l233l we obtain the relation 

tI.Tu - Tl^T*^ = 1, (243) 

which implies N'^ additional relations between the elements of the matrices Tu and T12. Other 
constraints are given by the first of Eq. J234t . 

TI1T12 - T^2T*u = 0, (244) 

which can be written in the form 

= A, where A^T^^Tu- (245) 



Equation ( 1245 > requires that the complex matrix A is symmetric. This requirement implies 
No{No — 1) relations between the elements of the matrices Tu and T12. Finally, we have that 
the matrices Tu and T12 are fully determined by N{2N + 1) real parameters. 



Following Ref. [178] let us look for the solution of Eqs. ( I243> in the form 

Tu = uLv and T12 = u'L'v', (246) 
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where u, u' , v, and v' are the unitary matrices and L, L' are the diagonal real matrices. Inserting 
into Eq. ( I243> we obtain that 

v^L^v- iv'f{L')^v')* = 1, 

(247) 

v^Lu'fu'L'v' - {v')^L'{u')^u*Lv* = 0, 

which can be solved with 

u' = u, v'^v*, {L'f = L^~l. (248) 

Since the complex unitary matrices u and v are determined together by 2 A^^ real parameters and 
the real diagonal matrix L needs N real parameters, we see that the solution ( I248t is consistent 
with our estimation of the number of free parameters of the transfer matrix. 

Taking L = vT+A we obtain that the transfer matrix can be parametrized in the form 

u o \ f VTTx Vx \ f V 



where u, v are the A'o x A'o unitary matrices. We also use the following parametrization of the 
diagonal elements AqI 

Aa = i [cOshXa- 1]. (250) 

It was shown in Ref. [178] that parametrization J249t is the most general one. Any transfer 
matrix of the orthogonal system can be expressed in the form given by Eq. ( I249I I. A similar 
parametrization can be derived for the transfer matrices with unitary and symplectic symmetry. 
For unitary symmetry, the transfer matrix has the form 



ui o \ f VTTa Va \ / vi 
U2) [ Va x/TTa j I V2 



(251) 



i. e., it is determined by four unitary matrices ui, U2, f 1 and V2 and by the diagonal matrix A. A 
detailed derivation of relations ( 1249125 if can be found in Ref. [178]. 

A.4 Transfer matrix vs conductance 

From the flux conservation ( I230> we obtain the following relations for the transfer matrix: 



and 



With the use of these equations we obtain 



[TtT]-.( J _;)TtT( J _M. (254) 
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Using Eqs. ( I254t and ( I234> we obtain 



[TtTl + TtT = f "2 t M ■ (255) 

Inserting into Eq. (I255t the matrix elements Tn = (/;+)^^^, given by Eq. J239> . and T22 = 
{t~)^^ we obtain the formula of Pichard [179] 

° ^ =i [TtT+(TtT)-i+2]"'. (256) 



{t-)H- I 4 



Note that relation ( I256t follows dkectly from the requkement of flux conservation. 

We can also use the parametric form of the transfer matrix, given by Eq. ( 125 1> . and we 
express the inverse of the r.h.s. of Eq. (I256t in the form 



Comparing the l.h.s. of Eq. ( I256> with Eq. ( I257K we obtain that 

t+{t+y =vl{l + X)~^vi (258) 

and 

=4(1 + A)~^W2. (259) 

Both equations are equivalent if time reversal symmetry is preserved. Indeed, (t^)'' — {{t^)^)* — 
{t~)* and V2 = vl. 

The matrices vi and V2, obey the relation of unitarity, v^^ = v'' . Using this relation we can 
express the Hermitian matrix {t^)H^ in the form 

{t^)h- ^V2'y^V2- (260) 

Thus, the eigenvalues of the matrix {t^)H^ are (1 + Xa)^^ and the unitary matrix V2^ contains 
in its columns the corresponding eigenvectors. 

Using the eigenvalues of the matrices {t^yt^ and (t+)tt+ we can express the transmission. 



(261) 



A.5 Open channels and evanescent waves 

In previous Sections, we assumed that there are only propagating states in the leads. This might 
not be true in real systems. For instance, for the 3D model, we see from the dispersion relation, 
given by Eq. ( I20L that the propagation along the leads, (which is the z direction in our notation), 
is determined by the z-component of the wave vector, k^, given by the relation 

2 cos kz = E — 2 cos k^ — 2 cos ky (262) 
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where the transversal wave vector is given by 

kn^ = -j—nx, Hx = 0,1, ■ ■ ■ , Lx - I (263) 
for the periodic boundary conditions and 

TT 

kn^ = 7 TT"^' "i: = 1, ■ • ■ , -^a: (264) 

Lx + 1 

for the hard wall boundaries. Similar relations hold for ky. 

Using the eigenfunctions (x) and $„ (y), related to the eigenvalues of kx and ky, we 
find that the wave function in the leads is 

*(f) = J2 (y) (265) 

^ V2« sm kz 



The wave vector k^ — kz{nx,ny) possesses N = LxLy values given by Eq. (I262> (not neces- 
sarily different from each other). The propagating solutions exist only for such values of k^, for 
which 

\2 cos kz{nx,ny)\ < 2, (266) 

i. e. when k^ is real. Otherwise, k^ is purely imaginary, and determines the exponential decrease 
of the wave function with distance form the left (right) boundary of the sample. Therefore, we 
have to distinguish between A^o and , where A^o is the number of propagating channels with k^ 
real, and N is the total number of channels. 

Since we assume that both leads are semi-infinite, it is clear that the evanescent waves, emit- 
ted from the reservoirs, are not able to reach the sample. Therefore, the size of the transmission 
and reflection matrices is x N^. The necessity to distinguish between N and complicates 
the calculation of the conductance. Indeed, the size of the transmission matrix is No but the 
size of the matrix T22 defined by Eq. (123 1> is Nx N. If we order the eigenvectors in the matrix R 
in such way that the eigenvectors with index 1 < a < A'o correspond to the propagating modes, 
and the remaining eigenvectors correspond to the evanescent modes, then the transmission is 
given only by the No x A'o sub-matrix [T22]ab, with a,b < No'- 



I L^22 
ab=l 



(267) 



Other matrix elements of T22 correspond to the scattering of the electron into the evanescent 
channels. Note, the composition laws for the transmission and reflection of the system given by 
Eqs. ( I228l229t are not valid in this case because evanescent waves between two subsystems can 
also contribute to the transport. 

It is important to underline that the analytical theories discussed in next two Appendices 
assume implicitly that there are no evanescent wave in the leads. 

In numerical simulations, we can avoid the evanescent waves for the band center (E = 0) by 
using anisotropic models, given, for instance, by Eq. (|9j. However, the evanescent waves cannot 
be completely excluded if we are interested in the transmission of the electron with energy close 
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Fig. 70. Small segment of length SL^ is added to the sample of length L^. 

to the band edge. Indeed, from the E dependence of k^, given by Eq. ( I262> . we see that iVo is 
maximal for the band center and decreases to zero when the energy E approaches the band edge 
of the unperturbed system. Therefore, we have to consider evanescent waves in the numerical 
calculation of the conductance. The calculation of the conductance is described in AppendixlEI 
Note, there are no propagating channels for the energy E outside of the unperturbed band. 
Therefore, we are not able to calculate the conductance for such energies. This constrain com- 
plicates the scaling analysis of the conductance along the critical line shown in Fig.|6| 

B Dorokhov-Mello-Pereyra-Kumar (DMPK) Equation. 

The transfer matrix of the disordered system is a statistical variable. One can ask if there is a 
possibility to calculate the probability distribution p(T), which is a joint probability distribution 
of all matrix elements of the transfer matrix, T. This problem was solved by Mello et al. [7]. For 
the transfer matrix, given by expression ( I249l l. they derived the equation for the joint probability 
distribution of the parameters A, which parametrize the transfer matrix in Eq. ( I249t . Since the 
derivation of the DMPK equation is rather difficult, we present here only main ideas and refer 
the reader to the original work, [7]. 

Consider a weakly disordered system of length L^, connected on both sides to the semi- 
infinite leads as shown in Fig. ^] There are N open channels in the leads, with N big enough 
to neglect all corrections of order N~^. The channels are equivalent to each other. In the leads, 
each channel is characterized by the same wave vector kz. The disorder is introduced inside the 
sample by random hopping terms between the channels. 

Suppose we have a sample of length and we know the probability distributionp^^ {T')dT'. 
Now, we add to the sample of length an additional random segment of length SL (Fig. YlQi . 
The resulting sample will have length + SL^ and is characterized by the transfer matrix 

T = T^-T'. (268) 

We suppose that the length SL^ of the additional segment is sufficiently small so that the reflec- 
tion coefficient of the segment is linear in SL^, 



i?+ =Tr(r+)tr+ oca^. 



(269) 
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In Eq. (I269> . the length £ is the mean free path and we put a = 1. We also neglect all terms 
proportional to higher orders of 6Lz- Equivalently, we can require that the trace, 



2N 



1 



6L, 

£ 



(270) 



is linear in SLz'' ■ 

We also need to know the probability distribution pslO^s) of matrix elements of the transfer 
matrix Ts- This distribution represents the main building block of the theory. In Ref. [7], the 
Ansatz 



PSL, (Taj oc exp Tr A5 



(271) 



was proposed. It represents the "most random" distribution of matrix elements of the transfer 
matrix Tg. We remind the reader that parameters A^- for the transfer matrix Tg are defined by Eq. 
( l249b . Clearly, Tr (A^) = NbL^jl. 

The probability distribution pL^+a^^ (T) fulfills the Smoluchovsky equation. 



(272) 



Inserting Eq. ( I27U into Eq. ( I272> leads, after rather complicated mathematical operations, 
to the equation for the joint probability distribution of eigenvalues A: 



5PL.({A}) 



1^ 



d 



d{Lz/£) f3N + 2-(3 J ^ dXa 
Here, J is the Jacobian, given by 



Xa{l + Xa)J 



dpLiiX}) 



dXa 



N 



J^\[\Xa 



X,\P 



(273) 



(274) 



a<b 



and (3 is the symmetry parameter, (3 — 1, 2 or 4 for orthogonal, unitary, or symplectic systems, 
respectively. Equation ( I273> is known in the literature as the DMPK equation. Note, the mean 
free path, £, is the only parameter which enters the theory. 

The DMPK equation, given by Eq, ( I273> does not contain any information about the matrices 
u and V. These variables were "integrated out" in the process of deriving Eq. J273> . The main 
assumption made in this process was that u, v and A are statistically independent. Next, it was 
assumed that the elements of the matrix v fulfill the relations 



{Kh'Vcd) = -^SabScd 



and 



[VcaVcb^daVdh) 



l+Sab 

N{N+1) 



Ocd, 



(275) 



(276) 



and similar relations hold also for the matrix u. These assumptions limit the validity of the 
DMPK equation to the systems with weak disorder. It is implicitly assumed that there are many 
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Fig. 71. The propagation of tlie electron througli the disordered sample. Left: If the disorder is weak, then 
there are many paths through the structure. We expect that the wave function of the electron on the opposite 
side of the sample is spatially homogeneous. This corresponds to the assumption <275l276t . Right: If the 
disorder is strong, there is in the best case only one path through the structure. Then, the wave function 
possesses sharp maxima at the opposite side of the sample; the positions of these maxima depend on the 
realization of the disorder in a given sample. It is clear that assumption <276t is not longer valid. Therefore, 
we do not expect that the DMPK equation, given by Eq. <273> . describes the transport in strongly disordered 
systems. 



paths for the electron to propagate through the sample. Then, the density of the electron on the 
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Fig. 72. Numerical verification of the validity of the DMPK equation. Left: The N = dependence 
of the {N + l)_ft'ii/2 for the quasi-ld systems L x L x AL and for the cubes . In the last case, also 
(A'^ + 1)K\2 is shown. The DMPK is applicable if these quantities converge to 1 in the limit of L — > oo. 
Right: The a dependence of Kaa for the quasi-ld systems and for the cubes. Data were obtained for the 
anisotropic Anderson model with t = OA with disorder W = 2 and W = A, which corresponds to the 
mean free path 1 = 9 and i = 1.8, respectively [167]. 
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opposite side of the sample is homogeneous, as shown in the left panel of Fig. ^\ 

Another constraint of the DMPK equation is given by the assumption of the equivalence of 
all incoming and outgoing channels. In the real world, the channels are never equivalent, since 
they are determined by the incident angle (or by the transversal momentum) inside the leads, as 
shown in Fig. ^] The assumption that all channels are equivalent removes from the theory any 
information about the topology of the leads. Further, in the derivation of the DMPK it is assumed 
that channels are equivalent also inside the sample, since disorder was introduced by the random 
hopping between any two channels. Thus, the DMPK equation represents the simplest "mean 
field theory". In spite of the above constrains, the DMPK equation is surprisingly successful in 
the description of the electronic transport in weakly disordered quasi-ld systems. The DMPK 
equation can be easily generalized to description of more realistic systems. For instance, instead 
of the Ansatz (127 1> we can use more complicated probability distribution, which reflects the 
properties of the system, inclusive anisotropy and topology [180]. However, such generalization 
does not bring any novel physical information. It only complicates the resulting formula for the 
transmission. 



B.l Numerical verification of validity of the DMPK equation 



In Fig. |72|we verify numerically the validity of the relation ( I276t . We consider the anisotropic 
Anderson model, given by Eq. (|9jl with t — 0.4 and calculate the matrix 



Kab ^ { y \UacnUbcr ) , (277) 




where, as usual, (. . .) means averaging over the statistical ensemble. Requirement ( I276> is equiv- 
alent to 

Kat = (278) 
with N = L'^. 

Our data show, that although the channels in the Anderson model are not equivalent to each 
other, the requirement j278> is perfectly fulfilled for the weakly disordered quasi-ld systems, as 
long as the mean free path, I, is comparable to the width, L, of the system Agreement is worse 
for 3D systems. This is natural, since the 3D systems are expected to possess slightly differ- 
ent conductance statistics than the quasi-ld systems (see Table [2 which compares data for the 
universal conductance fluctuations, var g. in the quasi-ld, 2D and 3D systems). When disor- 
der increases, the mean free path decreases and the transversal structure of the sample becomes 
important. This could be included into the theory if a more detailed model for the distribution 
p(^s) is considered. 



B.2 Conductance 

The DMPK equation can be directly used to calculate the first two moments of the conductance 
of the weakly disordered quasi-ld system. The method is described in Ref. [83]. When we are 

'Note, in the case of zero reflection, the matrices (i^ and t+ {t^)^ are diagonal and Tr T^T = 2N . 
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interested in the mean value of the function F{\), we multiply both sides of the DMPK equation, 
given by Eq. (I273> by F{X) and integrate over A. On the r.h.s., we obtain mean values of some 
other functions, Fi{X); repeating this procedure for Fi, we obtain a system of coupled equations 
for (F) and (Fi), which can be solved in the limit of cx). In this way, the mean conductance 
(in units of e^/ft.) was calculated as 



-4 



(279) 



The first term, Ni/L^ can be identified with the conductivity, a. Other terms represent the quan- 
tum corrections. Of particular interest is the diffusive regime, which is defined by the following 
relations between and £: 



£<s:L,<s: N£. 



(280) 



If the conditions ( 12 8 Oi l are fulfilled, then we can neglect the last term in Eq. ( I279> . Assuming 
also that the system is sufficiently long, so that 



£2 

Nj^ « 1, 



(281) 



we can neglect also the 2nd term on the rh.s. and we finally obtain the following estimation of 
the mean conductance: 



mi 



(282) 



The correction, 5g — —1/3, is the universal weak localization correction to the conductance of 
the quasi- Id system. This is a fully universal value, independent of any details of the system 
provided that the conditions ( 1280128 l( l are fulfilled. The numerical verification of the relation 
( I282t is shown in Fig. |28l 

In a similar way, Mello and Stone calculated in Ref. [83] the variance. 



varg=(5^)-(.g)^ = — . 



(283) 



The variance of the conductance in the quasi- Id system is universal, given only by the physi- 
cal symmetry (by the value of the parameter jj). The same result was obtained in Ref. [9] by 
perturbation Green's function analysis. 

The universal relations (I282> and ( I283t hold only when the system is long enough (as required 
by Eq. ( I281» . In this case, the electron is scattered many times inside the sample so that its wave 
function and phase are sufficiently randomized. As required by Eq. ( I280> . the length of the 
sample should not be too long, otherwise the effects of localization become important. This is 
analyzed in the next Section. 
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B.3 Limit > N£ 

In the limit of Lz ^ N£, the localization appears. As shown in Appendix|Dl the wave function 
decreases exponentially in the limit of infinitely long system. This must be true also for the solu- 
tion of the DMPK equation. Therefore, parameters Xa = (cosh cCa — 1 ) /2 increase exponentially 
with the system length. If we order the parameters A, 



Ai < A2,< . . . ,< Aat, 
then the Jacobian, given by Eq. ( I274t . reduces to 



(284) 



(285) 



and the DMPK equation splits into the TV independent equations for the probability distributions 

p(Aa), 



d{L/£) f3N + 2-f3 " dXa 



a+ldpUXa) 



dXa 



(286) 



In the derivation of Eq. ( 12861 we used that Aq ^ 1 so that AQ(Aa + 1) ~ A^. Substituting 

Xa — expxa, and using dp{Xa)/dXa — X^^dp{xa)/dxa, we obtain 



dpLA^a) 



dXa 



dx'i 



d{Lz/i) (3N + 2-(3 
which is the diffusion equation. The solution of Eq. (|2^} is 

1 {Xa-{Xa)f 



PL Ax a) = — ==CXp- 



/27ra 

with the mean value, 

, ^ l + /3(a-l)L, 



2a 



(3N + 2-P e ' 



and variance, 

cr = 2- 



(287) 



(288) 



(289) 



(290) 



f3N + 2-(3 I 

Thus, in the limit of L^/t ^ 1, the parameters Xa become statistically independent. They pos- 
sess Gaussian probability distributions with the same variance, a, and with mean values which 
increase linearly as a function of index a: 

13a 



(291) 



l3N + 2-(3 I 

The numerical data for the probability distributions p^^ [xa] are shown in Fig. ^3] 

When calculating the conductance, we can neglect all contributions from the higher channels 
We obtain 



9 



(292) 
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Fig. 73. The statistical distributions of the parameters Xa for the weakly disordered quasi-ld system with 
the box disorder W = 6. The size of the system is 10 x 10 x 200. Since disorder is weak, (the critical 
disorder Wc = 16.5 for the box distribution of random energies), we have (xa) ~ a{xi), in agreement 
with the DMPK equation. The length of the system, = 200, is sufficient to localize all electronic states. 
AH parameters Xa are large and posses the Gaussian distribution, as predicted by Eqs. <288l290l . 



which confirms that the conductance decreases exponentially when the system length increases. 
We also obtain that the logarithm of the conductance. In g = —xi + In 4, possesses the Gaussian 
probability distribution with the mean 

(lng) = -2 — +ln4 (293) 

^ ^' (3N + 2- f] £ ^ ' 

and with the variance, 

var .g = - 2 (In .g) - In 4. (294) 

The DMPK equation predicts the universal statistical properties of the conductance also in the 
localized regime. 

However, it is important to mention that the localization described by the DMPK equation 
does not correspond to the true insulating regime. Note, the disorder is still weak in the system. 
The localization appears only due to the increase of the length of the system, L^- We discuss 
in Sect. II It that the localization in the weakly disordered quasi-ld systems differs quahtatively 
from the locahzation in the strongly disordered 3D systems. 



C Random matrix theory 



Random matrix theory was developed 50 years ago for the description of the statistical properties 
of large matrices, which appear, for instance, in studies of the energy spectrum of big nuclei. 
The experimental data showed that the energy spectra possess common universal features - for 
instance, the repulsion of the neighboring eigenvalues, similar to that shown in the left panel of 
Fig. It was suggested that the Hamiltonian of a big nucleus can be approximated by a random 
matrix. Then, the qualitative properties of the energy spectra can be obtained from the general 
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properties of the random matrices. Recently, the random matrix theory for the transfer matrix 
was developed to explain some universal properties of the electron transmission in the diffusive 
regime. 

We present here the formula for the joint probability distribution p(A) of the eigenvalues of 
the random matrix. More details about the theory can be found in Refs. [44, 179]. 

Consider a real symmetric N x N matrix H with random elements Hat- Let us define the 
joint probability distribution of its elements, P(H)(iH. The measure cffl is given as 

N 

dH^Y[dHuY[dH,j. (295) 

1=1 i<j 

Our aim is to find the expression for the joint probability distribution of the eigenvalues of the 
matrix H. Consider only the orthogonal symmetry class. Then the matrix H is real and symmet- 
ric. It can be diagonalized like 

H = QAQ"\ (296) 
with the helps of the orthogonal matrix Q, which fulfills Q^^ = Q^. It is evident that 

QQ'^ = 1. (297) 
Then, each element Hij can be expressed as 

N 

i7y ^^A^Qfe^Qfej. (298) 
fc=i 

Since H is a real symmetric matrix, it is fully determined by the N{N + l)/2 real independent 
parameters, Hij, i < j. We want to express the measure dH in terms of the eigenvalues Aj and 
parameters Xa of the matrix Q (there are exactly N{N — l)/2 independent parameters Xa which 
determine matrix Q). We need to find the Jacobian, J, of the transformation 

N N 

dH^Yl dHu n dH,, = JAr({A}, {Q,,}) J] dK J] ^a- (299) 

i—1 i<j i—1 a 

The Jacobian matrix can be schematically written in the form 

/ dHu dH^ 

Ik IhI h (300) 



\ dxfj, dx^ 



where the upper left quadrant is the N x N matrix, (i counts columns and a counts rows), and 
the lower right quadrant is a square matrix of size N{N — l)/2 (i < j counts columns and /i 
counts rows of the matrix). 

With the help of Eq. (I298> we see that the first N rows of the matrix J do not contain A. Also, 
all elements in the lower N{N — l)/2 rows are linear in the eigenvalues. Therefore, J = det J 
is a polynomial of order of N{N — l)/2 in the eigenvalues A. In order to find the form of this 
polynomial one has to realize that if the matrix H possesses two degenerate eigenvalues, Xa — Xb 
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for a ^ b, then diagonalization of H is not unique and the matrix J~^, which determines the 
transformation inverse to (I296> . must be singular. Therefore, J = det J = for the degenerate 
matrix. Then, we easily find that 

J =Y[\Aa-Ai,\Fi{x}), (301) 

a<b 

where the function F{x) does not contain any information about the eigenvalues A. We obtain 

P(a)(M = P{{A})J{{A}dAP{x)dx. (302) 

Now, we can integrate out all parameters x to obtain the probability distribution for eigenvalues, 
A. 

In the case of unitary and symplectic systems, the Jacobian changes to 

Jp=(xll\Aa-Abf, (303) 

a<b 

where /3 — 2 {f3 — 4) for the unitary (symplectic) symmetry, respectively. 

The expression ( I303l l plays an important role in the theory of random matrices. It tells us 
that there is a zero probability to find degenerate eigenvalues. The level repulsion is a typical 
property of random matrices. We have shown already in Fig. that the spectrum of weakly 
disordered Hamiltonian exhibits level repulsion. 

To find an explicit form of the probability distribution P{A), we look for the probability 
distribution P(A) in the form 

P(A) = J^exp^F(A,). (304) 

a 

Choosing the expression ( I304t . we implicitly assume that F{A) is a function of only one eigen- 
value. We will see later that this assumption might not be fulfilled in some applications of random 
matrix theory. 

In the limit of ^ oo, we introduce the level density, o'(A) and we write the distribution 
( I304t in the continuum form, 

P(A) = e-f^" (305) 

where 

H = -^J dAa{A)F{A) -^J J dAdA'aiA)a{A') In |A - A'|. (306) 

To specify the function F, we must introduce some constrains. We can, for instance, require a 
specific form of the density, a{A) and look for the "most random" distribution P(A) by mini- 
mizing the "Hamiltonian" H. The condition 

SH 

— = (307) 

0(7 

leads to the following expression for F: 



F{A) = ~P J dA'aiA') In |A - A'|. (308) 
In the next Section, we apply the above formulas to the transfer matrix. 
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C.l Application to transfer matrix 

Assuming that the transmission matrix 

[thy^ ^ v^l + X)v (309) 

of the diffusive system belongs to the class of random matrices, we immediately obtain that the 
distribution of parameters A has a form 

P({A}) = e-^" (310) 

with 



poo 

H^-Y,ln\Xa-Xb\ + Yl / dAa(A)ln|A-Aa|. 

a<b a Jo 



(311) 



The first term is the Jacobian, and the second term is the one particle potential, given by Eq. 

ll308li . 

It is more useful to express the distribution ( I310l31l| l in terms of variables Xa instead of A. 
Inserting Aa = (coshxa — l)/2 into Eq. ( I31U we obtain 

N ^ N 

H{{x}) = — In I coshxa — coshxhl — — In | sinhxa| + V{xa). (312) 

a<b 'a a 



The second term on the r.h.s. of Eq. (I312t is the Jacobian of the transformation \ ^ x, and we 
have used a{X)dX — a{x)dx. The last term is the one particle potential. 



V{x) 



/>OC 

/ dx' CT(a;') In I cosha; — coshx'|. (313) 

The density <7{x) can be estimated from the numerical data. Fig. ^^shows that {xa) ^ a. 
This leads to the constant density 

'^^''^^ I otherwise. (314) 



We insert expression ( I314> into Eq. ( I313> . and we calculate the potential ^^(a;) following the 
method given in Ref. [181]. The first derivative of Eq. ( I313> reads 



dVix) Nl f°° sinha; , , 

^ - ' dx'. (315) 



dy, (316) 



dx Lz Jo cosh 
With the new variable, y = , we obtain 

dvjx) _m r 1 _ 

dx Lz Ji [y - e-^ y - 



with A ~ cxpLz/£ ^ 1. Since a; > 0, the first integral gives In \A — e~^\ — In |1 — e~^|. 
The function 1/(2/ — e^) possesses a singularity at e^. Since the (negative) contribution to the 
integral from 1 to cancels with the (positive) contribution from to 2e^ — 1, we obtain that 
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Fig. 74. The probability distribution p((5i2) of the normalized difference 12 — xi for three models in the 
diffusive regime. Two systems with orthogonal symmetry have the same distributional, given by Eq. j32L 
This confirms that the important property in the diffusive regime is not determined by the dimension of the 
sample, but by the randomization of the electronic wave function due to the multiple scattering. The third 
statistical ensemble, given by the 2D Ando model, possesses the symplectic symmetry (/3 = 4) The solid 
lines are the Wigner distributions. 



the second integral gives In | ^ — | — In | — 1 1 . For A >> the difference of the two integrals 
is In |e^ — 1| — In |1 — e^^| = a;. We finally obtain 

(317) 

ox Lz 

Thus, we have the quadratic one particle potential, 

N£ 

V{x) = —x\ (318) 

The probability distribution P{x) is similar to the statistical sum of the one dimensional 
Coulomb gas. In the Coulomb gas analogy, the parameters Xa represent the position of the ath 
particle confined in the quadratic potential V{x) and interacting with the other particles by the 
two particle logarithmic interaction. The parameter (3 is the "temperature". A detailed analysis 
of the probability distribution p{x) and of the consequences for the diffusive transport are given 
in a series of papers of Pichard et al. [179, 181-184]. 

The most simple consequence of the probability distribution p{x) is the "level repulsion". 
From random matrix theory it follows that the normalized difference 

5a.a+l = (319) 
{Xa+l - Xa) 

is distributed with the Wigner distribution pp. This was confirmed by numerical simulations 
[102, 181, 182]. In Fig. |74|we plot p{5i2) for three statistical ensembles for the 2D and 3D 
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orthogonal systems, and for the 2D Ando model. The data confirm that the distributions 73(^12) 
are indeed very similar to the Wigner distributions. 

Note, the term Insinha; in the "Hamiltonian" ( I312> can be interpreted as the interaction of 
the particle located at x with its "mirror image", i. e. the particle located at —x. Therefore, 
the distribution of p{xi) is given by the Wigner distribution, too. However, since the "interac- 
tion", Insinhcc, does not depend on (3, the distribution p(a;i) should not depend on the physical 
symmetry [181]. This is confirmed in Fig. 1751 

FigureESlconfirms that the spectrum of parameters Xa is linear in the diffusive regime, both 
in the 2D and 3D systems. 



(xa) = [l + (a-l)/3] (xi). 



(320) 



For /3 = 1, a more accurate estimation, which agrees also with results of numerical simulations, 
can be derived from random matrix theory, 

(321) 

where Xa are zeros of Laguerre polynomials Ln [90,91], Xa ~ [4(A^ + 1/2)] i^(a), where 
jo(a) is the ath zero of the Bessel function Jo{x). In the limit of a ^ 1, jo{a) w 7r(a — 1/4), so 
that Eq. (132 1> gives Xa oc a. 

The right Fig. ^6lshows the a dependence of parameters ga. 



ga = cosh ^ Xa/2. 



(322) 



Wigner surmise 
2D W=3 L=300 
• 3D W=6 L=14 
O 2D Ando model W=3 L= 1 00 




Fig. 75. The distribution p{xi) of the normalized parameter xi for three different models in the diffusive 
regime. The solid line is the Wigner distribution, given by Eq. <32t . Note, the distribution p{xi) does not 
depend on the physical symmetry of the model. 
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Fig. 76. The index dependence of the mean (xa) (left) and Qa (right) for the 2D and 3D weakly disordered 
systems. The data confirm the linearity given by Eq. <320L To calculate {xa), we have to diagonalize the 
matrix tH for each sample, extract Xa from the eigenvalue (1 + \a)~^. The mean value, (xa) is obtained 
by averaging over the statistical ensemble. Note, the data for large a are not present since cosh^(xa/2) 
exceeds the numerical accuracy of the computers for such large values of Xa- The conductance (g) = 16.1, 
4.3 for the 2D systems with W = 1 and W — 2, respectively, and {g) — 42.1 for the 3D system. In the 
last model, the anisotropy, t = 0.4 was used to avoid the evanescent channels in leads. 



Clearly, ga is the contribution of the ath channel and 

9 = Y.9a- (323) 

a 

Because of the rigidity of the spectrum, the fluctuations of Xa are of the order of the mean 
spacing, which is ^ Lz/Nl. We see that ga are very close to 1 for small a and exponentially 
small when a ^ N. There is Ncs channels for which < 1. Higher channels, with a > A^off 
give only negligible small contribution to transport Then, the conductance g « iVcff. The slope 
of the linear dependence, (xi), determines the conductance: if (xN^ff) — 1, then we have, from 
Eq. ( l320l l that 



(.,)=iVo. = l + ^(^-l). (324) 
This expression simplifies for (3—1: 

(g) = iVeff = (325) 
Since g — Ni/L^, we immediately see that 

(xi) = ^. (326) 

Thus, all parameters Xa increase linearly with the length of the system. The mean value, (xa), as 
well as the spacing between two neighboring parameters, (xa), decreases as N^^. 

The rigidity of the spectrum of parameters Xa was used for the explanation of the universality 
of conductance fluctuations in Ref. [92]. 
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Random matrix theory represents a powerful tool for the analysis of transport in disordered 
systems. More details are given in Refs. [35, 179, 185]. 



C.2 DMPK equation versus random matrix theory 

Both the DMPK equation and the random matrix theory were successfully applied to the transport 
in disordered systems. Originally it was believed that the theories are equivalent. This would 
mean that the probability distribution p{X), given by random matrix theory, solves the DMPK 
equation. 

However, these two approaches are not equivalent. The discrepancy between the two ap- 
proaches was observed when Beenakker [85] calculated the variance of the conductance from 
random matrix theory. The obtained result, var g — 1 /8(3 differs from the exact value, var g = 
2/15/3, obtained by the diagrammatic expansion [9] and from the DMPK equation [83]. This 
proves that random matrix theory is not exact. Beenakker and Rejaei [186] solved the DMPK 
equation for P — 2. Both in the metallic and localized regimes, the solution reads 



p{x) — exp — 



V{Xa) + ^ u{Xa,Xb) 



(327) 



1 a<b 

with the interacting term 

u{xa,Xb) = - i In I cosha;a - coshxfc| - ^ In \xl - xl\ (328) 
and with the one-particle potential, 

y{xa) = ttJ-xI - - \n{xa sinhxa). (329) 
ZL/z 4 

This solution is very similar to the probability distribution derived from random matrix theory. 
However, note, that it cannot be obtained from the Ansatz (I304> . since the two particle interac- 
tion, given by Eq. ( I328> . is not determined by the Jacobian only. Therefore, the one parameter 
function F{X), introduced in Eq. ( I304L is not sufficient for obtaining the distribution ( I327> . 

D Lyapunov exponent 

D.l One-dimensional case 

The Lyapunov exponent, 7(i?), of the one dimensional disordered system is defined by the rela- 
tion 

7(£;)- lim -^hi(vl/L+*L+i), (330) 
Note also, that the wave function 

l^-iJcxe-T^' (331) 
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decreases exponentially when increases. From Eq. ( I330> we see that the real part of the 
Lyapunov exponent determines the localization length, 

Re-f{E)^X-\ (332) 

Thouless [187] showed that the imaginary part is related to the density of states by the relation 

p{E)^-^lm-f{E). (333) 
TT oh 

In the ID systems, the Lyapunov exponent can be found analytically in the limit of weak 
disorder. For the Anderson model, defined by Eq. ( I75t . we find 

^(^) = ^^+ 2(4l'L^) ' ^^^^^ 

with E — 2cosfc. Note that the expression fails when close to the band edge, \E\ 2. 
A more detailed analysis [188, 189] showed that in this case 7 cx W^^^. In general, the weak 
disorder expansion exhibits a peculiar behavior in the neighborhood of energies E = 2 cos Trp/q, 
with p and q being integers. At the band center, E = 2 cos 7r/2 = 0, the fourth order term of the 
expansion diverges and gives rise to the correction of the 2nd order term [188, 189]. For instance, 
Re7 — for E = and the box disorder, defined by Eq. ( fTol i. since (e^) = 1/12. 

Correct expression, which takes into account the higher order terms of the expansion, gives 
Re7 = 14^^/105.4, which agrees with numerical data. 

The Lyapunov exponent plays an important role in the theory of localization. In the next 
Section, we generalize the ID case to the quasi-ld systems. Note that the Schrodinger equation, 

*n+i + (£« - + = 0, (335) 

can be written in the matrix form 



Then we can write 



(336) 



,(337) 



where we have introduced the transfer matrix, 

M(^^-) = nM„ = n(f-"" -J). (338) 

n=l n=l ^ ^ 

Thus, in the limit of 00, 7 is given by the eigenvalues of the matrix M*-'^^-'. Since M'^^-* 

contains the random energies e,, 7 is a statistical variable, too. 

Oseledec [190] proved that the probability distribution, ^(7) is Gaussian with the mean value 
(7) cx Lz and variance, var 7 cx (7) . Therefore, the variable 
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n-l n n+1 




Fig. 77. We cut the quasi- id system into the L x L vertical slices. The propagation of the electron in slice 
n is given by the Hamiltonian Tin- 



is a self-averaged quantity in the limit of long samples. 

Self-averaging of the Lyapunov exponent can be intuitively understood from the expres- 
sion J337> . The Lyapunov exponent represents the logarithm of the eigenvalue of the matrix 
(M^^^-*)^M'^^^-', which is the product of the random matrices. Naively speaking, 7 can be con- 
sidered as a sum of logarithm of eigenvalues of the random matrices M„, n ~ 1,2, . . . Lz- Then, 
applying the central limit theorem we expect, that both the mean value and variance of 7 are 
proportional to L^- 



D.2 Quasi-ld case 

Consider the quasi-ld system of the size L"^^^ x L^, and suppose that Lz S> L. For the purpose 
of numerical analysis, we divide the system into a set of vertical slices, shown in Fig. El^nd we 
write the Schrodinger equation in the form 

*„+i = - H„)*„ - ^'„_i (340) 

where Tin is the Hamiltonian related to the nth slice in Fig. El '^n is the vector which 
contains in its elements the wave function in sites of the nth slice. The length of the vector \I>„ 
isN = L'^-i. 

Equation ( I340> can be rewritten in the matrix form, 

"^"+1 ^ - M„ f ^ (341) 



^ n / \ ^ n — 1 

where the transfer matrix, M„, is given by the relation 

M„=(f-^" -J). (342) 

This relation is formally equivalent to Eq. ( I336> but now the matrix M„ has the size 2N x 27V. 
Note, det i\f„ = 1. Also, the eigenvalues of M„ appear in pairs, A and A^^. 
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Now we calculate the 2N x 2N matrix 

M(^^) = nM„ = n( f-^" -J). (343) 

n=l n=l ^ 

The matrix M'^^-* determines the behavior of the wave function at long distances. Since the 
system is effectively one dimensional, the wave function must decrease exponentially when 
the length of the system, L^, increases. This exponential decrease is given by the eigenvalues, 
exp —-fa, with 7a > 0. Clearly, the smallest 7^ determines the localization length. 

Since the matrices M„ contains the random variables, e„, aa are also statistical variables, 
and we need to know their probability distributions. Following Oseledec [190], we have that in 
the limit of 00 all the eigenvalues of the matrix 

[M^L,)^{L,Y/L, (344) 
converge to A a = e^" and Aat+q = e^^" where 

Ca = ^ X 7a. (345) 
The mean value, {(a), does not depend on L^, and the variance, 

varCa-(Ca>-(Ca)'(xi^. (346) 

^z 

Thus, Oseledec theorem states that the parameters (a are the self-averaged quantities in the limit 
of Lz — > oo. For a given realization for the disorder, the eigenvalues of the matrix ( I344> converge 
to their mean values. In numerical simulations, the length Lz is finite, so that the obtained 
numerical data, <^a, differs from the limited mean values, {(a)- Oseledec theorem enables us to 
estimate the typical difference, 

Ca - (Ca) - ^[TJ ■ ^^^^^ 

Thus, we can avoid the problem of statistical fluctuations in numerical calculations. It is sufficient 
to calculate the product of the transfer matrices, M*^^^\ for a sufficiently large system length, Lz, 
and then to calculate the eigenvalues, (a- If Lz is large enough, then the obtained values (a lie 
close to the mean values, {Ca)- This procedure is applied in the finite size scaling analysis of 
Lyapunov exponents, discussed in Section [T2.1l 

In what follows we order > (2 > ■ ■ ■ > Cn ■ That means that 

Ca = ZN+l~a, (348) 

where Za are discussed in Sect. El Clearly, the smallest Lyapunov exponent, zi, equals to the 
smallest eigenvalue Cn ■ 

In the numerical simulations, we need to calculate the smallest (in absolute value) parameter 
Cn = zi. This can be done by the following numerical algorithm. 
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1. Choose the required accuracy e of (n- Then, start with the following variables: first, we 
need the A*" x 2N matrix V which contains in N columns the mutually orthogonal vectors 
Va of the length 2N. We will also need the vectors da and Ca of the length N and put their 
stating values, da = (0, 0, . . . 0) and Co = (0, 0, . . . 0). Also, put = 0. 

2. Calculate the rij iterations 

V=Y[ M„V (349) 

3. Perform the Schmidt orthogonalization of vectors Ua, given by columns of the matrix U. 
Obtain new vectors Wa, which are already orthogonal to each other. 

4. Calculate the norm of vectors Wa, and add these norms into vectors d and e as follows: 

da = da + In \Wa\ (350) 



and 



5. Put 



ea = ea + {ln\wa\f. (351) 



Ua = Wa/\Wa\ (352) 

and the length, 

Lz = Lz + rii. (353) 
Calculate Ca = da/Lz, rja = ea/L^ and the accuracy 

VVa - Ca 



So 

6. If £jv > e then go to the step 2, otherwise stop. 



(354) 



The obtained value of Cat is the required smallest Lyapunov exponent zi. 

The length Lz, necessary for the calculation of the smallest Lyapunov exponent, depends on 
the required accuracy, s. For the 3D Anderson model, the length of the system, Lz was estimated 
inRef. [91]as 



Lz = ^L. (355) 



In numerical simulations, the accuracy e = 0.001 is usually required. Then, Lz ~ 500.000 x L. 
The number of iterations between two successive Schmidt orthogonalization is « 6 — 8 for 
the critical disorder, and smaller (larger) when W l^- Wc{W <C Wc), respectively. 
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E Calculation of the conductance 
E.l Id case 

It is convenient to write the Schrodinger Equation il5l in the form 

\ A* *n \ f E- En -1 \ f 



To calculate the conductance, consider the sample of the length sites^, connected to two 
semi-infinite ideal leads. This means that £„ ~ for both n <0 and for n > N ~ 1. 

Note that the transfer matrix M, defined in Eq. ( I356> does not have the structure of the 
transfer matrix T. Indeed, T connects the propagating waves on the left and right hand side of 
the sample, while M relates the wave functions in the the site representation. Both matrices are 
connected by the transformation. 



T„ = Q-^M„Q = ( ^ ^" ^ ) Q, (357) 



E-En -1 



where 

Q = ( l-^k ) • (358) 

Consider the electron coming from the right hand side of the system. Then, on the left hand 
side of the system, there is only a transmitted wave going to the left. Its wave function at the sites 
n = and n = — 1 can be written as 

The wave function on the right-hand side of the system is given by superposition of the incoming 
and reflected waves. We can use the transfer matrix, M, to express the wave function at sites 
and L, — 1: 



Now we multiply both sides of the of Eq. ( I360> by the matrix Q ^. We get 



Q^M ) -Q-^M^^QQ-M ]■ (361) 



We obtain that 



2ismk V -e"'**L, +*L,-i / 2ismk \ -e~'**o + * 



-1 



"For simplicity, we use the lattice constant, a as the unit length 
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where we used the relation J357> . between the transfer matrices M and T. Using the explicit 
form of and given by Eqs. ( I359> we find that the rh.s. of Eq. ( I362t equals to 



r(i.) 




1 

From Eq. ( 13621 we obtain 

1 / e*^*L, \ _f r+(i-)-i 



2isinfc V +*L,-i y V (*")"^ 

So we obtain the transmission coefficient in the form 



(363) 



(364) 



T=rr= ^ ^. (365) 



Equation ( I365I I is very useful for numerical calculations of the transmission T. 

E.2 Quasi-ld case 

To calculate the conductance in the quasi-ld case, we have to generalize the method of the 
previous Section. We again assume that our system is connected to two semi-infinite leads with 
zero disorder, W = 0. The electron is coming from the right and is scattered by the sample. The 
resulting waves either continue to the left in the left lead, or travel back to the right in the right 
side lead. 

In numerical simulations, we use the transfer matrix M, given by Eq. ( I343> 

M(^^) = nM„ = n(f-^" -J). (366) 

n=l n=l ^ ^ 

Similarly as in the case of the ID problem, discussed in the previous section, we have to transform 
this transfer matrix into the "wave" representation. In this representation, the transfer matrix in 
the leads is diagonal. Therefore, in the first step we have to diagonalize the transfer matrix 

Mo=(f~^° -J) (367) 

where Ho is Hamiltonian of the transversal slice without disorder 

In general, the transfer matrix has some eigenvalues with modulus equal to 1, i.e. A = 
exp ikz. The corresponding eigenvectors represent the propagating waves. Other eigenvalues are 
of the form A = exp ±k. They correspond to the evanescent modes. Note, if A is an eigenvalue, 
then A~^ is also an eigenvalue corresponding to the wave traveling in the opposite direction. 

As the transfer matrix Mq is not Hermitian, we have to calculate both the left and right eigen- 
vectors. Then, we construct four matrices: The N x 2N matrix i?icft (-Rright) which contains in 
its columns N right eigenvectors, for waves traveling to the left, (right), respectively. Similarly, 
matrices iioft and Ljight are 2N x N matrices which contains in N rows the left eigenvec- 
tors of the transfer matrix Mq which represent the waves traveling to the left and to the right, 
respectively. Then, by definition of the transfer matrix, 

rp _ / Til Ti2 \ f irightMi?right -LrightMi?ioft \ (368) 
\ T21 T22 J \ iloftMi?right LloftMi?ioft J 
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so that 



T22 = LicftTi?; 



left- 



(369) 



At this point we have to distinguish between the propagating and evanescent modes. If there are 
evanescent modes, then T22 7^ [^~] ~^ because the two matrices have different size. If we order 
the eigenvectors in the matrix R in such a way that the eigenvectors with index 1 < a < iVo cor- 
respond to the propagating modes, and the remaining eigenvectors correspond to the evanescent 
modes, then the transmission is given by the A^o x A^o sub-matrix [T22]a6, with a,b < Nq: 



T = 



ab=l 



22 lab 



(370) 



Other matrix elements of T22 correspond to the scattering of the electron into the evanescent 
channels. They do not contribute to the transmission since evanescent waves decay to zero in the 
semi-infinite leads. 

It seems that the relation ( I369l l solves our problem completely. However, the above algorithm 
must be modified. The reason is that the elements of the matrix (t^)^^ are given by their largest 
eigenvalues. We are, however, interested in the largest eigenvalues of the matrix . As the 
elements of the transfer matrix increase exponentially in the iteration procedure given by Eq. 
( I366t . any information about the smallest eigenvalues of will be quickly lost. We have 

therefore to introduce some re-normalization procedure. We use the procedure described in 
Ref. [72]. 

Relation ( I369l l can be written as 



T2: 

where we have defined the N x 2N matrices r^"^ rt = 0, 1, . . . as 

r(") =M„r("-i), and r'") = i?icft 
Each matrix r^") can be written as 

72 

with ri, r2 being the N x N matrices. We transform r as 

1 



r = r ri, 



(371) 



(372) 



(373) 



(374) 



and define r^"-* — M„(r')^" ^\ In contrast to the matrices ri and r2, all eigenvalues of the 
matrix r2r^^ are of the order of unity. The relation ( I37H can now be re-written into the form 



T2: 



loft (n) 



, -1 I '1 '1 



■ ' 1 '1 



(375) 



from which we get 



[T2; 



.(0) 



.(1) 



^Icft 



(376) 
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All elements of the matrices on the r.h.s. of Eq. (I376t are of the order of unity. 

References 

[1] P. W. Anderson: Phys. Rev. 109, 1492 (1958) 

[2] I. M. Lifshitz: Sov. Phys. Uspekhi 7, 549 (1965) 

[3] E. Abrahams et ah: Phys. Rev. Lett. 42, 673 (1979) 

[4] R. Landauer: IBM J. Res. Dev. 1, 223 (1957); R. Landauer: Phil. Mag. 21,683 (1970) 

[5] N. F. Mott, W. D. Twose: Adv. Phys. 10, 107 (1961) 

[6] A. MacKinnon, B. Ki-amer: Z. Phys. B 53, 1 (1983) 

[7] P A. Mello, P Pereyra, N. Kumar: Ann. Phys. (NY) 181, 290 (1988) 

[8] J. S. Langer, T. Neal: Phys. Rev. Lett. 16, 984 (1966) 

[9] P A. Lee, A. D. Stone: Phys. Rev. Lett. 55, 1622 (1985); P A. Lee, A. D. Stone, H. Fukuyama: 
Phys. Rev. B 35, 1039 (1987) 

[10] F. Wegner: Z. Phys. B 25, 327 (1976); ibid 35, 207 (1979) 

[11] S.Hikami: in [37] p. 429 

[12] B. L. Altshuler, V. E. Ki-avtsov, I. V. Lemer: Sov. Phys. JETP 64, 1352 (1986); JETP Lett. 43, 441 
(1986) 

[13] B. L. Altshuler: JETP Lett. 41, 648 (1985); B. L. Altshuler, D. E. Khmelnitskii: JETP Lett. 42, 

359 (1985); O. Tsyplyatyev etaJ.: cond-mat/0306613 

[14] S. Washburn, R. A. Webb: Adv. Phys. 35, 375 (1986) 

[15] A. B. Fowler, J. J. Wainer, R. A. Webb: IBM J. Res. Dev. 32, 372 (1988) 

[16] P W. Anderson, D. J. Thouless, E. Abrahams, D. S. Fisher: Phys. Rev. B 22, 3519 (1980) 

[17] A. MacKinnon, B. Ki-amer: Phys. Rev. Lett. 47, 1546 (1981) 

[18] D. J. Thouless: Phys. Rev. Lett. 39, 1167 (1977) 

[19] R. A. Serota, R. K. Kalia, P A. Lee: Phys. Rev. B 33, 8441 (1986) 

[20] P Markos, B. Kramer: Annalen der Physik 2, 339 (1993) 

[21] E. Abrahams, S. V. Kravchenko, M. R Sarachik: Rev. Mod. Phys. 73, 251 (2001) 

[22] C.M.SoukoulisetaJ.: Phys. Rev. Lett. 62,575 (1989) 

[23] M. C. W. van Rossum, Th. M. Nieuwenhuizen: Rev. Mod. Phys. 71, 313 (1999) 

[24] M.Stoychev, A. Z. Genack: Phys. Rev. Lett. 79, 309 (1997) 

[25] G. Bergmann: Phys. Rep. 107, 1 (1984) 

[26] D. J. Thouless: Phys. Rep. 13,93 (1974) 

[27] R A. Lee, T. V. Ramakrishnan: Rev. Mod. Phys. 57, 287 (1985) 

[28] M. JanBen: Phys. Rep. 295, 1 (1998) 

[29] M. JanBen: Int. J. Mod. Phys. B 8, 943 (1994) 

[30] B. Kramer, A. MacKinnon: Rep. Prog. Phys 56, 1469 (1993) 

[31] K. von Klitzing: Rev. Mod. Phys. 58, 519 (1986); B. Huckestein: Rev. Mod. Phys. 67, 357 (1995) 

[32] B. Kramer, T. Ohtsuki, S. Ketteman: Phys. Rep. 417, 211 (2005) 

[33] K. B. Efetov: Adv. Phys. 32, 53 (1983) 

[34] A. Mirlin: Phys. Rep. 326,259 (2000) 

[35] C. W. J. Beenakker: Rev Mod. Phys. 69, 731 (1997) 



Numerical Analysis of the Anderson Localization 



119 



[36] B. Kramer, G. Schon (Eds.) Anderson Transition and Mesoscopic Fluctuations Proc. of Workshop, 

Braunschweig 9-12 January 1990. Physica A 167 (1990) n. 1 
[37] B. Kramer (Ed.) Quantum Coherence in Mesoscopic Systems NATO ASI 254, Plenum Press NY 

and London (1991) 

[38] T. Brandes and S. Ketteman (Eds.) Tiie Anderson Transition and its Ramifications - Localisa- 
tion, Quantum Interference, and Interactions Springer Verlag "Lecture Notes in Physics" vol. 630, 
(2003) (ISBN 3-540-40785-5). Proceedings of the B. Kramer 60"^ Birthday Conference, Hamburg, 
September 4-6, 2002 

[39] E. N. Economu: Green's Functions in Quantum Physics 2nd ed. Springer (1990) 
[40] N. F. Mott, E. A. Davis: Electron processes in non-crystalline materials. Clarendon Press, Oxford 
1979 

[41] Y. Imry: Introduction to Mesoscopic Physics, Oxford Univ. Press (2000) 

[42] S. Datta: Electronic Transport in Mesoscopic Systems, Cambridge University Press, Cambridge, 
UK (1995) 

[43] Ping Sheng: Introduction to Wave Scattering, Localization, and Mesoscopic Phenomena, Aca- 
demic Press (2000) 
[44] M. L. Mehta: Random Matrices (Academic, NY 1991) 
[45] K. Slevin, Y. Asada, L. I. Deych: Phys. Rev. B 70, 054201 (2004) 
[46] J. M. Ziman: J. Phys. C 1, 1532 (1968); ibid 2, 1232, (1969); 2, 1704 (1969) 
[47] S. N. Evangelou, T. Ziman: J. Phys. C 20, L235 (1987) 
[48] T. Ando: Phys. Rev. B 40, 5325 (1989) 

[49] Y. Asada, K Slevin, T. Ohtsuki: Phys. Rev. Lett. 89, 256601 (2002); Phys. Rev. B 70, 035115 

(2004) 

[50] B. Huckestein, B. Kramer: Sohd State Comm. 71, 445 (1989); Phys. Rev. Lett. 64, 1437 (1990); 

B. Huckestein: in [37] p. 441. 
[51] J. T. Chalker, P D. Coddington: J. Phys. C 21, 2665 (1998) 
[52] B. Bulka, B. Kramer, A. MacKinnon: Z. Phys. B: Condens. Matt. 60, 13 (1985) 
[53] N. F. Mott: Philos. Mag. 26, 1015 (1972) 
[54] P Marko§, L. Schweitzer: J. Phys. A 39, 3221 (2006) 
[55] S. N. Evangelou: J. Phys. A 23, L317 (1990) 
[56] J. Brndiar, P Markos: cond-mat/0606056 
[57] P. Markos, L. Schweitzer, unpublished 
[58] B. I. Shklovskii et al: Phys. Rev. B 47, 11487 (1993) 

[59] A. G. Aronov, V. E. Kravtsov, 1. V. Lemer: JETP Lett. 59, 40 (1994); Phys. Rev. Let. 74 1174 
(1995) 

[60] D. Braun, G. Montambaux, M. Pascaud: Phys. Rev. Lett. 81, 1062 (1998) 
[61] 1. Kh. Zarekeshev, B. Kramer: Phys. Rev. Lett. 79, 717 (1997) 
[62] J. T. Edwards, D. J. Thouless: J. Phys. C 5, 807 (1972) 
[63] D. C. Licciardello, D. J. Thouless: J. Phys. C 8, 4157 (1975) 

[64] E. N. Economou, C. M. Soukoulis: Phys. Rev. Lett. 46, 618 (1981); ibid 47 973 (1981) 

[65] M. Azbel: Phys. Lett. A 78, 410 (1980) 

[66] D. S Fischer, P A. Lee: Phys. Rev. B 23, 6851 (1981) 

[67] D. C. Langreth, E. Abrahams: Phys. Rev. B 24, 2978 (1981) 



120 P. Markos 

[68] R. Landauer: Z. Phys. B 68,217 (1987) 

[69] D. Braun et aJ.: Phys. Rev. B 55, 7557 (1997) 

[70] A. Abrikosov, I. A. Rhyzkin: Adv. Phys. 27, 147 (1978); V. I. Melnikov: Fiz. Tverd. Tela 23, 444 

(1981); A. Abrikosov: Solid State Commun. 37, 997 (1981) 

[71] R D. Kirkman, J. B. Pendry: J. Phys. C: Solid St. Phys. 17, 4327 (1984); ibid 5705 (1984) 

[72] J. B. Pendry, A. MacKinnon, P. J. Roberts: Proc. R. Soc. London A 437, 67 (1992) 

[73] O.N.Dorokhov: JETPLett. 36,318 (1982) 

[74] P A. Mello: Phys. Rev. B 35, 1082 (1987) 

[75] B. Shapiro: Philos. Mag. B 56, 1032 (1987) 

[76] R Vagner et a/.: Phys. Rev. B 67, 165316 (2003) 

[77] P. Markos, C. M. Soukoulis: Wave propagation, Princeton Univ. Press (submitted) 

[78] M. Mosko etaJ.: Phys. Rev. Lett. 91, 136803 (2003) 

[79] D. Maily, M. Sanquer: in [37] p. 401 

[80] R Marko§, B. Kramer: Solid State Comm. 90, 615-617 (1994) 

[81] M. Kaveh, N. F. Mott: J. Phys. C 14, L177 (1981) 

[82] Y. Asada, K. Slevin, T. Ohtsuki: cond-mat/0509718 

[83] R A. Mello, A. D. Stone: Phys. Rev. B 44, 3559 (1991) 

[84] E. N. Economou. C. M. Soukoulis, A. D. Zdetsis: Phys. Rev. B 30, 1686 (1984) ibid. 31, 6483 
(1985) 

[85] C. W. J. Beenakker: Phys. Rev. B 47, 15763 (1993) 

[86] D. Braun etaJ.: Phys. Rev. B 64, 155107 (2001) 

[87] M. Riihlander, P Markos, C. M. Soukoulis: Phys. Rev. B 64, 172202 (2001) 

[88] L Travenec: Phys. Rev. B 69,2004 (033194) 

[89] M. C. W. van Rossum et al.: Phys. Rev. B 55, 4710 (1997) 

[90] K. A. Muttalib: Phys. Rev. Lett 65, 745 (1990) 

[91] R Markos: J. Phys.: Condens. Matter 7, 8361 (1995) 

[92] Y. Imry: Europhys. Lett. 1,249(1986) 

[93] P Markos, C. M. Soukoulis: Phys. Rev. B 71, 054201 (2005) 

[94] Th. M. Nieuwenhuizen, M. C. W. van Rossum: Phys. Rev. Lett. 74, 2674 (1995) 

[95] E. Kogan, M. Kaveh: Phys. Rev. B 52, R3813 (1995) 

[96] K. Tankei, M. Ikegami, Y. Nagaoka: J. Phys. Soc. Jpn. 63, 1090 (1994) 

[97] R Wegner: Nucl. Phys. B 316, 623 (1989) 

[98] S. Hikami: Phys. Rev. B 24, 2671 (1981) 

[99] U. Fastenrathetai.:PhysicaA 172,302(1991) 

[100] L. Schweitzer, I. Kh. Zarekeshev: J. Phys.: Condens. Matt. 9, L441 (1997) 

[101] B. Shapiro: Phys. Rev. Lett. 65, 1510 (1990) 

[102] R Markos, B. Kramer: Philos. Mag. B 68, 1993 (357-379) 

[103] A. Cohen, Y. Roth, B. Shapiro: Phys. Rev. B 38, 12125 (1988) 

[104] A. Cohen, B. Shapiro: Int. J. Mod. Phys. B 6, 1243 (1992) 

[105] R Markos: Phys. Rev. Lett. 83, 1999 (588) 

[106] M. Riihlander, R Markos, C. M. Soukoulis: Phys. Rev. B 64, 212202 (2001) 

[107] K. Slevin, T. Ohtsuki: Phys. Rev. Lett. 78, 4083 (1997) 



Numerical Analysis of the Anderson Localization 



121 



[108] P. Markos: in [38] 

[109] V. A. Gopar, K. A. Muttalib, P. Wolfle: Phys. Rev. B 66, 174204 (2002) 

[110] K. A. Muttalib et a/.: Europhys. Lett. 61, 95 (2003) 

[111] P Markos: Europhys. Lett. 26, 431 (1994) 

[112] K. Slevin, T. Ohtsuki: Phys. Rev. B 63,45108 (2001) 

[113] L Travenec, R Markos: Phys. Rev. B 65,114109(2002) 

[114] C. M. Soukoulis etal.: Phys. Rev. Lett. 80, 668 (1999); K. Slevin, T. Ohtsuki: Phys. Rev. Lett. 80, 

669 (1999) 

[115] K. Slevin, T. Ohtsuki, T. Kawarabayashi: Phys. Rev. Lett. 84, 3915 (2000) 
[116] M. Schreiber, H. GruBbach: Phys. Rev. Lett. 76, 1687 (1996) 

[117] R. Rammal, G. Thoulouse: J. Phys. Lett. (France) 44, 1983 (L13); R. Burioni, D. Cassi, S. Regina: 

cond-mat/97 10273 
[118] R Markos: Phys. Rev. B 65, 104207 (2002) 
[119] P J. Roberts: J. Phys.: Condens. Matt. 4,7795 (1992) 
[120] K. M. Slevin, J. B. Pendry: J. Phys.: Condens. Matt. 2, 2821 (1990) 
[121] J. Prior, A. M. Somoza, M. Ortuno: Phys. Rev. B 72, 024206 (2005) 
[122] A. M. Somoza, J. Prior, M. Ortuno: Phys. Rev. B 73, 184201 (2006) 
[123] L W. Kantelhardt, A. Bunde: Phys. Rev. B 66,035118 (2002) 
[124] S. L. A. de Queiroz: Phys. Rev. B 66, 195113 (2002) 

[125] K. A. Muttalib, R Wolfle: Phys. Rev. Lett. 83, 3013 (1999); R Wolfle, K. A. Muttalib: Ann. Phys. 

(Leipzig) 8, 753 (1999) 

[126] L. S. Froufe-Prez et al.: Phys. Rev. Lett. 89, 246403 (2002) A. Garcia-Martin, J. J. Saenz: Waves 

in Random and Complex Media 15, 229 (2005) and references therein 
[127] J.-L. Pichard, G. Sarma: J. Phys. C 14, L127 (1981); ibid, L167 
[128] R Markos: J. Phys. A: Math. Gen. 33, L393 (2000) 
[129] S. N. Evangelou: private communication 

[130] B. Kramer, M. Schreiber, A. MacKinnon: Z. Phys. B 56, 297 (1984); M. Schreiber, B. Kramer, 
A. MacKinnon: J. Pys. C: Solid State Phys. 17, 4111 (1984); B. Kramer, K. Broderix, A. MacKin- 
non, M. Schreiber: Physica A 167, 163 (1990) 

[131] P Cain, R. A. Romer, M. Schreiber: Ann. Phys. (Leipzig) 8, SI-33 (1999) 

[132] R. A. Romer, M. Schreiber: in [38] p. 3. 

[133] K. M. Slevin, T. Ohtsuki: Phys. Rev. Lett. 82, 382 (1999) 

[134] A. MacKinon, B. Kramer: Phys. Rev. Lett. 49, 695 (1982) 

[135] A. MacKinnon: J. Phys.: Condens. Matt. 6, 2511 (1994) 

[136] K. Slevin, R Markos, T. Ohtsuki: Phys. Rev. Lett. 86, 3594 (2001) 

[137] K. Slevin, R Markos, T. Ohtsuki: Phys. Rev. B 67, 155106 (2003) 

[138] B. L. Altshuler et ai.: Sov. Phys. JETP 67, 625 (1988) 

[139] E. N. Evangelou: Phys. Rev. B 39, 12895 (1989); ibid. 49, 16805 (1994) 

[140] M. Batsch et al.: Phys. Rev. Lett. 77, 1552 (1996) 

[141] A. Kaneko, T. Ohtsuki: Ann. Phys. (Leipzig) 9, 121 (1999) 

[142] K. Slevin, Y. Asada, T. Ohtsuki: Phys. Rev. B 73, (2006) 

[143] F. Evers, A. D. Mirhn: Phys. Rev.Lett. 84, 3690 (2000); A. D. Mirlin, F. Evers: Phys. Rev. B 62, 
7920 (2000); A. Mildenberger, R Evers, A. D. Mirlin: Phys. Rev. B 66, 033109 (2002) 



122 P. Markos 

[144] F. Evers, A. Mildenberger, A. D. Mirlin: Phys. Rev. B 64, 241303(R) (2001) 

[145] E. Cuevas: Phys. Rev. B 66,233103 (2002) 

[146] E. Cuevas et al.: Phys. Rev. Lett. 88, 016401 (2002) 

[147] D. VoUhardt, P. Wolfle: Self-Consistent Theory of Anderson Localization, in Electronic Phase 

Transitions, Ed. by W. Haake and Yu. V. Kopaev, Elsevier Sci. Publ. 1992 1992 

[148] R Marko§, M. Henneke: J. Phys.: Condens. Matter 6, L765 (1994) 

[149] I. M. Suslov: cond-mat/0 105325 2001; cond-mat/0106357 (2001) 

[150] 1. Suslov: J. Exp. Theor. Phys. 101,661 (2005) 

[151] V.N. Kuzovkov, V. Kashcheyevs, W. von Niessen: J. Phys.: Cond. Matt. 69, 1683 (2004) 

[152] R Markos, L. Schweitzer, M. Weyrauch: J. Phys.: Condens. Matt. 16, 2004 (1679-1681) 

[153] I. Zambetaki ef al: Phys. Rev. Lett. 76, 3614 (1996); Phys. Rev. B 56, 12221 (1997) 

[154] A. Croy, R. A. Roerner, M. Schreiber: in Lectiu'e Notes in Coinputationai Science and Engineering 

(Springer, Berlin), (cond-mat/0602300) 

[155] V. Dobrosavljevic, A. A. Pastor, B. K. Nikolic: Europhys. Lett. 62, 76 (2003) 

[156] A. MacKinnon: Z. Phys. B 59,385 (1985) 

[157] L. Schweitzer: in [38] 

[158] V. I. Fal'ko, K. B. Efetov: Europhys. Lett. 32, 627 (1995); Phys. Rev. B bf 52, 17413 (1995) 

[159] L. Schweitzer, P Markos: Phys. Rev. Lett. 95, 256805 (2005) 

[160] A. Hansen et al. Annual Rev. of Comp. Phys. 5, Ed. by D. Stauffer (World Scientific, Singapore, 

1997 (cond-mad/9703026) 

[161] S. Kivelson, D. H. Lee, S. C. Zhang: Phys. Rev. B 46, 2223 (1992) 

[162] X. Wang, Q. Li, C. M. Soukoulis: Phys. Rev. B 58,3576(1998) 

[163] Z. Wang, B. Jovanovic, D. H. Lee: Phys. Rev. Lett. 77, 4426 (1996) 

[164] T. Ohtsuki, K Slevin, B. Kramer: Physica E 22, 248 (2004) 

[165] K. A. Muttalib, J. R. Klauder: Phys. Rev. Lett. 82, 4272 (1999); K. A. Muttalib, V. A. Gopar: Phys. 
Rev.B 66, 115318(2002) 

[166] R Marko§ et al.: Europhys. Lett. 68, 867 (2004) 

[167] K. A. Muttalib, P Markos, P Wolfle: Phys. Rev B 72, 125317 (2005) 

[168] R Markos: Ann. Physik (Leipzig) 8,81 165 (1999) 

[169] R Markos: J. Phys. A: Math. Gen. 30, 3441 (1997) 

[170] L. I. Deych, A. A. Lisyansky, B. I. Altshuler: Phys. Rev. Lett 84, 2678 (2000); Phys. Rev. B 64, 
224202 (2001); L. I. Deych, M. V. Erementchouk, A. A. Lisyansky: Phys. Rev. Lett. 90, 126601 
(2003); Physica B 338, 79 (2003) 

[171] F. M. Izrailev, A. A. Krokhin: Phys. Rev. Lett. 82, 4062 (1999); F. M. Izrailev, A. A. Krokhin, 
S. E. Ulloa: Phys. Rev. B 63, 041 102 (2001) 

[172] F. M. Izrailev, N. M. Makarov: J. Phys. A 49, 10613 (2005) 

[173] Ch. Mudry, P W. Brouwer, A. Furusaki: Phys. Rev. B 62, 8249 (2000); P. M. Brouwer: et al.: 

cond-mat/0511622 (2005) and references therein 
[174] A. A. Asatryan et al.: Phys. Rev. E 71, 036623 (2005) 
[175] E. Larose et al.: Phys. Rev. Lett. 93, 048501 (2004) 

[176] R Cain etaJ.: cond-mat/0 106005 (2001); K. Slevin, T. Ohtsuki: cond-mat/0 106006 (2001) 
[177] K. Slevin, T. Nagao: Int. J. Mod. Phys. 21,2665 (1988) 
[178] R A. Mello, J.-L. Pichard: J. Phys. I 1, 493 (1991) 



Numerical Analysis of the Anderson Localization 



123 



[179] J.-L. Pichard: in: [37] p. 369. 

[180] P. A. Mello, S. Tomsovic: Phys. Rev. B 46, 15963 (1992) 

[181] J.-L. Pichard et a/.: J. Phys. France 51, 587 (1990) 

[182] N. Zannon, J.-L. Pichard: J. Phys. France 49, 907 (1988) 

[183] K. M. Slevin, J.-L. Pichard, K. A. Muttalib: J. Phys. France 1 3, 1387 (1993) 

[184] Y. Avishai, J.-L. Pichard, K. A. Muttalib: J. Phys. France 3, 1387 (1991) 

[185] T. A. Brody etal.: Rev. Mod. Phys. 53, 385 (1981) 

[186] C. W. J. Beenakker, R. Rajaei: Phys. Rev. Lett. 71, 3689 (1993); Phys. Rev. B 49, 7499 (1994) 

[187] D. J. Thouless: J. Phys. C 5, 77 (1972) 

[188] M. Kappus, F Wegner: Z. Phys. B 45, 15 (1981) 

[189] B. Derrida, E. Gardner: J. Physique 45, 1283 (1984) 

[190] V. I. Oseledec: Trans. Moscow Math. Soc. 19, 197 (1968) 



